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ABSTRACT 

Several multivariable frequency domain techniques that 
can be used to assist with the design of linear, feedback control 
systems are investigated and compared with other methods. The inverse 
Nyquist array, the characteristic locus, the direct Nyquist array, the 
optimal quadratic and the multiloop methods are evaluated, by 
experimental application to a double-effect, pilot plant evaporator. 

The direct Nyquist array method was derived by direct, 
intuitive extension of the conventional single-input, single-output, 
Nyquist design procedure and also based on the same theory that 
underlies the inverse Nyquist array and characteristic locus methods. 
It was concluded that the direct Nyquist array technique is a more 
practical and convenient design technique than the inverse Nyquist 
array and the characteristic locus methods and the intuitive development 
provides greater insight into the physical meaning of each step in the 
design procedure. The direct Nyquist array method was also applied to 
to the design of a control system for a binary distillation column 
model which contained pure time delays. 

The characteristic locus method was found to be slightly 
more general, but less practical, than either the direct or the inverse 
Nyquist array methods. 

The three frequency domain design techniques were also 
compared with a multiloop control scheme and the optimal regulator 
control method. All of the multivariable techniques were found to 
produce better control of the evaporator than the multiloop contro] 
scheme. The frequency domain techniques produce controllers that gave 


closed-loop performance equivalent to other multivariable schemes but 
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gave the designer better insight and opportunity for incorporating 
practical factors such as integrity, into the design, than optimal 


control techniques. 
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Schematic Diagram of the Double Pilot Plant 
Evaporator used for this work. 
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CHAPTER ONE 


INTRODUCTION 

During the last two decades considerable attention has been 
given to the design of multivariable control systems. Most of this 
attention has been focused on the development of new techniques which 
are based on the concepts of modern control theory and use a state- 
variable formulation. Unfortunately many of the state-variable methods 
have been found [16,26] to be of limited practical value in the design 
of multivariable control systems for chemical processes. The main 
reason for the limited success with these methods is due to the fact that 
modern control theory has arisen from the study of systems with 
Significantly different characteristics than chemical processes. 
Because of this, these methods solve only a small portion of the total 
design requirements of these complex processes and the design usually 
requires measurements (or estimates) of the state variables which often 
are not available because of physical or economical constraints. 

Furthermore, controllers designed using methods based on 
linear theory, but applied to systems that are inherently non-linear, 
usually require tuning when they are implemented in the field. The 
state-variable methods do not offer any guidance to the designer in the 
way the controller must be altered if an unsatisfactory solution is 
obtained. With most of these methods, the only alternative if a 
controller is found to be unsatisfactory is to redesign and retest it in 
the field. 

On the other hand single-variable frequency-domain techniques 
have been found to be very effective in designing controllers for single- 
input, single-output systems. These techniques are applicable to systems 
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that contain time delays which are commonly found in chemical processes 
and they provide the designer with information such as the gain margin, 
Phase margin, crossover frequency, etc. which are very useful for the 
tuning of the controller in the field. 

In spite of the success of single-variable frequency domain 
techniques, the extension of these techniques to multivariable control 
has received relatively little attention. Only in the last decade has 
this trend been redirected. It began with Rosenbrock [65] who developed 
the inverse Nyquist array method in the 1960's. Since then several 
techniques, which are essentially generalizations of the stability 
criteria introduced by Bode [7] and Nyquist [54] for single-variable 
systems, have been proposed [4,36,44,58]. Because of their recent 
development the experience with most of these methods is limited and, 
with the exception of the inverse Nyquist array method, experimental 
applications based on these metnods haye not been reported in the 
literature. Since most of these methods provide information similar to 
the single-variable techniques and produce a controller which is easily 
tuned in the field, it is foreseeable that they will play as important 
a role in the design of multivariable control systems as the single- 
variable frequency domain techniques did with single-variable processes. 
However there is still a need for additional study of these methods in 
order to improve them or to develop more practical methods which wil] 
make the design of control systems even easier. There is also a need 
for more applications of these methods to practical systems in order to 


encourage their use in industrial processes. 


' 


Ty 


' | 1 : 
: P ais ~ ae 5” ot : teal t 
Ft iA) oeie: . Tg * Ore nes F i 
is 4 : ® J ; 


= ™ 


1.1 Objectives of This Study 


Based on these considerations the study of the multivariable 
frequency domain design techniques was selected as the objective of this 
investigation. Because of the vastness of this field this study was 
limited, after a general review of this area, to a detailed evaluation of 
the methods which were considered to be the most promising techniques. 
This study has evaluated the inverse Nyquist array [65] and the character- 
istic locus method [4,42]. It has also developed and evaluated the 
"direct" Nyquist array method. Although the "direct" technique has been 
suggested [38,68,70,71] in the literature it has not been fully developed or 
applied. The direct Nyquist array method was derived, in this study, by 
two different approaches. First by modification and extension of the 
Same theory that underlies the inverse Nyquist array and cnaracteristic 
locus methods and secondly by extension of tne single-variable frequency 
domain techniques using an intuitive approach. The evaluation of these 
methods has been concerned with several different features of the method: 
the theoretical aspects, their computer implementation and the application 
of these techniques to the design of a control] system for the double- 
effect evaporator pilot plant located in the Department of Chemical 
Engineering at the University of Alberta. 

Although it is not included as part of this thesis this study 
has also included the development of a computer-aided design package 

[29,30,31,32]. This program assists with the design of controllers 
based on the inverse Nyquist array, characteristic locus and direct 
Nyquist array methods. The package has been written in a user 
oriented form and the controller can be designed in an interactive way 


using a CRT computer terminal. This package has been incorporated as a 
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part of the GEMSCOPE [14] package which operates on an IBM 360 series 
computer under the University of Michigan time-sharing system (MTS). 
This investigation has also covered the comparison of these 
techniques with the optimal control method and the multiloop design 
approach based on the experimental results obtained by applying these 


methods to the double-effect evaporator pilot plant. 


jecepotructupe of thealhesis 

This thesis is divided in eight chapters. Chapter Two contains 
a survey of the frequency-domain design procedures which have been 
proposed to date. The inyerse Nyquist array method js discussed jn 
Chapter Three and two design examples are included which illustrate the 
use of this method. Chapter Four deals with the characteristic locus 
method and presents a design example using this technique. 

In Chapter Five the attention is focused on the erereorent 
and application of the direct Nyquist array method. This chapter is 
divided in three main parts. In the first part the direct Nyquist 
array method is derived based on concepts which were introduced in the 
development of the inverse Nyquist array [65] and the characteristic 
locus method [4,42]. In the second part the direct Nyquist array 
method is shown to be a direct and intuitiye extension of single-input 
Single-output design techniques. This is a self contained section and 
jt is based on a report written by Fisher and Kuon [15]. In the last 
part the direct Nyquist array method is applied to design a contro] 
system for the double effect evaporator pilot plant, 

The experimental eyaluation of the controllers designed using 


frequency domain techniques is presented in Chapter Six and includes a 
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comparison with controllers designed using the optimal control method 
and the multiloop design approach. 

In Chapter Seven the direct Nyquist array method is used to 
design a control system for a multivariable process which contains time 
delays. A model [79] which represents the distillation column located 
at the Department of Chemical Engineering of the University of Alberta 
has been used in this example. Finally, in Chapter Eight, the overall 


conclusions of this work are summarized. 


1.3 Comments on Thesis Format and Use 

The following comments are included to assist the reader in 
identifying those parts of the thesis that are most directly applicable 
to his particular interests and needs: 
1) The description of the inverse Nyquist array, characteristic loci, 
and direct Nyquist array methods in Chapters Three, Four and Five respectively, 
were written so that they could be read, or referred to, either sequentially 
or independently. This required some redundancy in discussing certain 
basic concepts and was the reason why material, such as the description 
of the evaporator, was placed in appendices. Specific comments were 
included, where appropriate, to point out the similarities and/or 
relationships eR the different methods but these are not normally 
necessary to an understanding of any one method. | 
2) All the frequency domain techniques are interactive and require 
decisions and input at each step in the design procedure. These decisions 
are often subjective and each user would probably make different decisions 
and follow a slightly different path in designing a controller for a given 


application. Chapters Three, Four and Five contain rather detailed presen- 
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tations of the steps that the author took and comments about why these 
steps seemed appropriate. Each example is presented in a similar format 
to facilitate comparison of the "mechanics" of each method. If the 
computer aided design system is available locally some readers may prefer 
to skip over the details of each example and evaluate the methods by 
doing the examples themselves. The description of the examples is 
illustrative and not essential to the presentation of the theory and/or 
methods. 

3) The main conclusion of this work is that the direct Nyquist array 
method is the preferred design technique and should be used in place 

of the other methods. This conclusion is subjective and application 


dependent, but readers willing to accept it may prefer to consider 


Chapters One through Four as primarily historical and theoretical 


background to Chapter Five. 


4) Readers familiar with conventional single variable frequency domain 
techniques might prefer to read Section 5.3 before considering the details 
contained in the other chapters. (Section 5.3 is an intuitive explanation 
of the extension of single variable Nyquist techniques to the direct, 
multivariable Nyquist array approach. ) 

5) In order to handle practical problems, multivariable frequency domain 
design techniques require the use of interactive computer facilities to 

do the numerical calculations and plotting. This thesis is concerned with 
the design methods and the results of applying them to specific problems. 
The description of the computer programs, and instructions about how to use 
them, are contained in separate reports. (They are written so the user 
can apply them to practical problems without being familiar with all the 


theory and the "whys" discussed in this thesis.) 


6) All of the figures discussed in Chapters Two through Seven are located 

in the section of the thesis immediately preceeding the bibliography. This 
was done to facilitate multiple references to the same figure and comparisons 
between figures. The figures showing Nyquist diagrams etc. are direct, 


electronic copies of the display on the CRT of the computer terminal. 
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CHAPTER TWO 
LITERATURE SURVEY OF MULTIVARIABLE FREQUENCY DOMAIN TECHNIQUES 


2.1 INTRODUCTION 

Frequency domain techniques have proven to be efficient and 
effective design techniques for single variable feedback control problems 
Since they were first introduced by Bode [7] and Nyquist [54]. However 
it was not until the 1960's that any significant progress was made in 
the application of frequency domain techniques to multivariable control 
problems. A very important contribution was the development of the 
inverse Nyquist array method by Rosenbrock [65]. Other important 
techniques for the design of multivariable control systems using 
frequency domain techniques include: the commutative controller [33,36], 
the characteristic locus method [4,5,41,42], and the sequential return- 
difference method [45,46,47]. Several modifications and extensions of 
these techniques have also been proposed in the literature [e.g. 57,58]. 

Development of frequency domain techniques for the design of 
multivariable control systems has proceeded along two main paths. First, 
a number of investigators have tried to modify the problem formation, 
or develop a design approach, to permit the use of classical single 
variable design techniques. Secondly, others have assumed a true multi- 
variable formulation and developed new design techniques. However, 
almost all of the techniques can be interpreted as extensions and 
generalizations of the original work done by Bode and Nyquist. 

The following sections contain a review of frequency domain 
techniques presented in an order intended to emphasize the progressive 


development and inter-relation of the various techniques. This survey 
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is primarily concerned with the multivariable regulator problem and is 
restricted to techniques that can be used to assist with the design of 
linear feedback control systems which can be represented as shown in 
Figure 2.1. Unless it is otherwise stated the feedback transducer matrix 
will be considered to be an identity matrix since there is no loss of 
generality as far as the problem of determining system stability is 
concerned. All of the developments and equations are presented in terms 
of discrete systems but the analysis is also valid for continuous 
systems and most of them can be applied to systems represented by both 
Square, and non-square, transfer function matrices. The theory of some 
of the approaches can also handle the inclusion of a post compensator 
matrix, L(z), (which comes immediately after G(z) in Figure 2.1) but 
this is not usually of practical importance because it implies control 


1s exercised over a combination of the real plant output variables. 


2-2 MULTILOOP CONTROL SYSTEM DESIGN 

A common practice in the design of controllers for multivari- 
able systems is to neglect the interactions that occur between the 
system input and output variables and to proceed on the basis of a 
number of single variable control loops each of which can be designed 
independent of all bi others. This approach has been widely used and 
is frequently referred to as multiloop control. The main place that the 
multivariable aspects of the problem are considered is in the pairing of 
the system input and output variables and/or in the field tuning of the 
loops. A technique to assist with this problem of pairing system 
variables has been proposed [e.g. 8,9] and in most cases it can be 


interpreted as a simplification or special case of the more general 
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multivariable techniques discussed later. The actual controller design 
step is completed using any applicable design technique such as the 
classical Bode or Nyquist frequency domain methods. 

However in many systems the interactions between process 
variables are too severe to be neglected entirely. In such cases the 
application of conventional single variable techniques results in 
unsatisfactory control of the overall system. In some cases control 
engineers have intuitively designed controliers based on their previous 
experience followed by a tuning of all control loops in the field. 

The main reason for this lack of success is interaction j.e. 
the way in which a reference input r.(z) applied to input i affects the 
set of outputs (y;(z), j#i) [42]. Because of the interaction in a 
multivariable system the transfer function between the input r(z) and 
the output y;(z) will change when the other loops are closed and also 
when there are changes in the controller constants in the other loops. 
Thus in order for a single-loop technique to produce a suitable control- 
ler for loop i, the design must be based on the system transfer function 
of loop i that applies when loop 7 is open and the rest are closed. 

From a design point of view this is not very convenient since the 
transfer function involved in the design of loop i is a function of the 
controller designed and parameters in all other loops of the system. 
Thus, in the general case there is no obvious starting point and several 
assumptions must be made in order to proceed with the design of the 


first loop. 


Sequential Return-difference Method 


The Sequential Return-difference method developed by Mayne 
[44,45,46] in its most elementary form represents an efficient algebraic 
implementation of the approach that an experienced control engineer 
would probably take in applying single variable techniques to multi- 
variable problems. In this approach a controller k, (2) js designed for 
the first loop based on the transfer function 947 (2) between the input 
ry (Z) and the output y,(z) using a single-loop frequency domain 
technique. Then a new process transfer function matrix g!(z) is calcu- 
lated with this first loop closed and the rest open. A controller 
ky (z) for the second loop is then designed based on the transfer function 
dues (2) (element of g! (z) between input ry (z) and output yo (z)) using 
the same single loop design technique. A new process transfer function 
matrix 62 (z) is then calculated, this time with the first and second 
loops closed and the rest open. The procedure continues in a similar 
way for the rest of the loops. A schematic representation is shown in 
Figure 2.2. Mayne [44,45] has shown that a sufficient condition for 
stability of the overall multivariable system is that each of the indivi- 
dual control loops developed in this approach be stable. 

The algebraic updating of the transfer function is done in the 
following way: 

Let the controller matrix with loop i closed be K! = diag 
(k, (Zz), k,(z); we k;(z), 0... 0, 0) and assume that the first i 
controllers have already been obtained. The original process transfer 


function matrix G(z) is defined by: 


Viz) "eG (z)eutz) i (2.1) 


1] 


By closing the first i loops, and using negative feedback, u(z) is given 


by 


u(z) = -K'(z)y(z) + r'(z) t2e2) 
where r=K' Er, (2) 75 (2) iAP r(2)]. The new input-output relationship 


becomes 


y(z) = (1, + G(z)K"(z))7! @(z) r'(z) (2.3) 
or 


y(z) = F.7! G(z) r'(z) = @(z) r'(z) (2.4) 


The transfer function matrix @! (z) can also be easily 


calculated using one of Mayne's relationships [45,46] j.e. 


i a pt i-1_T> ,i-] 
Gazi) a= De k. (z) a e,] G@ ‘(z) (225) 
where: 
k.(z) = k,(z)/f, (2) (2.6) 
” i-l 

f,(z) =a let k; (z) 95 (z) amr) 
Tal a2 i-] 

g; = i column of G  (z) 

e! = vector whose ith element is unity and the 


remaining elements are zero. 
The selection of the controller for loop i is always based on 
the transfer function gii'(z) which is the (i,i1) element of aeay 
This is a very straightforward technique but it does not solve the 
basic problems associated with the design of multivariable systems. 


Specifically it does not provide a means of determining the best pairing 
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of input-output variables, nor does it give directions about the order in 
which the individual loops should be treated in the design procedure. 
Obviously the final controller will be influenced by. both decisions. 

More generally the technique does not provide a means for designing 
individual control loops (or picking individual controller constants) 
such that a performance criterion based on the overall multivariable 
system will be optimized. Thus while it might provide a convenient means 
of dealing with simple multivariable problems there is obviously a need 


for more rigorous multivariable techniques. 


2.3 MULTIVARIABLE CONTROL SYSTEM DESIGN 

Another approach to the application of single variable design 
techniques to multivariable problems has been to augment or modify the 
system by a pre- or post-compensator. If the compensator reduces the 
interactions between the input and output variables of the augmented 
plant then it will be possible to complete the actual controller design 
step using classical single variable techniques. The controller for the 
multivariable system then consists essentially of two parts: the multi- 
variable compensator and a set of single variable controllers. The 
concept of non-interaction is very appealing in some applications, but 
unfortunately in most problems it cannot be shown that non-interaction 
is a necessary condition for the design of an optimal controller. In 
fact in some cases it can be shown that a non-interacting controller will 
give poorer performance than other alternatives. [64] Thus in general terms 
one is trading off a guarantee of overall optimality for design conven- 


jence. 
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2.3.1 Non-interacting Design Techniques 


In this scheme the design of the controller is done in two 
steps. In the first step a pre-compensator K(z) is chosen so that the 
open loop transfer function matrix of the augmented plant, G(z)K(z) is 
diagonal. Since interaction has been eliminated by the compensator K(z), 
Single-loop theory may be used to design each loop separately. 

This design technique has not been widely used, especially on 

complex processes, because it often produces a very complicated 
controller which is expensive or impractical to implement. 

Also the selection of the pre-compensator matrix K(z) to make 

the system diagonal restricts the form of cumpensation that can 

be applied to each loop using single-loop theory. When the deter- 
minant of G(z) has zeros outside the unit circle this procedure gives 

a poor or unstable contrel [64]. This is easily shown by the following 
analysis. In order to obtain a diagonal open loop transfer function 
matrix for the system G(z), the pre-compensator matrix K(z) has to have 


the following form: 
K(z) = G7'(z) diag (4,(z), 4,(z) -.. 4,(2) (2.8) 


where 
Gy> Io» 43 = diagonal elements of the augmented plant G(z)K(z) 
The controller will then be unstable if the plant has right half plane 
ZeV0S. 
Non-interactive design techniques have been discussed by Tsien 


[75], Rae [62] and more recently by MacFarlane [40] and Berry [6]. 


2.3.2 Inverse Nyquist Array Technique 

The non-interacting design technique described above was a 
rigorous mathematical approach which could only be applied if the system 
met some specific mathematical condition so that the interaction could 
be eliminated completely. One would expect that if the objective was 
simply to reduce interactions rather than to eliminate them entirely 
then the "non-interacting" design approach could be applied to a wider 
range of problems and would involve compensators that were more practical 
to implement. The inverse Nyquist array method developed by Rosenbrock 
[65,67] essentially gives the control engineer a means of reducing the 
system interactions to a point where a controller design and overall] 
Stability analysis can proceed on the basis of classical single variable 
techniques. 

The inverse Nyquist array [65,67] is based on an extension of 
the Nyquist stability criterion for multivariable systems. The stabili- 
ty criterion is developed from the basic equation presented by McMorran 


[47], but which is equivalent to the relationship obtained earlier by 


Popov [61] 
in7'(z)| _ cLep a) 
ey | 
where: 


R(2Z)a=s Ces Q(z))"! Q(z) = closed-loop transfer function matrix 


(CLTFM) 
Q(z) = G(z)K(z) = open-loop transfer function matrix (OLTFM) 
CLCP = closed-loop characteristic polynomial 
OLCP = open-loop characteristic polynomial 
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The inverse Nyquist array works with the inverses of the open- 
loop and closed-loop transfer function matrices because of the simplicity 
of the transformation between these two matrices and hence can only be 
used to design control systems for square plants. 

The design of the controller with this method is done in two 
major steps. The first step is to design a compensator which makes the 
inverse of the open-loop and the closed-loop transfer function matrices 
"diagonally dominant". (Diagonal dominance is defined later but in general 
terms it means that the off-diagonal elements in the transfer function 
matrix for the augmented plant are small with respect to the diagonal 
elements.) A compensator, K(z), that will produce diagonal dominance 
can be developed in a systematic, interactive manner using tne inverse 
Nyquist array. (The inverse Nyquist array is a set of me Nyquist 
diagrams representing every element of the inverse of the OLTFM, where m 
is the number of the plant outputs or inputs.) Once a satisfactory degree 
of diagonal dominance has been obtained the next step is to use single 
variable frequency domain techniques to design a stable controller for 
each pair of input-output variables in the augmented plant. This method 


is discussed more fully in Chapter Three. 


2.3.3 Commutative Controller Method 

Another method which makes use of single-loop design techniques 
to design a multivariable control system is the commutative controller. 
This method, proposed by MacFarlane [33,36], can be very useful in the 
design of a control system for square plants which can be expanded into 


the following dyadic form: 
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g,(Z) = jth eigenvalue or characteristic value of G(z) 
Ww. = linearly independent, frequency independent, real 
eigenvectors of G(z) 
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dyadic matrix, vam, = Say, 
Single-loop techniques can be used to design a multivariable control 
system for this kind of plant. The design of the controllers is done in 
the modal or eigenvalue framework of the plant matrix G(z). Then by 
transforming back to the original set the controller matrix is 
obtained. 

In this method the controller matrix is assumed to be equal to: 
k.(z) wavl em 
J Jad 
Then, the transfer function matrix for the controller plus 


plant will then be equal to: 


m 
Ah i li 
Q(z) = G(z)K(z) = iby g, (2) k (2) WV. (ala) 
and the closed-loop transfer function matrix will be 
: =) - Gi m 9; (z)k.(z) T (2.13) 
R(z) = (1+G(z)K(z)) © G(z)K(z) = bh 149; (2), (@ eV 


The closed loop stability and transient response characteristic 


of the overall system are governed by the properties of the m modal ly 
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non-interacting closed-loop subsystems, 9; (z)k.(z)/(1+g, (2) k5(z)). 
Therefore each subsystem can be designed independently and treated as a 
Single-loop system. 

Unfortunately most systems cannot be expanded in the dyadic 
form given by Equation 2.10 because their eigenvectors are frequency 
dependent. The use of single-loop techniques is not suitable for these 
systems because the closed-loop stability and transient response will be 
governed not only by their eigenvalues but also by their eigenvectors. 
However the commutative controller method can still be applied to these 


systems if its transfer function matrix is approximated by a dyadic 


transfer function matrix using the procedure recommended by Owens [56]. 


2.3.4 Characteristic Locus Method 

A very important contribution in the design of multivariable 
control systems using multivariable frequency domain methods has been the 
generalization of the scalar return-difference quantity introduced by 
Bode [7] to a "return difference matrix". (This approach makes the design 
of multivariable systems equiyalent to the design of single-loop systems.) 
The retuyn-difference and the return-ratio matrices can be calculated 
[38] by assuming that all the loops are broken at one point as shown in 
Figure 2.3. If a signal transform vector a(z) is injected into the feed- 
back control system, as in Figure 2.1, the transform of the returned signal 
at a' is given by -G(z) K(z) H(z) a(z). The return difference matrix F(z) is 
defined as the coefficient matrix of the difference between the injected 


and return-signal i.e. 


peaeemo 2 om (lereG(z)K(2 Hz) Matz) F(z) 7atz)* (2.14) 
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where: 
F(z)= return-difference matrix = I + G(z)K(z)H(z) 


The return-ratio matrix is defined by 


T(z) = G(z)K(z)H(z) (2515)) 
It follows directly that 
E(7jeoe let 2) (2.16) 


and wnen the feedback transducer matrix H(z) is unity 
Tz ee Q(Z) (2.17) 


Two of the methods proposed to date are based on the use of the 
return-ditference matrix to determine the stability of the system. These 
are the characteristic locus method [4,5,41,42] and the sequential return- 
difference method [44,45]. The extended Nyquist stability criterion used 
by these techniques is based on the following relationship proved by 


Popov [61]. 


[F(z)| = ares :Y (2.18) 


The development of the characteristic locus method by Belletruti 
and MacFarlane [4] is based on the use of the characteristic values 
(analogous to the eigenvalues of a matrix of numeric constants) of the 
return-difference matrix. They also developed a more practical approach 
using the characteristic values of the return-ratio matrix [4,42]. 

This method is more general than the inverse Nyquist array 


method [65] since it does not require the system to be diagonally dominant 


be" . 
we sth rye 
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to determine its stability. But the selection of the controller matrix 

is much more difficult because the stability of the system is established 
from the characteristic loci of the return-ratio matrix and it is extremly 
difficult to find a relationship between the elements of the controller 
and the characteristic values of the return-ratio matrix that can be used 
as a guide in the design stage. The design of the controller matrix could 
be done in a more systematic way if a scheme recently developed by Owens 
[59.60] or Kouvaritakis [508] is used. With the characteristic locus method 
it is possible to reduce interactions at high frequencies by aligning the 
characteristic directions or eigenvectors [41,42], of the OLTFM with the 
Standard basis set of vectors. At jfow frequencies interaction is reduced 
by increasing the magnitude of the characteristic loci of the OLTFM [41, 
42]. The method has a provision to check the integrity of the control to mea- 
Surement transducer, error monitoring and actuator failures [4]. (Earlier 
it was pointed out that the performance of a conventional single variable 
feedback loop when applied to a specific input-output pair of variables on 
a multivariable system is a function of the controllers and the other 
elements in other parts of the multivariable system. Thus it is possible 
that a component failure or other change in one part of a multivariable 
system could effect the performance and stability of a anime applied 
to other input-output variable pairs in the same interacting system.) A 
system is said to have high integrity if it remains stable under all 
likely failures [4]. The commutative controller can be considered aS a 
special case of the characteristic locus method that arises when a 
controller is designed to change the characteristic values. of a system 
without changing its characteristic directions (eigenvectors). 


In the characteristic locus approach the controller is calculated 


in stages. At each stage a controller is specified to modify the 
characteristic loci and the characteristic directions of the return-ratio 
matrix or the OLTFM so that it will have the desired properties. The 
final controller is equal to the product of the controllers used in each 
Stage. 
- It should be pointed out that prior to the development of the 
characteristic locus method by Belletrutti and MacFarlane [4] the 
characteristic values of the return-ratio matrix were also used by Bohn 
[502,503] to determine the stability and to design a multivariable 
control system. 

The characteristic locus method is discussed in more detail 


in Chapter Four. 


2.3.5 Direct Nyquist Array Method 

One approach that has not been fully discussed in the literature 
[36,38,70,71] to date is the "direct Nyquist array method". This is 
another technique which gives the control engineer a means of reducing 
the system interaction to a point where a controller design can proceed 
on the basis of classical single variable techniques. This method can be 
considered as a combination of the characteristic locus and the inverse 
Nyquist array techniques. The stability criterion is derived from the 
theory associated with the former method and with the concept of diagonal 
dominance and Nyquist arrays introduced in the latter technique. 


This method requires that the return-difference matrix, i.e. 
B(Ze-=1 +4Q (z+) (22 19)) 


be diagonally dominant. However it only uses the Nyquist diagrams of 
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of the open-loop transfer function matrix, Q(z). 

The design procedure is similar to the inverse Nyquist array 
method, i.e. the selection of the controller is done in two major steps. 
In the first step a simple controller is design to make the return- 
difference matrix diagonally dominant. In the final step a single-loop 
controller is designed for each loop in the system. This method is 


discussed in more detail in Chapter Five. 


2.3.6 Sequential Return-difference Method 

The sequential aetanmmedlik evans method [44,45,46,51,57] in its 
general form is not much difyerent from the sequential method discussed 
in section 2.1 to design a multiloop control system. The main difference 
is that in the general case it is possible to design a pre-compensator at 
each step simultaneously with the controller for each loop. The pre- 
compensator is restricted to perform elementary column operations and, at 


the ith stage of the design procedure, has the following form: 


@)(z) = Je(z)l=1 (220) 


where X defines the column operation to be performed. Then the input 


vector when the first i loops are closed is equal to: 


u(z) = -G.(z) KY(z) y(z) + r'(z) (2.22) 


and the closed-loop transfer function matrix is given by 


22 
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y(z) = (I, + G(z) G(z) K'(z))7! G(z) ri(z) (2.22) 


or 


Viz) =F, we(z) riz) = Gi(z) rz) (2.23) 


j 


In this case the stability of the system is established from the return- 


difference matrix 


F(z) = 1, + 8(z) @.(z) K'(z) (2.24) 


or: 
Fi(z) = 1, + G(z)[6, Ge... G1] K'(z) (2.25) 


The sequential choice of gl does not affect F(z) for j > i. The design 
procedure is as follows. On the ith step the controller k,(z) is designed 
based on the TREUNz) element of g'-! (2). Then a new transfer function 
matrix Gh) is Pecuated with a new compensator matrix @ and the first 
ij loops closed. A new controller kKi4, (2) is designed for the i+] loop 
based on the Te eka element of g'(z) such that this loop is stable. 
The procedure continues in a similar way for the rest of the loops. The 
matrix 6" (z) at the ith step can be calculated using an algorithm 
developed by Mayne [45,46]. | 

The use of a pre-compensator does not solve the main problem of 
the sequential return-difference method which is the order of loop 
closing. Any arbitrary choice of ordering may limit the achievable 
performance of the control systems. Owens [57] has proposed a scheme to 


solve this problem by expanding the plant transfer function into a dyadic 


form. However his scheme requires the use of a post-compensator which 
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is highly undesirable in most practical cases. 

The sequential return-difference method has flexibility in that 
it can be used in combination with any other design method. Any method 
can be used to design the pre-compensator G(z) and the sequential return- 
difference method can be applied to select the final controllers in each 
loop. A very practical combination appears to be the direct Nyquist array 
and the sequential design method. 

A modification of the sequential design method has been 
proposed by Owens [58]. The main difference between these two schemes is 
that Mayne's method [44,45] produces control signals ui (z), us(Z), vers 
u. (z) based on the feedback error e,(z), e,(z) Siig e. (z) respectively. 
2.3.7 Qwen's Modification of the Sequential Return-difference Method 

Owen's modification [58] introduces sequentially at each stage 
the effect of each controller u,(z), u,(Z) Ay u(Z) while retaining full 
output feedback. A schematic diagram of this method is presented in 
Figure 2.4. Using the latter approach the plant and the controller are 


expressed in the following form: 


T | 

G (z) = 2, 9; (2) es (2.26) 

) A: 

gi) = a e sk (2) (2.27) 
= 


where 


g, (z) = jth column of the plant transfer function matrix 


k 5(z) jth row of the controller matrix 


er 
J 


column vector whose jth element is unity and the remaining 


elements are zero. 


For the ith step, let r'(z) = K'Ery(z)s ro{z) ... ri(z), 0... O] 


then the plant transfer function matrix is given by 


y(z) = 6'(z) u(z) (2.28) 


For a negative feedback control: 

u(z) = -K' (2) y(z) + r'(z) (2,29) 
and the closed loop transfer function matrix is 

y(z) = (I, + 61(z) K'(z))7! Bl(z) ¥!(2) (2.30) 


The stability of the system is determined from the return-difference 


matrix, 7.e. 
+ G'(z) K'(z) (22319) 


or 


= Ae 
F(z) = Fy_y(z) +.95(2) ky(2) = Fy_y(2)(Iy * Fpaygy (20K (21) (292) 
Thus the determinant of the return-difference matrix is equal to: 


IF, (2)| Sa ee iz) ec Lee Fy! (z) g.(z) k.(z)| (2.33) 


and applying some formulas for partitioned determinants given in 


Gantmacher [19]: 
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[Fy(2)1 = legey(2)] 11 + ki(z) Fy h yg, (2) | (2.34) 


Owens' method is based on Equation (2.34). It is similar to Mayne's 
method [44,45] in that the design of the controller is done sequentially. 
At each step the m components of the row vector k; (z) are chosen by 
examining the scalar term (1 + k; (z) Fr 4g, (z)) which must satisfy the 
Nyquist stability criterion. The weakness of this method is that it does 
not solve the basic problems of the sequential return-difference method 
j.e. order of loop closing, and selection of the controller constants to 
Satisfy an overall system performance criterion. The modification 
introduces extra design flexibility but it also makes the selection of the 
controller much more difficu.t. The only guideline the user has is that 
the scalar term (1 + k.(z) Foy (2) g.(z)) must satisfy the Nyquist - 


Stability criterion. 


2.3.8 Another Modification of the Sequential Return-difference Method 

When comparing Owen's modification to Mayne's Sequential Return- 
difference method [44,45] it becomes obvious that instead of using the full 
output to generate sequentially us (z), u,(Z) ae u,(2Z)» the feedback 
error e, (z) (i=1, ... m) can be used to simultaneously generate all tha 
control signals. Since the sequential methods are not discussed elsewhere 
in this thesis the modification is presented here. A schematic diagram 
of this modification is shown in Figure 2.5. The control signals 
generated by every feedback error are sequentially accumulated. 


For this case the controller matrix can be represented in the 


following form 


k .(z) e} (2.35) 


where: 


k. = jth column of the controller matrix. 


At the ith stage the closed loop system is given by 
u(z) = K'(z) (y(z) - r(z)) + u'(z) (2.36) 
and 
y(z) = (I, + Gz) K'(z))"! @(z) ul(z) + (Iq + 82) K'(2)P 
@(z) k'(z) r(z) (2.37) 


The stability of the system will depend upon the return- 


+ 
k ;(2) e, (2.38) 


or 


re -] if 
F; (z) = Fs _(z) (I, + ae G(z) k. e. ) (2.39) 
The determinant of the return-difference matrix is 
lpe(z)| = |Fs_q(z)| [t. + Fil Glz) k, el] (2.40) 
j i-| m i-] eal ; 


Similarly applying some formulas for partitioned determinants [19]: 


IF, 
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This modification can be developed from Equation (2.41). The 


design is done sequentially. At the ith stage the m elements of column i 
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are chosen such that the scalar term 1] + er Ba G(z) ks satisfies the 
Nyquist stability criteria. This modification has the same flexibility 


as the Owen's approach [58] but it also has the same problems. 


2.4 CONCLUSIONS 

Of all the methods discussed in this study the following appear 
to be the more promising ones: the inverse Nyquist array [65], the 
characteristic locus method [4,5,41,42] and tne direct Nyquist array 
method [developed in Chapter Five]. (With the exception of the inverse 
Nyquist array method, these methods can also be used to design a control 
system for non-square plants but the selection of the controller is more 
difficult.) Because of their flexibility and convenience it appears that 
these multivariable frequency domain design techniques will be widely 
used in the future. Therefore these tecnniques are described more fully 


in the following chapters. 


CHAPTER THREE 
THE INVERSE NYQUIST ARRAY METHOD 


3.1 INTRODUCTION 

The inverse Nyquist array method first proposed by Rosenbrock 
[65] has received considerable attention recently [1,17,20,21,22,23,40, 
50,67] and several applications of this method have been presented in 
the literature [3,41,48,49,52]. Experience with experimental implemen- 
tations using this method is very limited and only a very small number 
of experimental cases have been reported in the literature [3,39]. For 
these reasons and in order to compare this method with other multivari- 
able frequency-domain techniques, the evaluation of the inverse Nyquist 
array method was undertaken as part of this study and included: 

- the theoretical aspects 


- the computer implementation of this method 


the design of an actual control system using this technique 


the implementation of the resulting controller on process 
equipment 

This chapter deals specifically with the theoretical aspects 
of the inverse Nyquist array method and with the design of a control 
system for a double-effect evaporator pilot plant. The experimental 
results are presented in Chapter Six. 

The inverse Nyquist array method for continuous systems has 
been discussed in detail by Rosenbrock [65]. This chapter discusses the 
most important aspects of the inverse Nyquist array method for discrete 
systems. The inverse Nyquist array method is considered here to be a 


special case of a more general method: the inverse Nyquist method. The 
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objective in doing this is to emphasize the similarities with the 
characteristic locus method [4,42] and especially with the direct 
Nyquist array method discussed in Chapters Four and Five respectively. 
Control systems for the double-effect evaporator pilot plant 
were designed using two different order state-space models of the same 
pilot plant. The results of the experimental implementation of the . 


resulting controller in this chapter are discussed in Chapter Six. 


Oe eae HEeINVERSESNYQUIST: METHOD 

The inverse Nyquist method and its special case the inverse 
Nyquist array method can be applied to the design of a multivariable 
control system of the form shown in Figure 2.1. Since the feedback 
matrix H(z) can always be made part of plant transfer function matrix 
G(z), without loss of generality, it will be assumed that the feedback 
matrix, H(z) is a unity matrix. 

McMorran [47] has proved that for systems which can be 
represented by Figure 2.1 the following relationship holds: 


In! (z)| _ cLeP. (3.1) 
io t(zy] Ober ; 


where: 
CLCP = closed-loop characteristic polynomial 
OLCP = open-loop characteristic polynomial 
Q(z) = open-loop transfer function matrix = G(z)K(z) 
R(z) = closed-loop transfer function matrix = (I + Q(z))7!o(z) 
This was developed independently by McMorran [47] but it is 


equivalent to an equation derived earlier by Popov [61] and later by 
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Hsu and Chen [24]. The main reason for using the inverses of the open- 
and closed-loop transfer function matrices is the simplicity of the 


transformation between these two matrices i.e. 
-] * -| 
Pe) = 1 eae (3.2) 


The use of Equation (3.2) instead of the equation which relates the open- 


and closed-loop transfer function matrices i.e. 
R(z) = (I, + Q(z))"'Q(z) (3.3) 


simplifies, as will be seen later, the design of a control system because 
the eigenvalues of Rl (2) may be obtained by simply "shifting" the 
eigenvalues of q7!(z). 

A general design method can be developed in terms of the 
determinants of the inverses of the closed and open loop transfer 
function matrices by using the following result. Let D be a contour in 
the complex z-plane as shown in Figure 3.1 consisting of a unit circle 
and a circle of radius a which is large enough to ensure that any zero 
of 1a !(z)| and ip! (z)| outside the unit circle is inside the contour 
Deeeretaa (sc ia (z)| map D into a and [uueal map D into iv encircling 


the origin counterclockwise n_ and Ne times respectively (clockwise 


0 
encirclements are assumed to be positive). Using Equation (3.1), the 
principle of the argument [66] and the well known fact that a closed 
loop multivariable discrete system is stable if and only if the closed- 
loop characteristic polynomial does not contain any zero outside the 


unit circle centered at the origin it can be shown [47] that the system 


is closed loop stable if and only if 


n.-n.= (3.4) 
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where Po is the number of zeros outside the unit circle of the open-loop 
characteristic polynomial. 

Although Equation (3.4) can theoretically be used to design a 
multivariable regulator it has no value as a practical design method. 
The controller cannot be designed in a systematic way because of the 
difficulty of predicting the effect of specific elements of the control- 
ler on the values of 1a! ¢z) | and in” '(z)J. 

A more practical approach can be developed by using the 
eigenvalues or the characteristic values of go! (2) and R! (2) to determine 
the stability of the system instead of their determinants. The term 
"characteristic values" has been used by MacFarlane [38] to emphasize 
the fact that they are not constant values (as the familiar eigenvalue 
is) but are functions of z. In this case if every characteristic value 
of a7! (z) and Ro! (2) maps D into ae and ie (Sls2anee mi encine hing 
the origin counterclockwise N64 and na times respectively, it can be 


shown [4] that 


i 
n= ” Ws (3.5) 
and 


m 
n= i nN. (3.6) 


m m 
L CCl fe ole (3.7) 


The approach based on Equation (3.7) can be simplified further 


because the stability of the system can be determined by using only the 
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characteristic values of the inverse of the open-loop transfer function 
matrix (q.(z), j= 1...m). By applying the eigenvalue shift theorem 
[35] to Equation (3.2), the characteristic values of the inverse of the 
closed-loop transfer function matrix (r.(z): j = 1... m) can be expressed 


as: 
r.(z) eae q, (Z) (3.8) 


Then, ustng Equation (3.7) it can be stated that a necessary and 
suffietent condition for the closed-loop system to be stable ts that the 
number of enetrelements of the critteal potnt (-1,0) minus the number of 
enetrelements of the ortgin by the characteristic loci (Nyqutst diagrams 
of the charactertstic values) of matrix ata) be equal to niga 

This stability criterion is the basis of a design method which 
is equivalent to the characteristic loci method proposed by Belletrutti 
and MacFarlane [4] and discussed in Chapter Four. The main difference 
between these two methods is that the characteristic loci method uses 
the characteristic values of the open-loop transfer function matrix or 
in general of the return-ratio matrix while this approach uses the 
Characteristic values of the inverse of the open-loop transfer function 
matrix. Both methods have similar properties and similar disadvantages. 


These are discussed in more detail in Chapter Four. 


3.3. THE INVERSE NYQUIST ARRAY METHOD 

Rosenbrock's inverse Nyquist array (INA) method [65]-is a 
simplification of the general inverse Nyquist method. The INA method js 
based on a more practical stability criterion which can be obtained for 


systems represented by Figure 2.1 when the inverses of the open-loop, q7! (z) 
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and the closed-loop, Rp! (2), transfer function matrices are row or column 
diagonally dominant. The matrix 7! (z) is said to be row dominant if for 


all the z on the contour D the following equation is satisfied: 


, age 
\a5g(Z)1- 2 1 945 (Z)| > 0 (3.9) 


where [455] is the (i,j) element of gy! (2). Similarly the matrix Rg! (2) 
is row dominant if Equation (3.10) is satisfied, i.e. 


m 
[rs (z)| 3 if Ir; 5 (2) ra 0 (3510) 


jf 


where [r; 4] is the (i,j) element of plz), Equivalent equations can be 
written when a7! (z) and Rg! (2) are column dominant. 

When matrices 7! (z) and RH (z) are diagonally dominant the 
Stability of the closed-loop system can be established using only the 
diagonal elements of matrix a7! (z). Rosenbrock [65] has shown, using 
Gershgorin's theorem [43] that if the diagonal elements, q..(z), of 
g7!(z) (not the characteristic values as in the general case) map D into 


[_.. and the diagonal elements, raa(z), of RH (2), map D into r 


O11 Gia 


and ae times respectively, 


11 


encircling the origin counterclockwise n).. 


then 
m 
pea va a eae (ea) 
and 
m 
th. ee oe bean (312) 


Thus using Equations (3.4), (3.10) and (3.11) it can be shown that a 
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suffietent and necessary conditton for a dtagonally dominant closed- 


loop system to be stable ts 
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Because r.. (2) =] + q.4(z), the only difference between the Nyquist 
diagrams of r;.(z) and q,,(z) is a shift in the imaginary axis. Thus the 
system is closed-loop stable if the number of encirclements of the (-1,0) 
point minus the number of encirclements of the origin by the Nyquist 
diagrams of {q. 5: EE 2 owes Ns pale: is equal to -p,. 

Based on this stability criterion Rosenbrock [65] has developed 
a very practical technique to design a multivariable control system. 
Rosenbrock [65] has also introduced the concept of the inverse Nyquist 
array which is a set of mé Nyquist diagrams (m is the number of inputs or 
outputs in the system) that correspond to every element of a7 lz). The 
inverse Nyquist array is a very useful graphical tool to observe the 
effect of each controller matrix in the control system and to select the 
type of control which is necessary. 

Diagonal dominance in the system can also be established 
graphically by drawing circles for different frequencies on the Nyquist 
diagrams of the diagonal elements of gy! (z). The centres of these 
circles are located on the Nyquist loci of the diagonal elements, q.,(z)s 
of q7!(z) and their radii are equal to the sum of the magnitudes, for a 
given frequency, of the off-diagonal elements in a row or in a column 
depending whether row or column dominance is to be established. 


For row dominance the radii of the circles are equal to 


idee tz) (3.14) 


m 
5 ae, I ACaR (3.15) 
7 


The bands formed by these circles are usually referred as the 
Gershgorin bands. The inverse of the open-loop and the closed- 
loop transfer function matrices are diagonally dominant if the Gershgorin 
circles do not enclose the origin and Pe (-1,0) point respectively. 
The Gershgorin bands have two more important properties. The first one 
has been proved by Gershgorin [43] and states that the characteristic 
values of the matrices q7!(z) and Ro! (2) are located inside the 
Gershgorin bands. This property was used by Rosenbrock [65] to prove 
the stability criterion on which the inverse Nyquist array method is 
based. It is because of this property that the Gershgorin bands are 
also used to determine the stability of the system. 

The second property is a very important one because it gives 
some physical insight into this method. Rosenbrock [67,69] has proved 
using a theorem due to Ostrowski [55] that when the inverse of the 
open-loop transfer function matrix 7! (z) is diagonally dominant the 
Gershgorin bands are the limiting bounds of a narrower band which is 
usually referred as the Ostrowski band. This latter band also has the 
property of being a limiting bound of the region where the inverse of 


the transfer function between input i and output 71, when loop 7 is open 
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and the rest are closed, is located. The Ostrowski bands are generated 
by circles whose centers are located in the Nyquist diagram of the 


diagonal elements, q,;(z) of a! (z) and whose radii are equal to the 


radii of the Gershgorin circles multiplied by a factor a or a These 


factors are equal to 


ds 
a, = max ~— +s (silo) 
+eq. 
J aa | 
if row dominance is used, or 
6. 
a; = max ——J— (3.17) 
“a [1 + q54| | 
ie 


if column dominance is used. 
The latter property of the Ostrowski bands can be represented 


by the folitowing equations: 


1a; 5 (2) - hy! (z)| < od. ace Gor 


or 


is -| 1 
Iq. 5 (Z) 7 h. (z) | < wee < OE (3.19) 


where, h, (z) is the transfer function between input 7 and output 7 when 
the loop i is open and the rest of the loops are closed. 

Since the width of the Ostrowski bands depends upon the values 
of the controller constants in the other loops they are not generally 
used in the design of the controller except when it is necessary to 
define more precisely the region of stability of each loop. Its use 


requires a trial and error procedure. However the property of the 
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Ostrowski bands of being located inside the Gershgorin bands is very 
useful because once the system is diagonally dominant the design of the 
controller can be done using single-variable techniques on each of the 
Gershgorin bands of the inverse of the open-loop transfer function matrix. 
The Gershgorin bands can also be considered as a measure of 
the degree of interaction in the system. The narrower these bands are 
the less interacting the system is. 
The design of a multivariable regulator using the inverse 
Nyquist array method involves three steps: 

1) Pairing of the input-output plant variables. 

2) Design of a pre-compensator to make the inverse of the open-loop, 
Q-1(z), and the closed-loop, R-!(z), transfer function matrices 
diagonally dominant. 

3) Design of a multiloop control system using single-variable 
frequency domain techniques. 

In the first step the inputs are reordered in such a way that the 
interaction between the input-output systems fu, (z): y, (z) licibe cepa ets 
minimized or that the control of the output y,(z) is done by the input 
that has the greatest influence on it, u.(Z). This step is usually done 
by graphical comparison using the plant inverse Nyquist array which is 
also useful to determine if a compensator is required to make the system 
diagonally dominant. 

In the second step a pre-compensator is selected in order to 
make matrices g7! (2) and R7! (2) diagonally dominant. This compensator 
is usually designed by performing successive elementary row operations 
on the inverse of the plant transfer function matrix. Most of these 


operations can be selected by visual inspection of the corresponding 


inverse Nyquist array. 


Acker 


eye 


The final control matrix will be equal to the product of the 
control matrices found in each step. Each step is an iterative and 
interactive process. However the design of the control matrix is always 
done in a systematic way. 

Once the matrix gq”! (z) 1s diagonally dominant a multiloop 
control system is designed by using the single-variable frequency domain 
techniques which apply to inverse polar plots. The only difference 
between the design of single-loop control systems and the design of a 
controller for each loop in this step is that the Nyquist diagram of 
the inverse transfer function of each loop is replaced here by a band 
(Gershgorin band). | 

This three step procedure also allows the designer to solve 
several of the control system requirements like stability, reduction of 
interactions, integrity against transducer and error-monitoring failures 
and the performance of the system. In the first two steps a system is 
given good non-interacting and high integrity properties against trans- 
ducer and error monitoring failures by selecting a controller which 
makes the inverse of the open-loop and closed-loop transfer function 
matrices diagonally dominant. For example, if a failure in the transducer 
of any loop in the system occurs all the remaining loops in the system 
will remain stable. Although such a failure will affect the transfer 
function of each loop when it is open and the remaining loops are closed, 
the Nyquist plots of these transfer eae will still be located 
inside their respective Gershgorin bands. Since the design was.based on 
the Gershgorin bands the stability of the loops will not be affected 
( although the gain margins may be changed). The same occurs when there 


is an error-monitoring failure. 
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In the last step of the design the feedback gains {k,,i=1,..m } 
are chosen so that the closed-loop system is made stable and the 
performance of the system is satisfactory. This makes the inverse 
Nyquist array method a very useful and practical technique. 

There are however some shortcomings in the use of the inverse 
Nyquist array method. These are: 

1. It is not always possible to find a simple controller matrix 
which will make gt (2) and Rv! (2) diagonally dominant. For the 
case when the matrix plz) has been made diagonally dominant and 

an analytical expression of G(z) (or Q(z)) is given it is 
possible to overcome the requirement that gq (2) be diagonally 
dominant at the cost of additional calculations. It is always 


possible to go back to the general stability criterion given 


by Equation (3.4), i.e. 
) 


Since the inverse closed-loop transfer function matrix, R*(z), 
is diagonally dominant, No will still be given by Equation 
m 
(3.12) (n, = &2n_.; ). The encirclements of the origin, n 
Co sa Cli 
can be obtained directly by examining the Nyquist diagrams 


aye 


corresponding to the mapping of the contour D by |Q(z)| or 
HG Gz) te 

There is also the alternative of modifying Equation (3.4) 
by taking into consideration the fact that Wa 
where Zo is the number of zeros of |Q(z)| outside the unit 


circle. Equation (3.4) then will be equal to: 
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or 


From this equation it can be seen that the closed-loop 
system will be stable if and only if the number of encirclements 
of the (-1,0) point by the Nyquist diagrams of {q,;(Z), j=1,..m} 
is equal to the number of zeros of |Q(z)| outside the unit 
circle. 

In some cases g7 (2) can be made diagonally dominant by 
designing a post-compensator, L(z), in combination with a pre- 


compensator, K(z), such that gt (2) is equal to: 


ghz) = K4*(z)e@72(z)L 4 (z) (3.20) 


From the practical point of view this alternative is questionable 
because the control would then be exercised over a combination 
of the output variables instead of the output themselves. 

Another alternative has been used by Hawkins [S02 505 |e 
A feedback controller, K(z), can always be found that will 
make the closed-loop sytem containing Q(z), diagonally dominant. 
This closed-loop system can then be treated as the new system 
matrix Q, (2) and the design continued using the standard 
inverse Nyquist array method. 


Since diagonal dominance is not a necessary condition for 
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non-interaction [42] or good control, some design feedom 
(and/or control energy)” may be wasted. making the system 
diagonally dominant. 

3. Integrity against actuator failures cannot be checked in this 
method. 

4. It can only be applied to plants which have a square transfer 


function matrix. 


3.4 USE OF A BILINEAR TRANSFORMATION 
Althougn all the results presented here for discrete systems 
have been obtained using the z-transform, the design of a multivariable 


regulator for a discrete system is, for practical reasons, usually done 


using a w- or r-bilinear transformation, i.e. 


rh Sa omc rz) 
or 
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A ame (eee) 


A bilinear transformation maps the unit circle in the z-plane 
into the left half plane in the w- or r-plane. Its use makes the 
computation of the Nyquist diagram easier because the Nyquist contour 


used in the w- or r-plane is identical to the Nyquist contour used for 
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continuous systems in the s-plane as it is shown in Figure 3.2. This 
also simplifies the programming to implement this technique because the 
Same programs can be used to design a continuous or a discrete control 
system. The Nyquist diagrams obtained using a bilinear transformation 
are exactly the same as the ones obtained using a z-transform. 


In this study the w-transformation was arbitrarily chosen. 


3.5 AN ALTERNATIVE TO THE INVERSE NYQUIST METHOD 
Similar stability criteria can be obtained using the open- and 
closed-loop transfer function matrices instead of their inverses by 


rewriting Equation (3.1) in the following form: 


R(z)t ~~ OLCP (3.23) 


Using this equation equivalent stability criteria can also be 
obtained in terms of the eigenvalues of Q(z) and R(z) and in terms of 
their diagonal elements when these matrices are diagonally dominant. 
However the methods based on these stability criteria are less practical 
than the inverse Nyquist array discussed in this chapter and the direct 
Nyquist array methods discussed in Chapter Five. This is because the 
methods based on Equation (3.23) require the calculation of the closed- 
loop transfer function matrix using Equation (3.3) to determine the 
Stability of the closed-loop system. This calculation is not necessary 


in the inverse and the direct Nyquist array methods. 


3.6 DESIGN OF A MULTIVARIABLE REGULATOR FOR THE THIRD ORDER EVAPORATOR 
MODEL 


In this section a multiyariable regulator for the double-effect 
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evaporator pilot plant, described in appendix A, is designed using the 
inverse Nyquist array method. The evaporator was represented by a 
a discrete, three state, linear, time-invariant model with normalized 
perturbated variables. (The coefficient matrices for the state-space 
model are shown in Table A.3) This model was obtained by reduction, 
using Marshall's method [77], of a ten-state linear, time-invariant model 
derived by Wilson [76]. | 

The z-transform of this model can be represented by the 


following equation: 


W, (z) S(z) 
Wo (Z) = G(z) B, (z) (3.24) 
C,(z) B, (Zz) 


The outputs of the system are the first-effect hold up, Wis the second- 
effect hold up, Wo» and the product concentration, Cos The inputs are 
the steam to the first-effect, S, the bottoms flow from the first-effect 
B> and the bottoms flow from the second-effect, B,. The plant 

z-transfer function matrix G(z) was obtained from the state space model 
using the algorithm of Souriau-Frame-Faddeev [66]. G(z) is presented in 
Table 3.1. For reasons explained earlier the design was carried out 
using the w-bilinear transformation of G(z). The plant w-transfer 
function matrix of the third order evaporator model, G(w), is shown in 
TaDilem3e2. 

The poles of the z- and w-transfer function matrices or the 
zeros of the open-loop characteristic polynomial of the third order 
double-effect evaporator are shown in Table 3.3. 


The evaporator model is marginally stable since it has two 
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TABLEas3 32 


w-transfer Function Matrix Derived from the 


Third Order Evaporator Model 


-0.01628 -0.04055 0 
W W 
G(w) = (1-w) -0.01890 0.04275 -0.0203] 
W W W 
0.02690 -0.02252 0 
(w+0.0203) (w+0.0203) 
TABLES irs 


Poles of the Tnird Order Evaporator Model 


z-plane w-plane 


then) 0.0 
aC 0.0 


0.9602 -0.0203 


TABLE 3.4 


Inverse of the w-transfer Function Matrix Based 


on the Third Order Evaporator Model 


-15.444 w 0 27.811 (w+0.0203) 
pean (ay -18.508 w 0 -11.650(w+0.0203) 
204.518 w -49.150 w -49 225 (w+0.0203) 
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poles located on the unit circle in the z-plane. 
Since the inverse Nyquist array works with the inverse of the 


plant transfer function matrix, 67! (w) is presented in Table 3.4. 


Step one: Pairing the input-output variables of the evaporator 

The inverse Nyquist array of G(w) is shown in Figure 3.3 for 
the frequency range of 0.1 - 10.0 radians/second. The labels in this 
figure indicate the numerical value at the end of each axis. The 
pairing of variables has been done by graphical comparison between the 
elements of the Nyquist array of onl (w). It can be seen in Figure 3.3 
that none of the rows or columns of 67! (w) are diagonally dominant. This 
can be seen more clearly in Figure 3.4 where the diagonal elements of 
om! w) with the Gershgorin bands for row dominance are displayed. Thus 
the plant in its present form is highly interactive. 

In order to reduce interaction and to have more effective 
control of the outputs a renumbering of the inputs is necessary. From 
Figure 3.3 it can be seen that the best configuration is obtained when: 

row 2 is made row 1 
row 3 is made row 2 


row 1 is made row 3 
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This operation is accomplished by pre-multiplying ar lw) by the 


controller matrix ? 
0 
ts Da, tel (3.25) 
0 0 


Note, when using the interactive design program developed at the 
University of Alberta the user simply identifies what elementary 
operations he wants done and the computer generates the appropriate 
compensator, e.g. a in Equation (3.25). 

This means physically that if a muZtiloop control system is to 
be designed the least interacting scheme occurs when the first effect 
hold ups Wie is controlled by Bis the second effect hold up> Wos is 
controlled by B., and the product concentration C,, is controlled by 
steam flow, S. This is in agreement with previous results obtained by 
Newell [53] and Jacobson [27] when designing a multiloop control] system 


for the evaporator using conventional techniques and experimental data. 


The new open-loop transfer function matrix Q, (w) is given by 
-] Beal nal 
Q, (w) = K G (w) (2376) 


Its inverse Nyquist array is shown in Figure 3.5. The diagonal elements 
of CD with their Gershgorin bands are presented in Figure 3.6. From 
these figures it can be observed that only rows 1 and 3 of ar! (w) are 
diagonally dominant. 

Because row 2 of aw) is not diagonally dominant there is a 
need to proceed further with the design procedure,i.e. to design a 


controller which makes the system diagonally dominant in all three rows. 
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Step two: Design of a pre-compensator to make matrices a7! (w) and 
po! w) diagonally dominant 
The Nyquist array of 07! (w) in Figure 3.5 helps in the design of 
compensator K(z) that will make Q(z) diagonally dominant. For example, by 
comparison of the magnitude q,, (2) with qo, (2) (and/or by trial and error) 
it can be seen that subtracting from row 2 1.325 times row 1 will make 
row 2 diagonally dominant. Simultaneously, diagonal dominance can be 


improved in row 3 by subtracting from it 0.830 times row 1. Then the 


inverse controller is given by: 


1.0 0.0 0.0 
5 {pepe | i 0.0 (oez7p 
DE OA leo 


iS 


The inverse Nyquist array of Q,(w) (where 05 | (w) = G7! (w)) 
is presented in Figure 3.7 and the diagonal elements of a! w) with the 
Gersngorin bands are shown in Figure 3.8. All the rows in Q5'w) are 
now diagonally dominant. However diagonal] dominance can still] be 
improved in row 1 and 2. This can be done simultaneously by adding 0.3 
times row 3 to row 1 and 0.930 times row 3 to row 2 respectively, i.e. 
Q5'(w) must be pre-multiplied by: 
a0 0.0 0.300 


K, = 0.0 20 0.930 (2828) 
0.0 0.0 hot) 


to obtain Q3'(w). From the inverse Nyquist array of a5! (w) shown in 
Figure 3.9 it can be observed that the compensated system is almost 
diagonal, i.e. almost totally non-interacting. 


Since the system is diagonally dominant the stability of the 
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closed-loop system can be investigated using the inverse Nyquist array 
method. If the open-loop system is assumed to be stable (poles at the 
origin assumed to be stable), then the contour shown in Figure 3.2 can 
be used and Poe the number of zeros of the open loop characteristic 
polynomial in the right half plane is equal to zero. According to 
Rosenbrock's [65] stability criteria the closed-loop system will be 
Stable if the diagonal elements of the inverse of the open-loop transfer 
function matrix encircle the origin counterclockwise the same number of 
times as they encircle the point (-1,0). 

To satisfy this encirclement criterion it is necessary to 
change the signs in loops |] and 2. This can be done using the following 


controller: 


ine eee a 
ie 2 i. Sager ae (3.29) 
0 a LG 


The open loop transfer function matrix at this stage is Qy(w) >» 
where CD = Kn Q3|(w). The inverse Nyquist diagrams of its diagonal 
elements are shown in Figure 3.10. The pre-compensator which has been used 
to make the open-loop transfer function matrix diagonally dominant and stable 


is given by Kas where: 
Ky (3.30) 


and 


Ke = = (5930 Zs090 gee 160 (ors) 


or 


5 | 


-0.830 0.0 Ory oy 
Kp = =|) i088 0. 0s -02300 (aro2) 
plo 2 49 ee 1 0net—e327 


Step three: Design of a multiloop control system for the open-loop 


transfer function matrix 
The final step in the design procedure is the selection of the 
final controllers in each loop. Because the transfer function between 
the output and the input for each loop is represented by a first order 
system the use of only proportional plus integral controller in each 
loop should be sufficient to produce good control. Physically, this 
means that a controller matrix K. (z) should be selected such that the 


final multivariable regulator is given by 
K(z) = Kp K. (2) (3233) 
where K.(2) is given by Equation (3.34). If only proportional.control is to 


be used: 


a 0 ko 0 (3.34) 


When proportional and integral control is desired K.(z) is given by 


Equation (3.35) i.e. 


K-(z) = 0 ko ti 0 (3.35) 


Je 


The selection of the controller constants of matrix Ke can be 
either done using simulated runs or by direct tuning of the controller 
constants in the field. Although the selection of the controller constant 
was, in this case, done by experimental tuning, the ee of simulated runs 
could be helpful in some cases. 

A very useful piece of information to guide the selection of the 
the final controller constants in matrix Ke is the estimation of the 
stability gain region for each loop (i.e. the gain margins). This can be 
obtained graphically from the diagonal elements of the inverse of the 
open-loop transfer function matrix, Qy(w) in Figure 3.10. The closed-loop 


System will be stable and diagonally dominant for the following gains: 


igefeyss I ca 1/k, Ss WARS @p i) < k, < 18.5 (Se 
Loopy2: se > 1/k, | (490 eee Oye Oe ky < 49.2 (33 (5) 
(FOOpmG ramcm I/k3 >ml/S720anor BO. < k, aye ae) (e38)) 


Note that this system has "high integrity" against measurement 
transducer or error monitoring component failures because it remains stable 
and diagonally dominant for a wide range of gains (including zero or "open- 
loop"). According to the Ziegler-Nichols rules the recommended value for the 
controller constants k,, ko and ky is 50% of the ultimate gains given in 
Equation (3.36), (3.37) and (3.38) respectively. These values were found 
to be too high when the regulator was experimentally implemented although 
the simulated runs had indicated the opposite. For this reason a decision 
was made to tune the loops in the field. 

Figure 3.11 presents a simulated response of the closed-loop 
third order, evaporator model for a 20% step change in the feed flow 


disturbance using different gains for Ky > ko and k3- The values used for 
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these gains are 50% (controller: FD0350), 30% (controller: FDO330) 20% 
(FD0320) and 10% (FD0310) of their respective ultimate gains. It can be 
seen that there is an improvement of the closed-loop system response 
when the gains are increased. This improvement is mainly in the control 
of the first-effect hold-up, Wy. 

Figure 3.12 presents the response of the controlled system for a 
20% step change in the feed flow when the gains used in each loop are 
equal to their ultimate gains. As expected the evaporator is on the 
verge of instability. 

In order to evaluate the degree of interaction in the control 
system Figure 3.13 presents the response of the control system for a 10% 
step change in the setpoint of the product concentration, Co. when the 
controller constants k,> Kos Ky are 50% (controller: FD0350) and 20% 
controller: FD0320) of their respective ultimate gains. In this case 
there is not much improvement in the response of the double-effect 
evaporator when the gains Ky > Ks and k3 are increased. However the 
system is, for all practical cases, non-interacting regardless of the 
gains used. The use of the controller (FD0350) for this step change 
violates the physical constraints in variable B.. 

Tne selection of the controller constants for the proportional- 
integral control case was also done by experimental tuning. In this 
case it was not possible to apply the Ziegler-Nichols' rules to select 
the integral constants because the crossover frequency for each loop is 
zero. 

The results of the experimental tuning of the controller 


obtained in this section are presented in Chapter Six. 
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3.7 DESIGN OF A MULTIVARIABLE REGULATOR FOR THE FIFTH ORDER EVAPORATOR 
MODEL 

The design of a multivariable regulator for the double-effect 
evaporator pilot plant, described in appendix A and represented by a 
fifth order model, is presented in this section. 

The z-transform of the system used in this example is also 
given by Equation (3.24). The only difference with respect to the 
example presented in the last section is that in this example G(z), which 
is shown in Table 3.5, was obtained from a discrete, five-state (rather 
than three state) linear, time-invariant model with normalized perturbation 
variables. This model was obtained by Wilson [76] and it is shown in 
Table A.2 in Appendix A. As in the previous example the design was 
carried out using the w-bilinear transformation of G(z). The w-plant 
transfer function matrix, G(w),is shown in Table 3.6. 

The poles of G(z) and G(w) or the zeros of the open-loop 
characteristic polynomial of the fifth order double-effect evaporator 


model are shown in Table 3.7. 


VES S7/ 


Poles of the Fifth Order Evaporator Model 


z-plane w-plane 
ibelé 0.0 
120 0.0 
0.9604 -0.0202 
0.9216 -0.0408 
0.4328 -0.3905 


The inverse of the w-plant transfer function matrix is presented 


ant ath sega a a 


wre ie data ee er 


: UV 


os 


vray ag a patiently softies sl a 
ind Heian ae FO atendssy att _: 
ie aa? }«1S8 Se) oot 
ais SCE TUALARRTS: 4 ici tovete. wor 
ee) a ae athneagah at sit 
‘ ee inetiei bo 47 at ent at) tuo 


nd \“Gginty 2R wee hel stron 98 


itp LA a 1 apis 


he 
at a 


_ a 
ae 


Ne 
: t pene: 4 ‘Zz? ‘Mabta 7 ‘S 


qo, 824! nay Ay esas he 
& — oe Oe ew. & ° t ¢ shies we 


ce wean gia spear cy 


: 


“a lagi ed 45 


wit sion tne 0 » note alt 


25) 


(p8Er°0-Z)(9LZ6°0-Z) (€096°0-Z){(0°1-Z) = dd10 


2 


sq uaAtB SL XLujzew UOLQOUN} UaysuUeu? SLUZ JO LeLWOUALOd OL4SLUazZOeUeYD aUL 


(€096°0-2Z) (v8Ev" 0-2) (9126 °0-Z) (€096°0-Z) 


0 8zer0'0- (8688 °0-Z SZ€10°0 
(0° L-z) (0° 1-2) 86h"0-Z) (O° L-Z ibe 
1S6v0 0- 180 °0 $29/°0+Z) SLELO'0- 

; (oatez) (y8er'0-z 


69L80°0- (8Z9Z°0+Z) L6LLO’O- 


el 


L2pow 4OZeUOdeAR Uapug UZsL4 BY WOUJ PAALUAG XLUzZeW UOLZOUNY Uajsueuy-Z 


Gmc mela 


“oe a P. - 
HN comes if 
iyo 


= aa A - 
aa his oo a — . : ae - ae rr 
_— eyti i ot? move bavined «eet norzomer 
: eye ; — - 
ay : a = = ad : a 
ron t i, a (= - 
AE _ 
“ ae ~ a 
: S ; 
7 ; 
_— , - J 
: ‘on 4 - a 
| e ’ 
\ 
4 ' 
. f y a ‘ af ‘ 
; a, uh 
F , x 
; | ii 4 
; ; Badd Ae 
- : i r 7 
i 
47 y) 
jj i 
by , : 
= 


M 


62020 0 


(2020 °0+) 
80220°0 


M 
becvUs0s 


ver Na 
S800 °0 


LOpow A0ZeUOdLAZ UapugQ YIJFL4 GY WOUJ paALUaq XLUQeW UOLJOUNY AassueA.-M 


(2020"0-+M) (8010 *0+M) (GO6E*0+M) 
(€850°O+") (ZEv°Z+M) ZPLLOO'O 


GO6E “O+M)M 
(2p'L+) YELL00°0 


GO6E "O+M)M 
(Zep °Z+™) 286000°0 


Pts, SU 


A 


eh Rm cm om a teenie 


en 


0 as iy not 


in Table 3.8. The plant in this example is also marginally stable since 


it contains two poles on the unit-circle in the z-plane. 


Step one: Pairing the input-output variables of the evaporator 


The inverse Nyquist array of G(w) presented in Figure 3.14 


Shows that there is a strong interaction between the loops in the system. 


This is also shown in Figure 3.15 which presents the diagonal elements 

of G7! w) with their Gershgorin bands. In this example interaction can 
be reduced by an appropriate pairing of variables or by renumbering the 
inputs. 


As in the previous example the best combination is obtained if 


W is controlled by B, > We 1s controlled by B. and Cy TS controlled by S: 


This change can be implemented by pre-multiplying the matrix 67! (w) by 


eS where: 
lace 
e > 4 (3.39) 
0 


The inverse Nyquist array of CD = ea Gain) is shown in 
Figure 3.16 and its diagonal elements with their Gershgorin bands are 
Shown in Figure 3.17.. From these figures it can be observed that the 
matrix a7! (w) is not diagonally dominant. Therefore it is necessary to 


continue with the second step in the design procedure. 


Step two: Design of a pre-compensator to make the matrices Vy and 


Rl (w) diagonally dominant 
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Step 2a) Design of a constant pre-compensator 
From the Nyquist array of a7! (w) in Figure 3.16 it can be 


observed that rows 1 and 3 of a! (w) are slightly diagonally dominant at 
low frequencies. Row 2 can be made diagonally dominant at low frequencies 
by subtracting from this row 2.0 times row 1 and by adding 0.3 times row 
3 to row 2. This operation can be performed by pre-multiplying a (w) 
by the following controller: 
] 0 0 
-] _ 


Ko = <7. 0 arg OSS (3.40) 
0 0 120 


The new open-loop transfer function matrix is given by 
sl inyet ae ele 
Qo (w) = Ky Qy Ww) (3.41) 


The Nyquist array and the diagonal elements with the Gershgorin bands of 
05 (w) are shown in Figures 3.18 and 3.19 respectively. Diagonal 
dominance can still be improved in row 3 of a5 | (w) by subtracting 2.0 


times row 1 from row 3, i.e. a5! (wr) must be pre-multiplied by 


0 mansG 
Ka ioe ha 1 0 (3.42) 
iy 


The improvements on diagonal dominance in a3! (w) = Re Q5'(w) 
can be observed from Figure 3.21 which shows the diagonal elements of 
a5! (w) with the Gershgorin bands. The Nyquist array of Q; is also 
presented in Figure 3.20. 

Diagonal dominance can be improved marginally in row 1 or 


a5! (w) by adding to this row 0.15 times row 3, or by pre-multiplying 
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a5 (w) by the controller: 


ie Omen Om. ChAT 
ee lhe ao (3.43) 
0 0 1] 


The Nyquist array of the new open-loop transfer function matrix a (w) = 
a Q; is shown in Figure 3.22 and the diagonal elements of this matrix 
with their Gershgorin bands are presented in Figure 3.23. 

In order to satisfy the stability criterion in each loop the 
Sign has to be changed in rows 1 and 2 by pre-multiplying ij by the 


controller matrix: 


Sid 
a nmi ts) 1-15 eto (3.44) 
0 oO 


The diagonal elements of Qe! (w) are shown in Figure 3.24 for a 
frequency range of 0.01 - 1.0, in Figure 3.25 for a frequency range of 
0.01 - 0.5 and in Figure 3.26 for a frequency range of 0.1 - 5.0. From 
these figures it can be concluded that the OLTFM and the CLTFM are only 
diagonally dominant at low frequencies. Consequently the closed-loop 
Stability and the stability region cannot be determined at this stage 
and the design cannot. be continued. In order to make the inverses of 
the open-loop and closed-loop transfer function matrices diagonally 
dominant a dynamic pre-compensator has to be used. 

Several proporticnal-integral controllers were tested but 
none of them make the inverses of the OLTFM and CLTFM diagonally 
dominant. Thus a more complicated controller was necessary to complete 
the design of a control system for the double-effect evaporator. This 


is discussed in the next section. 
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This example is a good illustration of one of several 
pitfalls of the inverse Nyquist method. It is not always possible to 
find a simple controller which will make the inverses of the OLTFM and the 
CLTFM diagonally dominant. Furthermore if diagonal dominance is not 
obtained the design cannot be completed using this design method. (In 
Chapter Five it will be shown that the direet Nyquist array method can 


be applied to this example, to design a proportional control matrix.) 


Step 2b) Design of a dynamic pre-compensatur 

From the Nyquist array of a! (w) and the Nyquist diagrams of 
its diagonal elements in Figures 3.16 and 3.17 it was observed that only 
row 3 was diagonally dominant. It is obvious from Figure 3.16 that 
diagonal dominance at low frequencies could be obtained in row 2 by 


subtracting row 1 from row 2, i.e. by pre-multiplying the matrix Q (w) 


by 
eee eO 
Cape RS) cet (3.45) . 
0 0 


The Nyquist array of (Q(w))| = (ols Q, (w) is shown in 


Figure 3.27 and the Nyquist diagrams of its diagonal elements with their 
respective Gershgorin bands are shown in Figure 3.28. A careful 
examination of the poles and zeros of the matrix (Q(w))7| and of Figure 
3.27 shows that row 3 could be made perfectly diagonally dominant if a 
phase-lag compensator is used and the following elementary row operation 


is performed: 


row 3 = row 3 - 19.2 (MEO) x row (3.46) 
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This operation can be accomplished if (Q3(w))7! is pre-multiplied 


by: 
] 0 0 
1 elo 
(K, (w)) = 0 ] 0 (3.47) 
SS 
The new OLTFM is given by: 
(a3 Qw))7! = (KY(W))7! (QS (w)) | (3.48) 


The Nyquist array of (Q3(w))7! 


and the Nyquist diagrams of its diagonal 
elements with their Gershgorin bands are shown in Figure 3.29 and 3.30 
respectively. Figure 3.29 and an analysis of the poles and zeros of 
(Q3(w))7! indicates that row 1 can be made perfectly diagonally dominant 


if a phase-lead compensator is used and the following operation is 
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performed: 
: w+ 7.400_ 
row 1 = row 1 + 0.0165 oF 0.390 yEXECOWes (3.49) 
This operation can be accomplished by using the following controller 
matrix: 


w+ 7.400 
] OR OS0).65 (+0300 7 
(Kw)! = | 0 0 (3.50) 
0 0 ] 


The Nyquist array of CHO = (Ky ()) | (Qh (w)) and the 


Nyquist diagrams of its diagonal elements with their Gershgorin bands in 


Figures 3.31 and 3.32 show that rows 1] and 3 are almost diagonal. 
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However row 2 is only diagonally dominant at low frequencies. An improve- 
ment in the diagonal dominance of row 2 can be obtained if another phase- 


lead compensator is used and the following operation is performed: 
: w+ 7.400 _ 
row 2 = row 2 + 0.056 i 0-39 ) x row 3 (32s) 


To perform this operation it is necessary to use the following 


controller matrix: 


] 0 0 
“| -] _ w + 7.400 
(Ke (w)) = 0 ] 0.056(7 + 0.390 ) (3.52) 
0 0 1 


The Nyquist diagrams of the diagonal elements of matrix 
(Qe (w)) | = (Ke (w))7| CHO in Figure 3.33 indicate that this matrix 
is almost diagonal. The Gershgorin circles are very small but the 
diagonal elements (loops) 1 and 2 are unstable. To correct this 


situation it is necessary to pre-multiply (Qe (w))7| by 
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The Nyquist diagrams of the diagonal elements with their 
Gershgorin bands of the OLTFM (Qe(w)) | = (ie (Qe (w)) | are presented 
in Figure 3.34. 

The pre-compensator used to make the inverses of the open- and 


closed-loop transfer function matrices diagonally dominant is then given 


by: 


peteaeey Sf5 Dre baru 2t.9oseem 


eryeyvess 24s tore of — Si 
=. - 7 
: . eon wel boa 
-~ ; Wt 
4 7 1 a 
| 
t — 
ae in . 
= a - y ad wr 


3 _ icy \ Pyer 97) 4 ‘ : ft joy ty in " 

= - ‘ - : 2 ‘. Fi J; J 7 sl 

risora yi ie. we LAM ga? ((«) 0 
: : = ~~ 


_ _ 
- 


ga ing ous sifogho. WW Yastivae ah _Abeeieth State gE 
7 


3 sen Wiel aS tee, 7 Cen consie {Segake 


or 
w + 0.390 
-19 G + 7.400 ) 0 0.668 
- b re w + 7.400 
Kp (w) ] 0 0.0165 & $0.390 ) (3.55) 
, i: _ w_+ 7.400 
] ] We0725 Ma $0390 ) 
The z-transform of pre-compensator Kp (w) is equal to: 
z - 0.4383 
-3.144 +7619) 0 0.668 
5 Zatm0e76 19 
Kp (z) = -] 0 -0.0997 (—9-9383) (3.56) 
ZEtoOsLOlg 
-| -] -0.438 (91383) 
Step three: Design of a multiloop control system for the open-loop 


transfer function matrix: 5 (w), 


In this step the controller matrix Ke must be selected such 


that the multivariable regulator is equal to: 


iS & Kp Ke 


(Beol)) 


If only proportional control is desired Ke will be given by: 


kK, 0 0 
Ke = 0 ky 0 
0 0 ky 


(2550) 


When proportional plus integral control is required Ke is given by 


Equation (3.35). 


In order to select the final gains in each loop it is very use- 


ful to obtain the gain margins of each loop for which 


the system will 
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remain closed-loop stable. This piece of information can be graphically 
obtained from Figure 3.34. The control system will be stable for the 


following set of gains: 


Loop 1: 15.6 < kK, eo OG Un kK, 2 NE (3.59) 
] 

Loop 2: Roepe ehoec mort One ky < 49.2 (3.60) 
] 

Loop 3: 0.0 < k, a (ig M0) < k, < 50.9 (S261) 


These gain margins are almost the same as the gain margins 
obtained in Section 3.6 using the third order model with the exception 
of the third loop in which there is a difference of less than 20%. This 
makes the comparison between the controller matrices obtained in this 
section and Section 3.6 easier. 

The simulated response of the fifth order modei for a 20% step 
change in the feed flow and using the controller obtained in this 


section (continuous line) and the controller obtained from the third 


order model (broken lines) are presented in Figures 3.35, 3.36 and 3.37. 


Theatinali gains used for each) controller in Figuress3.35.,03-sorandss.3/ 
are approximately 10% (controllers FD0310, NADY0510), 20% controllers 
FD0320, NADY0520) and 30% (controllers FD0330, NADY0530) respectively of 
the corresponding ultimate gain of each loop. As expected, in every 
case the dynamic regulator gives slightly better results than the 
proportional controller designed based on the third order model. 

Figure 3.38 presents the simulated response of the fifth 
order model for a 10% step change in the product concentration, Co» set- 
point when the dynamic compensator obtained in this section (continuous 


line) and the constant compensator obtained from the third order model 
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(broken lines) are used. The final gains used for each loop in this 
case are approximately 20% of their corresponding ultimate gain. It can 
be observed that the dynamic compensator produces a less interacting 
control system than the controller obtained from the third order model. 
Since similar problems to those described in the previous 

section were expected in the experimental implementation of the dynamic 
regulator no decision about the final proportional and proportional- 
integral constants was made at this point. The selection of the final 
gains was done by experimental tuning and the results are presented in 


Chapter Six. 


3.8 CONCLUSIONS 

From the examples presented in the chapter and from the 
experience obtained in implementing and applying this method the follow- 
ing conclusions were drawn: 

1) The inverse Nyquist array is a simple technique which provides 
a systematic means of incorporating several typical control system 
design requirements, e.g. stability, interaction, integrity against 
transducer and error monitoring failures, and a satisfactory transient 
response. 

2) A direct and easy interpretation of one of the theorems due to 
Ostrowski [55] is possible when using this method. This theorem permits 
the determination of the location of the inverse of a transfer function 
of a loop when it is open and the rest are closed [67]. With this theorem 
it is possible to determine the exact stable gain space more accurately 
[69] and to find out when the system is under tight control. 


3) The requirement of diagonal dominance jis the most severe short- 
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coming of this method. Some design freedom could be wasted in making the 
system diagonally dominant. (Conversely, changes to the design method 
that would generate a narrower band than the Gershgorin band would mean 
greater accuracy and, in many cases, less control action to produce 
diagonal dominance.) In some cases, like in the fifth order evaporator 
model, it is not possible to fulfill this requirement with a proportional 
or proportional plus integral controller and consequently it is not 
possible to design a simple controller for this type of system. In other 
cases this problem can be overcome at the cost of additional calculations 
using the procedure indicated in Section 3.3. 

4) Since it is the controller matrix and not its inverse which 
will be implemented there is no direct indication in the inverse Nyquist 
array method of the complexity of the controller when a dynamic compen- 
Sator is being designed. This problem is more critical for discrete 
systems where as a matter of convenience it is desirable to use a 
bilinear transformation. 

5) The method does not have the flexibility to design a control 


system of high integrity for actuator failures. 


67 


CHAPTER FOUR 
THE CHARACTERISTIC LOCUS METHOD 


4.1 INTRODUCTION 

Most of the frequency domain design methods for multivariable 
control systems which have been proposed recently [e.g. 41,45,57,58] are 
based on a generalization of the scalar return-difference quantity 
introduced by Bode [7]. Included in this group is the characteristic 
locus method [4,5,41,42] which was developed by Belletrutti and 
MacFarlane [4] and has been discussed in detail by MacFarlane and 
Belletrutti [5,41,42]. 

Because of its more recent development the experience with this 
metnod is more limited than with the inverse Nyquist array method [65]. 
Only a very limited number of examples using this method have been 
presented in the literature [e.g. 4,42]. However this method is 
receiving more attention recently [60,72,76]. 

In order to compare the characteristic locus method with the 
inverse and direct Nyquist array methods the former method has also been 
included in this study. The investigation of the characteristic locus 
method has covered the same aspects covered in the study of the inverse 
Nyquist array method presented in Chapter Three. 

Since the characteristic locus method for continuous systems 
has been discussed in detail by MacFarlane and Belletrutti [4,5,41,42] the 
discussion of theoretical aspects of this method in Section 4.2 of this 
chapter is restricted to the design of discrete control systems. In 
Section 4.3 an example of the design of a multivariable regulator for a 


double-effect evaporator pilot plant represented by a discrete three- 
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state linear time-invariant model is presented. The results of the 
experimental implementation of the controller developed in this chapter 
are presented in Chapter Six. In the final section of this chapter some 
conclusions about the characteristic locus method based on the experience 


obtained in using this method are presented. 


4.2 THE CHARACTERISTIC LOCUS METHOD 
ree tteebdS 1 Ce NeOny 

The characteristic locus metnod was developed to design 
multivariable feedback control systems of the type represented by Figure 
2.1. For this type of system Popov [61] and later Hsu and Chen [24] have 


proven that 


2) ets =f (4.1) 


where: 
F(z) = bi + T(z) = return-difference transfer function matrix 


T(z) = return-ratio transfer function matrix 


CLCP = closed-loop characteristic polynomial 
OLCP = open-loop characteristic polynomial 
The return-ratio matrix T(z) will depend on the way the loops 
are broken in the control system represented by Figure 2.1]. When the 


loops are broken after the plant transfer function matrix, T(z),is equal 


to 
T(z) = G@(z) K(z) H(z) (4.2) 


If the loops are broken before or after the control matrix K(z) the 


return-ratio matrix T(z) is given by Wellas or ne respectively where 
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alle) = H(z) G(z) K(z) (4.3) 
and 
T,@) = K(z) H(z) G(z) (4.4) 


Since the feedback transducer matrix, H(z), can always be made 
part of the plant or control matrix without loss of generality it will be 
assumed that the matrix H(z) is equal to a unity matrix. In this case 
Equation (4.2) is equal to the open-loop transfer function matrix, Q(z). 

Using Equation (4.1) and the principle of the argument, 
Belletrutti and MacFarlane [4] have shown that, for systems represented 
by Figure 2.1, if the determinant of the return-difference matrix |F(z)| 
maps the contour D, shown in Figure 3.1, into Te encircling the origin Ne 
times clockwise, tnen the closed-loop multivariable system will be stable 


if and only if 
Ne = wl (4.5) 


where Po is the number of zeros outside the unit circle of the OLCP. 
The use of this stability criterion as a basis for a design 
method is impractical because it is not possible to design a multivari- 


able regulator in a systematic way, 1.e€., it is not possible to 


correlate the values of |F(z)| with the elements of the controller matrix. 


A more useful stability criterion has been obtained in terms 
of the characteristic values or the eigenvalues of the return-difference 
matrix, F(z). It has been shown [4] that if the characteristic values 
F.(z) (i=1,2 ... m) of F(z) map D into le, encircling the origin ne; 


times clockwise then 


70 


a a 
eee 


i 


a 
ete sfty Pot 
, ary b are S hw > 9% 
Tet ea ear -, al 
iL. 


m 
flys = RY Nes (4.6) 


and the system will be stable with all its loops closed if and only if: 
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Li 
is 


: tye (A573) 

Although this stability criterion is more useful than the 
previous one its use as a basis for a design technique is still very 
limited because it does not provide a systematic way to design the 
control matrix. 

A more practical method can be obtained if the characteristic 
values, t, (z) (i=1,2 ... m) of the return-ratio matrix T(z), of the 
system shown in Figure 2.1 are used. Since the return-difference matrix, 
F(z), and the return-ratio matrix, T(z), are related by the following 


equation: 
F(z) = iF eZ) (4.8) 


the relationship between the eigenvalues of F(z) and T(z), by applying 


the eigenvalue shift theorem [35], is given by: 
F.(z) =] + t.(z) Geb oon 4} (4.9) 


From Equation (4.9) it can be observed that the difference 
between the Nyquist diagram of fF. (z) (characteristic locus of F(z)) and 
the Nyquist diagram of t.(z) (characteristic locus of T(z)) is simply a 
shift in the imaginary axis. Consequently, if t, (z) maps the contour D, 
times 


shown in Figure 3.1, into f,. encircling the point (-1,0), 
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clockwise it follows that: 


Nes = Nye (4.10) 


and the multivariable system shown in Figure 2.1 will be stable with all 


its loops closed if and only if: 


m 
4 yap ie ; (4.11) 
Any of the return-ratio matrices given by Equations (4.2), (4.3) 
and (4.4) can be used to determine the stability of the closed-loop 
system using Equation (4.11). However the return-ratio matrix given by 
Equation (4.2) is usually preferred because the Humor matrix is 


equal to the open-loop transfer function matrix, Q(z). 


4.2.2 Design Objectives: Stability, integrity, non-interaction and 


accuracy 
Based on the stability criterion given by Equation (4.11) 


Belletrutti and MacFarlane [4] have developed the characteristic locus 
method. This is a general method which includes provision for taking the 
following design objectives into consideration: stability, integrity, 
interaction and accuracy. 

A control system is said to have good integrity if it remains 
stable under all combinations of a stipulated set of failure conditions 
[4]. A control system of high integrity, against the failure in the 
output transducer of one loop, can be designed using the characteristic 


locus method by ensuring that the characteristic loci (Nyquist diagrams 


of the characteristic values) of the principal sub-matrix obtained by 
deleting the corresponding column and row from the return-ratio matrix 
given by Equation (4.2) satisfy the stability criterion given by 
Equation (4.11). Other types of failures can also be considered with 
this method by requiring that the respective principal sub-matrix of the 
appropriate return-ratio matrix satisfy the stability criterion given by 
Equation (4.11) [4]. 

For systems represented by Figure 2.1 the closed-loop transfer 


function matrix, R(z),is expressed in the following form: 


R(z) = [1 + Q(z)]”!Q(z) (4.12) 


where: 


CZ Ze ec c Jaki) 


When the eigenvalues or cnaracteristic values of Q(z) are 
distinct and are associated with linearly independent characteristic 
vectors or directions denoted by wi (2), (j=1,2 ... m) the open-loop 


transfer function matrix can be expressed in the following way [42]: 


Q(z) = W(z)[diag{t;(z)}] V(z) (4.13) 

where: 
W(z) = [w,(z), Wo (z) aes w,(Z) J (4.14) 
V(z) = wl(z) (4.15) 


t,(z) = jth characteristic value of Q(z) or T(z) 
Then from Equation (4.12) it can be shown [42] that the closed- 


loop transfer function matrix is equal to: 


EeeZ) 
R(z) = W(z)[diag (F->E>qp)] Vz) (4.16) 
7 


It can be observed from Equation (4.16) that the characteristic 
values, r,(z); of the closed-loop transfer function matrix are given by: 


t,(z) 


vs (z) 2 eae e (2) alec oem aie 


and by comparing Equations (4.16) and (4.13) the characteristic 
directions (eigenvectors) of the closed-loop transfer function matrix 
are equal to the characteristic directions of the open-loop transfer 
function matrix, Q(z). 

Using this information the characteristic locus method can be 
used to design a control system of low interaction by aligning the 
eigenvectors of the open-loop or return-ratio matrix with the standard 
basis set of vectors, i.e., by designing Q(z) = G({z) K(z) such that the 
matrix W(z) approaches the identity matrix. 

At low frequencies interaction can be reduced when using the 
characteristic locus method [42] by ensuring that the gain of the 


characteristic loci, t.(z), of Q(z) are high enough such that: 


(t,@) is 
ie Bie) | lees. om (4.18) 


From Equation 4.16 it can be seen that both approaches wil] 
tend to make the closed-loop transfer function matrix, R(z), diagonal. 

The accuracy of the control system, defined as the degree to 
which the actual outputs follow the desired outputs [42], can also be 


considered when using the characteristic locus method. High accuracy can 
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be obtained in a system when the magnitudes of the characteristic loci of 
the open-loop transfer function matrix are relatively large at low 
frequencies. When this is the case Equation (4.18) will be satisfied and 
the closed-loop TFM will tend to be equal to a unity matrix. Consequently 


the output will follow closely any change in the setpoints. 


4.2.3 Design Procedure 


The characteristic locus method can be applied to square as 
well as non-square plants but it is relatively easier to apply to square 
plants. 

The design procedure using the characteristic locus method 
involves four phases whose order depends upon the relative importance of 
each property. These four phases are 

- stability phase 
- integrity phase 
- interaction phase 
- performance phase 

In each phase a basic controller, K-(z), is designed jn an 
interactive and iterative way to give the system the desired properties. 
The final controller is giyen by the product of the basic controller 
matrices found during the different phases, i.e. 

K(z) = K. (z) (4.19) 

Three types of diagrams are commonly used with this method. 

The most important ones are the characteristic loci diagrams (part a of 


Figures 4.1 through 4.15) which are the Nyquist diagrams of the character- 


= 
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istic values (eigenvalues) of the open-loop or return-ratio transfer 
function matrix, Q(z). These diagrams are used to determine graphically 
the stability of the closed-loop system by applying to them the stability 
criterion represented by Equation (4.11). Another diagram (part b of 
Figures 4.1 through 4.15) used is the magnitude of the characteristic loci 
as a function of frequency. This plot is useful to select the gains which 
are desirable in order to have low interaction and high accuracy at low 
frequencies. From this diagram it can be readily observed when Equation 
(4.18) will be satisfied and consequently when the system will have non- 
interacting and high accuracy properties. The third diagram (part c of 
Figures 4.1 through 4.15) plots the alignment of the characteristic 
directions (eigenvectors) of the matrix Q(z) or R(z) with the standard 
basis yectors as a function of frequency. These diagrams are useful to 
determine the degree of interaction of the closed-loop and open-loop 
System at any frequency. When the characteristic directions are properly 
aligned they will be orthogonal, i.e.,matrices W(z) and R(z) in Equation 
(4.16) will be diagonal. Consequently the closed-loop system will be 


non-interacting. 


4.2.4 General Comments 

Different nee of basic controller matrices haye been suggested 
[5,42] to give the control system specific properties. Several of these 
controllers have also been recommended for use with the inverse Nyquist 
array method [65,67]. The idea of the commutative 
controller can be incorporated as part of this method. The commutative 
controller method [36] can be used when it is required to change the 


characteristic values without changing the characteristic directions 
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(eigenvectors) of the open-loop transfer function matrix. 

Unfortunately the design of a controller using the character- 
istic locus method, at least as originally proposed, is a difficult task 
especially for a plant whose transfer function matrix is of order greater 
tnan two. In higher order systems it is difficult to establish the 
relationship between the characteristic loci and the elements of the 
controller matrix and hence the systematic selection of a controller 
becomes the major problem. The designer also does not get any physical 
insight into the system because the characteristic loci do not have a 
"pnhysical meaning" such as can be attributed to the direct Nyquist 
diagrams. 

Recent papers have proposed more systematic approaches to the 
to the design of controllers using the characteristic locus method [60, 
508] One approach proposed by Owens [60] involves the selection of a 
controller such that the characteristic values of the open-lcop matrix, 
Q(z), take specified values at a given frequency. This is done by expanding 
the transfer function matrix in a dyadic form [59]. The other scheme 
developed by Kouvaritakis [508] designs the controller by breaking the 
frequency range of interest into three regions: low, intermediate and 
high. At high frequencies a real controller, K, > is designed to improve 
the aligment between the eigenvectors of G(z) and the standard basis 
vectors. The resulting open-loop transfer function matrix Q, (2) = G(z)K,, 
is then considered as the new plant transfer function matrix. At inter- 


mediate frequencies a new controller, K = W diag(k;(z)) V. is designed 


m 
to compensate (about the critical point) the characteristic values of 


Q, (2) (where W, and V., approximate at intermediate frequencies the 
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eigenvector and the reciprocal eigenvector structure of Q, (2). The 
controller See must tend to a unity matrix at high frequencies so that 
it does not upset the eigenvector alignment affected in the first stage 
of the design. The new open-loop transfer function matrix Q5(z)=G(z)K,K (z) 
now becomes the new plant transfer function matrix. Finally at low 
frequencies a proportional-integral controller, Ky is designed. The 
proportional matrix will be a unity matrix and the integral matrix will 
be equal to oWjdiag(k,)V5 where W and Vy approximate, at low frequencies, 
the eigenvector and the reciprocal eigenvector structure of Q5(z). The 
gains, k., are selected to adjust the magnitude of the characteristic 
values of Q, (Z) and a determines the duration of the integral action as 
the frequency increases. This controller must also tend to the unity 
matrix at high frequencies. 

The design of multivariable control systems using the 
characteristic values of the return-difference and return-ratio matrix 
to determine the stability of the system was proposed by Bohn [502 ,503] 
previous to the development of the characteristic locus method by 
Belletrutti and MacFarlane [4]. In his work [503] which was mainly 
directed toward the diagonalization of the return-ratio matrix, Bohn 
pointed out some of the difficulties of working with characteristic values. 
These difficulties have also been commented on more recently by Rosenbrock 
and Cook [72] and by MacFarlane [506]. The difficulties arise because 
the characteristic values (eigenvalues) of the return-difference matrix 
are not necessarily rational polynomials as in the case when the closed-loop 
characteristic polynomial is factored by conventional methods. Therefore 


the characteristic values may not be single-valued analytical functions 
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of z at all point within the closed contour D. If the characteristic 
values are not single-valued then the Nyquist criterion cannot be applied 
directly to each individual characteristic value unless branch cuts 

are introduced in the Nyquist contour [503]. These branch cuts are 
important in determinig the number of encirclements of the critical point 
of each individual characteristic value. However, they do not need to be 
considered when determining the total number of encirclements, Ne = : n 


i=l 
of the deteminant of the return-difference matrix because the branch cut 


ti’ 


locus is always traced through in opposite directions for each conjugate 
pair of characteristic values [503]. Thus, in practical applications, 

the test for system stability based on the total number of encirclements, 
ne, is not changed by dia (paasnnel existence of irrational characteristic 
values. 

One important restriction in the characteristic locus method 
arises because the characteristic values are a property of the system 
matrix rather than being associated with individual feedback gains 
{k,,7=1,..m}. Therefore, the use of characteristic values to determine 
the stability region of a system requires that the same gain be used 
for each loop in the system [503]. Recently Rosenbrock and Cook [72] 
have examined the characteristic locus method for the case when different 
gains, ks are used in each loop. It has been found [72] that if different 
gains, {k,,i=1,..m}, are used then further restrictions are required 
on the return-ratio or open-loop transfer function matrix, More 
specifically, when different feedback gains tk, ,i=1,..m} are used the 


characteristic locus method is valid only when: 
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1. The open-loop matrix Q(z) is a normal matrix for all values of 


z around the contour D, i.e. 
Q(z) Q'(z) = al(z) Q(z) (4.20) 


where g!(z) is the transpose of Q(z). Then the stability of the 
closed-loop system can be determined by application of a 
theorem due to Freeman et al. [18]. 

2. The return-difference matrix, F(z), is diagonally dominant. The 
stability of the system in this case can be determined using a 
theorem due to Rosenbrock [70] and Cook [11] ( This theorem 
is the basis of the direct Nyquist array method discussed in 
next chapter.). 

The same restrictions apply for the case discussed in Chapter 
Three when the characteristic loci of the inverse of the open-loop trans- 


fer function matrix are used to design a control system. 


4.3. DESIGN OF A MULTIVARIABLE REGULATOR FOR THE THIRD ORDER EVAPORATOR MODEL 
A design of a multivariable regulator for the double-effect 

evaporator, described in Appendix A, using the characteristic locus 

method is developed in this section. The model] used in this example js 

a discrete, three state, linear, time-invariant model with normalized 

perturbation variables. This model has been obtained from a tenth order 


model using the discrete Marshall's reduction technique [78]. The tenth 
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order model has been obtained by Wilson [77]. This same model has also 
been used in Section 3.6 to design a multivariable regulator using the 
tnverse Nyquist array method [65]. 


The z-transform of the model is given by: 


W (z) S(z) 
W (Z) = | G(z) B, (z) CA) 
C,(z) B,(z) 


where G(z) is the plant z-transfer function matrix and is shown in Table 
3.1. The outputs of the double-effect evaporator are the first effect 
hold up, Ws the second effect hold-up, Wo s and the product concentration, 
C., The inputs of the system are the steam flow to the first effect, 
S, the bottoms flow from the first effect, B,, and the bottoms flow from 
the second effect Bo. 

For reasons explained in Section 3.4 the design of the 
multivariable regulator for the double-effect evaporator has been carried 
out using a w-bilinear transformation of the z-plant transfer function 
matrix. Thus G(w) has been used instead of G(z), where w = ae The 
w-transfer function matrix, G(w) is presented in Table 3.2. The poles 
of this transfer function matrix are shown in Table 3.1. 

The system is marginally stable because two of its poles are 
at the origin in the w-plane. To determine the stability of the system 
the Nyquist contour shown in Figure 3.2 has been used. Using this 


contour the number of right-half-plane zeros of the open loop character- 


istic polynomial is equal to zero, i.e. 


Py = 0 (4.22) 
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The design of the controller is based on the characteristic 
loci diagrams of the open-loop or return-ratio transfer function matric 
The three characteristic loci of G(w) are shown in Figures 4.1, 4.2 and 
4.3. Each figure consists of three graphs. Graph (a) shows a character- 
istic locus of G(w). The labels indicate the numerical value at the end 
of each axis. The magnitude of this characteristic locus for different 
frequencies is shown in graph (b). The magnitudes are given in decibels 
and the frequency in radian/sec. Finally the alignment of the corre- 
sponding eigenvector with a standard basis vector at different 
frequencies is presented in graph (c). In this graph the angles are 
given in degrees. From these figures it can be observed that the closed- 
loop system without compensator will be unstable since they do not 
satisfy the stability criterion given by Equation (4.11). It also shows 
that there is a considerable amount of interaction in the open-loop 
transfer function matrix because eigenvectors 2 and 3 are both aligned 
with the same standard basis vector. 

The most important design factor is the absolute stability 
followed by the interaction of the system. The performance 
and integrity of the system are considered later in the design. 

The characteristic loci diagrams in Figures 4.1, 4.2, and 4.3 
provide an indication of how well the control system will perform but 
they do not give any specific idea about the type of controller which 
must be used. 

The selection of the controller is usually done by trying a 
specific type of controller matrix which is known to give the system 
certain properties. Then a decision about the value of this type of 


control is made by observing the characteristic loci diagrams. 
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A controller to reorder the inputs of the system was tried 
first. This could reduce the interaction and improve the relative 
Stability of the system. From previous experience with the double-effect 
evaporator the best control scheme has been found to be when first effect 
hold ups We is controlled by bottoms flow By the second effect hold up, Wos 
is controlled by bottoms flow, B, and the product concentration, Cas is 


controlled by steam flow, S. This can be accomplished by using the following 


control matrix: 


0 0 a0 
K = Ae) 0 (4.23) 
0 Wo 
The new OLTFM is then given by 
Q, (w) = G(w) K, (4.24) 


The effect of Ky on the characteristic loci of Q, (w) RSedihalcuLcecd 
predict but can be observed in Figures 4.4, 4.5 and 4.6 which show the 
three characteristic loci of Q, (w). Although two eigenvectors of Q, (w) 
are aligned with the same standard vector the system is less interacting 
in this case as shown by Figures 4.4c, 4.5c and 4.6c. The second 
eigenvector is perfectly aligned with the second standard vector as shown 
by Figure 4.5c. The siete is however unstable because it does not 
satisfy the stability criterion given by Equation (4.11). It is necessary 
to change the signs in the characteristic loci 1 and 2 to make the 


system stable. This is done by using the following control matrix 


iS 0 =| 0) (4.25) 


From the characteristic. loci of Q (w) = Q, (w) K, in Figures 


4.7, 4.8 and 4.9 it can be observed that the closed-loop system is stable. 


This suggests that the characteristic loci 1 and 2 are strongly 
influenced by columns 1 and 2 respectively. In this case the three 
eigenvectors of Qo (w) are aligned with the second standard basis vector 
which indicates that Q, (w) is more interacting than Q, (w). 

The closed-loop system is already stable. However the system 
does not have the desired properties, i.e., interaction, integrity, 
performance with the controller matrix Kiko. Therefore the controller 
matrix Ky K, was eliminated and another control matrix was tested instead. 

One type of controller which is very useful and very often 
recommended [5,42,67] is one which makes the open loop system diagonal 
at very low and/or at very high frequencies. At very low frequencies the 


recommended [42] controller is given by: 
K, = & '(0) (4.26) 


For the double-effect evaporator 7! (0) does not exist because 
it has two poles at the origin. However the controller which makes the 


system diagonal at high frequencies [42] is equal to: 


K = G | («) 
or 
EO WH Sach 
K = | 18.508 0.0 11.650 (4.27) 


24.518 49.150 49.225 


In order to apply this controller to the plant jt is necessary 
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85 
to adjust the gains in each column. Then the control matrix will be 
given by 
Ky = KD, (4.28) 


where D, is a diagonal matrix that can be chosen arbitrarily. In this 


case D, was chosen to be equal to: 


1/24.5 0 0 
Dy = 0 1/49.15 0 (4.29) 
0 Oem /2i-5 


so that the controller designed in this section is equal to 
the controller obtained using the direct Nyquist array in 


Section 5.4, 


Then the new open-loop transfer function matrix is equal to 
1 Fe ] 
Q, Cw) = GQw)Ky (4.30) 


The characteristic loci of a! (w) are presented in Figures 4.10, 4.11 and 
4.12. From these figures it can be observed that the system is almost 
diagonal (the eigenvectors are properly aligned and consequently 
orthogonal) not only at very high frequencies but it also is diagonal for 
all the frequency spectrum. For this reason the characteristic loci can 
be used without any problem. From these diagrams it can also be 

observed that the system is unstable. This can be easily corrected by 


using the following controller: 


4) a Gr ee een (4.31) 
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Figures 4.13, 4.14 and 4.15 show the characteristic loci of 
the new open-loop transfer function matrix Q|(w) = OF Kee The 


closed-loop system with the controller Kp = eel, Tees 


=O 030 0.0 1.000 
Kp S| #017705 0: 0 -0.420 (4.32) 


-1.000 =56 -].777 


has very good properties besides stability. It is almost non-interacting 
and has high integrity against different transducer and error-monitoring 
failures. 

The final phase in the design is the performance phase. It is 
a requirement in the design that the system must possess good performance 
characteristics. Good accuracy (in the sense of following changes in the 
setpoints) is achieved for setpoint changes whose frequency content is in 
the range in which the moduli of all the characteristic loci of the open- 
loop transfer function matrix are large. These moduli can be increased by 
increasing the final gain for each column. 

From the characteristic loci of Q, | (w) in Figures 4.10a, 4.1]a 
and 4.12a it can be observed that the system can be represented by three 
first-order sub-systems. Thus the offsets of the system can be reduced 
by increasing the moduli of the characteristic values in Figure 4.10b, 
4.11b and 4.12b. This means that proportional control only jis sufficient 


to have a good control system. In this case the final control will be 
K = Keke (4.33) 


where, Kp is given by Equation (4.32) and Ko is equal to 
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kK, 0 0 
Ke = (pee ko 0 (4.34) 
0 0 k 


Before selecting the final gains in each loop it is convenient 
to determine the stability region or the stability margins of the closed- 
loop system. In general this is done by selecting the controller Ko Such 


that: 


5 (4.35) 


where I, is the identity matrix and k is a parameter which is used to 
Study the stability of the system. From Figures 4.13a, 4.14a and 4.15a 
it is found that the system will remain stable for the following values 


Otek 
Pare SOIT = 24.5 (4.36) 


In this particular example it is possible to expand this 
stability region by using different gains in each loop, i.e. by using 
Equation (4.34). This is because the system has very good diagonal 
dominance (almost diagonal) as will be seen in the next chapter ( cf. 
Figure 5.12) and the Gershgorin bands are disjoint [72]. This ensures 
that the characteristic values of Q5 (w) are single-valued for every 
value of w on the contour D and that a different gain can be applied 
to each characteristic value [72]. The expanded stability region for 
which the closed-loop system will remain stable can also be obtained 


from Figures 4.13a, 4.14a and 4.15a and it given by: 


NS hp ae (4.37) 


1 0.0407 
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Once the stability region has been determined the design 
procedure can be continued by selecting the final gain in each loop. 

The selection, in this case, was also done based on the stability margin 
of each loop. However for reasons similar.to the ones explained’ in Section 
3.6 the selection of the final gains was done by experimental tuning. 

By selecting ky» k, and k, in Equation (4.34) to be 50%, 30%, 
20% and 10% of their respective ultimate gain the controllers FD0350, 
FD0330, FDO320 and FDO310 were obtained respectively. These control 
matrices are almost equal to the ones obtained in Section 3.6 using 
the inverse Nyquist array method. The simulated responses of the 
evaporator when these controllers are used and a 20 % step change in the 
feed flow is introduced in the system are presented in Figure 3.11. The 
simulated responses of the evaporator for a 10% step change in the setpoint 
of the product concentration when controller FD0350 and FD0320 are used 
are shown in Figure 3.13. 

Integral action can also be added to the control system by 
designing, in this case, the matrix Ke given by Equation (3.35). The 
selection of this controller was also done by experimental tuning. 

The experimental results obtained by implementing the controllers 


designed in this section are presented in Chapter Six. 


4.4 CONCLUSIONS 


The characteristic locus method [4,42] is a technique which can 
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be used to design a stable control system with good integrity properties, 
low interaction and good performance. This method is more general than 
the inverse Nyquist array method [65] since it does not require the 
System to be diagonally dominant. 

On the other hand the characteristic locus method is a less 
practical method than the inverse Nyquist array technique. The design of 
Simple control matrices with the characteristic locus method is very 
difficult because the characteristic loci show how good the control is 
but they do not give any specific indication about the way the controller 
must be modified to improve the control system. This means that the 
design of the controller with this method cannot always be done in a 
systematic way. However the characteristic loci diagrams are a very good 
tool in the analysis of a control system. A controller can be easily 
accepted or rejected according to inspection of the characteristic loci 
diagrams of the open-loop transfer function. Nevertheless it is also 
possible to design a control system in a systematic way with the 
characteristic locus method if a scheme recently developed by 


kouvaritakis [508] is used. 
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CHAPTER FIVE 
THE DIRECT NYQUIST ARRAY METHOD 


5.1 INTRODUCTION 

A very practical and powerful design technique which has been 
suggested [38,68,70,71] but not fully developed or discussed in the litera- 
ture is what can be called the "direct" Nyquist array method. This method 
arises from the combination of principles from the inverse Nyquist array 
[65] and the characteristic locus methods [4,42]. More specifically, it 
combines the requirement of diagonal dominance used in the inverse 
Nyquist array method [65] with a stability criterion based on the 
encirclements generated by the return-difference transfer function 
ae used in the characteristic locus method [4,42]. The direct 
Nyquist array method can also be derived by intuitive extension [15] of 
the conventional single-input, single-output Nyquist design procedures 
to the design multivariable control systems. 

The direct Nyquist array method has been evaluated in this 
study and the following discussion focuses attention on the development 
and applications of this method. 

The following discussion of the direct Nyquist array method is 
restricted to discrete systems. However, this discussion also applies 
directly to continuous systems. The basic theoretical aspects of the 
direct Nyquist array method are presented in Section 5.2 and they build 
upon what has been said about the characteristic locus method in Chapter 
4. In Section 5.3 the development of the direct Nyquist array method 
by extension of conventional techniques is presented as a self-contained 


section. The application of the direct Nyquist array method to the 
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design of a control system based on discrete three- and five-state, 
linear, time invariant double-effect evaporator models is presented in 
Sections 5.4 and 5.5, respectively. The models used in these examples 
are the same as the ones used in Chapters Three and Four. 

In the final section of this chapter some conclusions about 
this method are presented based on the experience obtained in applying 
this method and from comparing this technique with the inverse Nyquist 


array method and the characteristic locus method. 


5.2 THE DIRECT NYQUIST ARRAY METHOD 
eee leotabi lity Criterion 

The characteristic locus method which has been discussed in 
Chapter Four is based on the Nyquist criterion of encirclement of the 
Origin by the characteristic values of the return-difference matrix, 
F(z). It has been shown [4,42] that if the characteristic values 
{f.(z), i=l, ... m}-of the return-difference transfer function matrix 
of a system represented by Figure 2.1 map the contour D, shown in 
FAGUregse|, hinto Te encircling the origin Ne times clockwise then the 
multivariable system will be stable with all its loops closed if and 


only if 


m 


e Ne; = —Po (Sv 
;= 


] 


where p. is the number of zeros outside the unit circle of the open-loop 
fe) 
characteristic polynomial. 
It has also been mentioned in Chapter Four that this stability 


criterion was not very practical because it did not allow the design of 
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the controller to proceed in a systematic way. However if the concept 
of diagonal dominance used in the inverse Nyquist array method [65] is 
used then a more practical stability criterion can be obtained. By 


requiring the return-difference matrix, F(z) to be diagonally dominant, 


i.e. that 
m 
|f.5(z)| it If 55 (2)| (5.2) 
ake 
for column dominance or 
m 
egal > E 1Fg 2) Ga) 
j#i 


for row dominance, where ee represents element (i,j) of F(z), a 
stability criterion, which involves only the diagonal elements, F..(z), 
of F(z) can be obtained. Using a theorem due to Gershgorin and a proof 
Similar to the one used by Rosenbrock [65] in the inverse Nyquist array 
method it can be shown [11,38,70] that if the return-difference matrix 
is diagonally dominant :for all the z on the contour D and the diagonal 


! <a : le: ae ibe 
elements f.- (2) map D into Vee encircling the origin Nes, times cloc 


wise, then the system will be stable with all its loops closed if 


m 


Neiq = “Po (5.4) 


i= 


Since the return-difference matrix, when the transducer feed- 


back matrix is an identity matrix, is equal to 


MBS Wea nlea (5.5) 
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the relationship between the elements of F(z) and the open-loop 


transfer function matrix, Q(z),is given by: 


f,4(z) = 1+ 4,;(z) (5.6) 


and 
£552) i q; 5 (2) nN Ge al (OF) 


From Equation (5.6) it can be observed that the difference 
between the Nyquist diagram of f;.(Z) and q;; (Z) is a shift in the 
imaginary axis. Thus using Equation (5.4) it can be shown [38,70,71] that, 
when the return-difference matrix of a system represented by Figure 2.1 


is diagonally dominant for all the z on the contour D, the closed-loop 


system will be stable if and only if the number of encirclements of the 
critical point (-1,0) by the Nyquist diagrams of the diagonal elements, 
{q.,,i=1,..m}, of Q(z) is equal to Epa where p, is the number of poles 


of Q(z) outside the unit circle in the complex plane. The direct Nyquist 


array method is based on this stability criterion. 


5.2.2 Design Procedure 


The design procedure used in the direct Nyquist array method 
is very similar to the one used in the inverse Nyquist array method. 

The direct Nyquist array method makes use of the Nyquist 
array concept introduced by Rosenbrock [65] in the development of the 
inverse Nyquist array method. The Nyquist array of Q(z), which isa mxm 
set of Nyquist diagrams corresponding to each element in the matrix 
Q(z), is a very useful graphical tool for determining the effect of the 


controller matrix onthe controlled system and for selecting the type 
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of control which is necessary to satisfy the design requirements. 
Diagonal dominance of the return-difference matrix can also 

be determined graphically by drawing the Gershgorin circles with their 

;; (2), 


OfR(Z). amiinesradia d. of the Gershgorin circles are equal to the sum 


centers located on the Nyquist loci of the diagonal elements, q 


of the magnitudes of the off-diagonal elements, i.e. 
dj = 2 [4;;(2)| (5.8) 


if column dominance is used. When row dominance is used the radii of 


the Gershgorin circles, are equal to 
m 

d. : } 19, 5 (2) egd)) 
if 


The return-difference matrix ts dtagonally dominant tf the 
Gershgorin ctrcles do not tnelude the point (-1,0). 

Rosenbrock [69] has shown that these Gershgorin bands have the 
property of containing a narrower band, the Ostrowski band, which is the 
boundary of a region where the transfer function between input 7 and 
output i is located, when loop 7 is open and the rest are closed. The 
Ostrowski circles have their centers on the Nyquist loci of {q;; (2), 
1a se ewem sand their radii are equal to the radius of the Gershgorin 
circles multiplied by the factor as where 

d. 


= J ; = 
O ae T+ 4557 J lees 0) (5.10) 
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The design of a control system using the direct Nyquist array 


method involves three steps: 


1) Pairing the input-output plant variables. 

2) Design of a simple pre-compensator which makes the return- 
difference matrix, F(z), diagonally dominant and gives Q(z) 
the desired properties. 

3) Design of a multiloop control system based on the diagonal 


elements of the open-loop transfer function matrix, with their 
respective Gershgorin and Ostrowski bands. 


9.2.3 Advantages of the Direct Nyquist Array Method 


In this three step procedure the designer can take into account 
in an interactive and iterative way the following design requirements: 
Stability,interaction, integrity against different failures, accuracy 
and the transient response of the closed-loop system. The design of the 
controller using this method is done in a systematic way and at each 
step the designer has some indication as to how the controller should be 
modified to achieve the desired performance. 

The main difference between the direct and the inverse Nyquist 
array methods is that the direct method is based on the "conventional" 
Nyquist diagrams of Q(z) with which the practicing control engineer is 
familiar, while the inverse Method is based on the Nyquist diagrams of 
gl (z). Another major difference is that the inverse method requires 
the inverses of the closed- and open-loop transfer function matrices to 
be diagonally dominant while the direct method requires only the return- 
difference matrix to be diagonally dominant. This means that in the 
inverse method in order to determine the stability of the system the 


Gershgorin circles must not include neither the origin nor the critical 
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point (-1,0) while in the direct method it is only necessary that the 
Gershgorin circles do not enclose the critical point (-1,0). As can be 
seen there is more flexibility and freedom in the selection of the 
controller with the "direct" approach. When Q(z) is known algebraically 
it is also possible to design a control system with the inverse Nyquist 
array method without requiring diagonal dominance from the inverse of 
the open-loop transfer function matrix. However this. requires. some 
additional calculations. 

It is also possible to derive the direct Nyquist array design 
procedure based on an intuitive extension of conventional single-input, 
Single-output Nyquist design procedures. Such an approach has the 
advantage of requiring no new theory and all the concepts and mathemati- 
cal operations involyed are familiar to most control engineers. 
Furtnermore the intuitive development provides greater insight into the 
physical meaning of each step in the design procedure and also results 
in a more practical and more powerful design technique than has been 


obtained to date. This development is presented in the next section. 


5.3 DEVELOPMENT OF THE DIRECT NYQUIST ARRAY METHOD BY EXTENSION OF 

CONVENTIONAL DESIGN TECHNIQUES 

The problem of interest is feedback control of a multi-input 

multi-output, linear system such as shown in Figure 2.1]. It is assumed 
that the system of interest is fully described by the mxg transfer 
function matrix, G(z), and that the compensator matrix, K(z), is to be 
designed to give stable, closed-loop control. In order to emphasize 
the two main steps in the proposed design procedure the feedback system 


will be treated as shown in Figure 5.] in which the compensator is 
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divided into two parts: 

1) the dynamic component K(z) which is designed first so that the 
augmented plant Q(z) = G(z)K(z) is square and has the desired 
dynamic characteristics. 

2) the static component, K, which is a diagonal matrix of constant 
controller gains {k.., i=] ... m}. AS shown later, K can be 
designed on the basis of approximations which bound the true 
value of the system transfer functions and permit a modified 
Nyquist stability analysis to be completed. 

The final step in the proposed design procedure generates m conventional 
Nyquist diagrams each of which is an exact representation of the 
relationship between one input-output pair of the augmented system when 


all the other loops are closed. The necessary and sufficient condition 


for stability of the closed-loop multivariable system, when the return- 

difference matrix is diagonally dominant for all the z on the contour D, is 

Simply that these diagrams all satisfy the Nyquist criterion (cf. Appendix D). 
For the sake of completeness the following discussion includes 

some concepts and definitions that have been discussed elsewhere or are 


already familiar to most control engineers. 


5.3.1 Formulation of the Design Problem 


System Model 


Multivariable frequency domain design techniques assume that 
a suitable mathematical model of the system of interest can be deriyed 


and expressed in the form: 


y(z) = G(z)u(z) (Sal) 
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Equation (5.11) shows that the relationship between the m-dimensional 
vector of plant output variables, y(z), and the 2-dimentional vector of 
plant input variables, u(z) is defined by the system. transfer function 
matrix, G(z). It is written in z transform notation because multivari- 
able control systems are normally implemented using digital computers 
and hence the discrete formulation is most appropriate. However, the 
following development also applies directly to continuous systems which 
can be represented by Equation (5.11) written as a function of the 


Laplace transform variable s. 


Control Strategy 
It is assumed that the system described by Equation (5.11) is 


to be controlled by the multivariable feedback control system shown in 


Figure 2.1 where: 


Wa) = Tae) Se aa (5.12) 


K(z) = gxm compensator matrix 


r(z) = an mxl vector of setpoints or inputs. 


It is easily shown that the controller matrix K(z) can be expanded, 


without loss of generality into the form 
K(ziy= iK(z) Ki (5.13) 


where: 


K is a mxm diagonal matrix of constant controller gains, ks 


K(z) is a 2xm matrix of dynamic elements 
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It will be seen that the use of a general controller matrix K(z) as 
defined in Equation (5.13): 

| ] is frequently necessary to make the design procedure practical. 

2) provides flexibility in the pairing of input-output variables. 


) 
) 
3) permits non-square systems, G(z), to be handled. 
) permits the specification of dynamic feedback controllers. 
) 


5) provides a means of achieving specific design objectives such 
as non-interaction, or specified system dynamics. 

6) allows the designer to consider properties that are unique to 
multivariable systems such as the “integrity" (e.g. stability 
of each loop when one or more components in the other loops 
Palle 

7) implies no loss of simplicity since setting K(z) equal to the 
identity matrix is equivalent to working directly with G(z). 

Figure 5.1 shows the general feedback control system of Figure 2.1 with 
the controller K(z) broken into two parts: the static component, K, 

and the dynamic component, K(z). The following discussion will be 
concerned first with determining the elements kK. of K and hence it will 
be assumed that the design of the dynamic part of the controller is 
complete and that it has been included as part of the augmented system 


Q(z) where 
Q(z) = G(z) K(z) (5.14) 


The specification of the dynamic compensator, K(z), such that the 
augmented system Q(z) has the desired properties is discussed later. 
(Note that in Section 5.3 q,5 (2) refers to an element of Q(z) while in 


the rest of the thesis it refers to an element of Q(z) = G(z)K(z) K = 
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Q(z)K. The use of Q(z) rather than Q(z) emphasizes the role of the 
feedback gain elements Kas of the matrix K but in latter sections it will 
be seen that the usual practice is to use Q(z) to define the augmented 
plant at all stages in the design procedure regardless of how K(z) is 
defined. The meaning is normally clear from context. ) 

The above approach implies that although G(z) and K(z) may be 
non-square, Q(z) is always square and that the diagonal 
elements of the K matrix can be interpreted as constant feedback 
proportional controller gains, Kiss associated with the control of the 
input/output pair u;(z):y;(z). One of the advantages of this approach 
is that it associates a single gain element, Ks with each output, Yas 
which is convenient both in the design stage and for field tuning of 
the resulting control system. 

The specification of the individual controller gains, Kia 
obviously depends on the transfer functions relating each input-output 
pair of the augmented system. Therefore the following section discusses 


the derivation of these transfer functions. 


5.3.2 Development of the Transfer Function Relating Each Input-Output 
Bot ’ 
The "open-loop" relationship between the input to the augmented 
plant, u;,and output Yi» when the proportional feedback control loops 


are closed around al] the other input-output pairs, is defined as: 
y;(z) = hy (z) us (2) (5.15) 


This relationship is illustrated in Figure 5.2 where the large shaded 


block represents the augmented system Q(z) with inputs u: and outputs 


10] 


y;(z). Obviously h,(z) must be derived from the OLTFM of the augmented 
plant and the following subsections consider four cases in order of 


‘increasing generality (and complexity). 


Diagonal OLTFM 
The simplest case is when Q(z) is diagonal since h,(z) is 
then equal to the diagonal element q.;(z) and the selection of the 
feedback controller gains reduces to m independent, single-input single- 
output design problems. The design can be completed using any 
applicable design approach such as conventional Nyquist stability 


analysis. 


Ze 2a OLEM 

In the more general case where the off-diagonal elements of 
Q(z) are non-zero there are interactions between the components of the 
input and output vectors, and in general h, (Z) in Equation (5.15) 
becomes a function of several elements of Q(z) and also of the feedback 
controller gains, K-45 1n all the other loops. This is readily apparent 
when the general feedback system shown in Figures 5.1 and 5.2 is expanded 
into the more familiar block diagram form. Figure 5.3 shows the block 
diagram representation of a 2 x 2 system of the type shown in Figure 5.1] 
from which it follows directly that the transfer function relating the 
second input-output pair (assuming that the other loop is closed) is 
given by: 
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which illustrates the dependence or interaction, between loop 2 and loop 
1. This approach to the analysis of a 2 x 2 system is developed more 


fully in a separate research report [30]. 


Saxe oe Ure 
The analysis of higher order systems is considerably more 
complicated but a straightforward block diagram analysis of a 3 x 3 


system will show that hy (Zz) WomaSeshownein EGuationa(>. 17). 


m x m OLTFM 

The transfer functions defined by Equation (5.16) and (5.17) 
consist of the appropriate diagonal element q44 (2) plus m-] terms which 
represent the effect of interactions with the other feedback control 
loops. Each of these "interaction terms" has the form of a coefficient 
multiplying an off-diagonal element in column i of Q(z). This 
representation generalizes to the form given by Equation (5.18). The 
relationship given by Equation (5.18) was derived independently and has 
a different structure, but gives the same h.(z) (when the men loop is 
open and the other m-1 loops are closed) as given by Rosenbrock [60]. 
Once the relationship between each input-output pair has been determined 


it is possible to proceed with the specification of the gains, Kk... 
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goo! (2) _ cofactor of element (iy iQeor F(z) 
cofactor of element (i,i) of F" (z) 


6" (2) _ cofactor Gj)" of Grice 


ij cofactor (ii) sof lew! 
F(z) . = Return Difference Matrix = (1 + Q(z)) 
F" (z) = Normalized Return Difference matrix obtained 


by dividing each column/row of F(z) by the 


dtagona] element in that column/row 


(5212) 


(5.18a) 


(5.18b) 


(5.18c) 


(5.18d) 


(5.18e) 


(5. 18f) 


(5.18g) 
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9.3.3 Development of the Multivariable Frequency Domain Design Method 
5.3.3.1 Direct Application of Nyquist Techniques to Each Input-Output 

The design problem reduces to the specification of the 
individual feedback gain kK. associated with each input-output pair of 
the augmented plant. Equation (5.18) can be used to analyse the system 
Stability directly, but even with the use of computers, this approach is 
impractical for destgn purposes and hence resort is made to graphical 
techniques analogous to conventional Nyquist analysis. 

Assuming the availability of an interactive computer program 
based on Equation (5.18) and a CRT terminal with graphic capability it is 
a relatively simple matter to display the Nyquist loci of the transfer 
function th, (z): i=] ... m}, relating each input-output pair. The gains 
ks can then be adjusted such that each of the loci satisfies the Nyquist 
Stability criteria and has the characteristics that imply satisfactory 
closed-loop performance. This approach represents an intuitive 
extension of frequency domain analysis techniques from single-input 
single-output (SISO) to multiple-input multiple-output (MIMO) systems. 
However it is only valid when the return-difference matrix, F(z), is 
diagonally dominant for all the z on the contour D (cf. Appendix D). 
NOTE THAT NO MODIFICATION OF THE "SISO" NYQUIST STABILITY ANALYSIS IS 
SUGGESTED, OR REQUIRED. THE ONLY DIFFERENCE BETWEEN "STSO" SYSTEMS AND 
"MIMO" SYSTEMS FOR MOST PRCTICAL APPLICATIONS IS IN THE DERIVATION OF 
THE TRANSFER FUNCTION RELATING EACH INPUT-OUTPUT PAIR. 

The stability analysis procedure outlined above is NOT 
practical for design purposes because, for example, a careful analysis 
of Equation (5.18) will shown that in addition to the requirement 


of diagonal dominance it is necessary to have specific numeric 
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values for m-1 of tne control parameters ks in order to evaluate any one 
transfer function, h(z). In general, Equations (5.16), (5.17) or (5.18) 
are too complicated to give much guidance for the selection of a 

feasible set of control parameters. Similarly even if the designer comes 
up with an initial estimate for each controller gain, k;;5 Equation (Sac) 
and its associated Nyquist diagram, provide little practical guidance 
about how the controller gains should be modified to improve the overall 
system performance. The next section will show how a practical, conven- 
jient design method can be developed from the application of Nyquist 


analysis to modified forms of Equation (5.18). 


5.3.3.2 Development of a Practical Frequency Domain Design Method 

The form of Equation (5.18) is relatively simple but the 
coefficients 65, (2) are in general functions of the controller gains in 
aldecnes ocvhere | oopsa andvhences itmaisedifticultytesdesignsthe first.m=) 
loops. The following subsections discuss various approximations for the 
interaction terms, L.(z), in Equation (5.18) which can be used to make 


the application of Nyquist techniques more practical. 


Cases with Zero Interaction (L. z) = 0 


A common, and practical approach for systems with very weak 
interactions, is simply to complete the control system design assuming 
that h. (z) = q.,(Z), which implies that the interactions are neglected. 
The same effect could obviously be obtained by setting L;(z) ='0 in 
Equation (5.18). 

An alternatiye approach is to specify K(z) such that the OLTFM, 


Q(z) is diagonal. This approach is normally referred to as “the non- 
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interacting design approach" and in terms of Equation (5.18) would result 


in values of zero for all the interaction terms L.(z). 


Using Limiting Values for L,(z) 


Instead of simply neglecting the interactions or specifying 
K(z) so that h,(z) = q,;(z) it is possible to find limiting values of 
L(z) which can be used in Equation (5.18). For example if it were 
possible to define an upper limit, Lamsee) on the magnitude of the 
interaction term L;(z), then jt follows that the magnitude of the transfer 


function relating the ‘hah input-output pair would be bounded by: 


fla;q(2) - U?%C2)p = Img (2d] < Alagg(2d] + LX(2)b (5.19) 


i 


In terms of a Nyquist diagram, Equation (5.19) shows that the Nyquist 
locus of the transfer function h, (z) lies within a band centered on the 


locus of q..(2z), and including the region defined by adding and subtract- 


Ws Gale 


ing, at each frequency, the value L. More specifically the band is 


defined by circles, with radii equal to be 


aA) centered on the loci 

of q..(z), for all frequencies. An example is shown in Figure 5.5. 
Although this approach may sound awkward it leads to a very practical 
design procedure which is in fact equivalent to the Direct Nyquist Array 
(DNA) Method discussed earlier in this chapter. The following sections 
deal with concepts such as Gershgorin circles that are familiar from the 


proceeding development of the DNA method. This is followed by extensions 


based on Equation (5.18). 
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Gershgorin Circles: It can be shown (see Appendix C) that if the return 
difference matrix, F(z), is diagonally dominant, (as defined below) then 


o.:(z) $ 1 for all values of controller gains and for all frequencies. 


iJ 
Thus it follows from Equation (5.18) and (5.19) that the Nyquist locus of 
nh. (z) lies within a band centered on the locus of q,;(z) and with a 


width defined by circles with the following radi: 


radius of Gershgorin circles = bea (2) =r F1q\2 


This band is in fact the "Gershgorin band" defined in Equations (5.8) and 
(OoeR (q; 


equations) It has the obvious advantage of being independent of the 


(z) does not contain ke and is defined the same in all 


controller gains and hence once the OLTFM is defined the Nyquist plots 
with the Gershgorin bands can be plotted for all diagonal elements, q,;(z). 
Since the locus of the exact transfer function, h.(z), must lie within 
the Gershgorin band the stability analysis can proceed as in the conven- 
tional Nyquist procedure but using the band in place of the locus. 
However if the critical point falls within the Gershgorin band then no 
decision can be made concerning the closed loop stability of that input- 
output pair. 

This approach clearly shows that the Gerhsgorin band represents 
a region that will contain the locus of h. (z) regardless of what gains 


k.. are used in the other loops - a fact that will be important to the 


ii 
understanding of system integrity. Also it is clear that the width of 
the band is a direct indication of the potential interactions of loop j 


with all the other loops in the multivariable system. 
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Diagonal Dominance: By analogy to the definitions used in the INA 
method, the return-difference matrix F(z) is said to be (columnwise) 


diagonally dominant if: 


.(z)| a apa (5.21) 


where f; 5 (z) represents element (i,j) of F(z). Since f,5(2) = q5; (2) for 
non-diagonal elements it follows directly from the physical significance 
of Equations (5.20) and (5.18b) that diagonal dominance is a guarantee 
that the critical point will not lie within the Gershgorin band and hence 
that this design approach is feasible. In practice the most convenient 
way of checking a system for diagonal dominance is to display the 
Gershgorin band associated with the Nyquist loci of the diagonal elements 
of Q(z) and see for what range of gain they satisfy the Nyquist stability 


criteria with respect to the critical point CUE OE 


Ostrowski Circles: As shown earlier in this chapter derivation of the 
direct Nyquist array method using theorems of linear algebra leads to the 
definition of the Ostrowski band. In terms of Equation (5.18) and (5.19) 


the Ostrowski band implies that when F(z) is diagonally dominant: 


k m 
max max 
LA) = reer Wo Gall 2h ys m (5.22) 
1 Dame iat Kop Z j= ji | 
we Die 


The Ostrowski band defined by Equations (5.19) and (5.22) can 
be deduced from Equations (5.16) through (5.18). It is narrower than 


the Gershgorin band but requires values of the gains, Kiss in all the 


108 


3 e «5 iT. 
i ; \ yak q 
~ 
4 
= 
xt 
. 
a >. - 
fog l Sie 
{ . 
o 7 “ = 
> 8a meth gig, 
, * 
c¢) enotteupay 
‘i e 4 
Se eae = 
4 i Vy ott ar 
: ae, ' = 
; } ; oh 
‘d ~ tj2 t+ 
{ pe 2Ths Dus 
- 


badnteoess based abyopd 


‘ 
> > 7 


" va" 
- 


—" 4 7 _ sak 
J Mae ‘ ssa ond 70 | nats g 


109 


other loops. As discussed later it is recommended that the "exact" 
relationship defined by Equation (5.18) be used rather than the 
Ostrowski approximation and hence Equation (5.22) is not discussed 


further. 


9.3.4 Extensions to the Direct Nyquist Array Method 

The DNA equivalent of the INA method developed by Rosenbrock 
involves a Nyquist analysis based on the diagonal elements of the 
augmented plant, q,,(Z), plus the bands defined by the Gershgorin and 
Ostrowski circles. However the intuitive development that led to 
Equation (5.18) gives greater insight into the physical significance of 
the DNA method and permits several extensions which make the method more 
practical as well as more accurate. 


For example once (m-1) of the gain elements kK. have been 
th 


determined, the m~ transfer function h,(z) can be evaluated directly 
from Equation (5.18) rather than using Gershgorin, Ostrowski, or other 
approximations. Also Equation (5.18) makes it easier to recognize cases 
where the interaction terms, L.(z), are zero and to calculate limiting 


m 


values Lez) that are more accurate than the Gershgorin circles. Some 


of these are discussed in the following sections. 


Cases witn Zero Interaction 
If the augmented system, Q(z), is dtagonaZ then, as pointed out 
previously, the radius of the Gershgorin circles is zero and the Nyquist 
analysis can be done on the basis of q;; (2). 
A less demanding approach is to design the matrix K(z) so that 


the augmented system OLTFM is triangular in form. This would mean that 
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all the interactions, L.(z), as defined in Equation (5.18) would be zero 
and hence the Nyquist analysis could again proceed on the basis of q.(z). 

The importance, to the design procedure, of zero elements in 
Q(z) does not seem to have been fully realized or utilized in previous 
proposals or implementations. Since Q(z) = G(z) K(z), it is possible 
in many practical applications to generate zero elements in specific 
elements of Q(z) by proper specification of K(z) ( Note that K(z) must 
be physically realizable and its poles and the zeros of det K(z) must be 
located inside the unit circle in the complex plane). If some of the 
interaction terms, L.(z), in Equation (5.18) can be made zero then this 
reduces the number of loops that must be designed by approximation and/or 
iteration. If enough zeros are introduced to make Q(z) triangular then 
the design of the MIMO system reduces to that of m SISO problems. 

The Gershgorin circles, as defined previously, can be an overly 
conservative and misleading bound on system interactions because they do 
not take into account the existence of zero elements in Q(z). For example, 
most triangular OLTFM would result in Gershgorin circles of nonzero 
radius which would imply uncertainty, or interactions, when in fact none 


exists. Therefore it is desirable to develop better limiting values. 


Limiting Values for the Interactions, L; (2) 


The Gershgorin circles result from the fact that the magnitude 
of the coefficients, 654 (2)> as defined by Equation (5.18d) is always 
less than one when F(z) is diagonally dominant. However when Q(z) 


contains zeros then some of the terms that contribute to the interactions, 


L.(z), will be zero regardless of the gains, kiss in the other loops. 


max 


In many cases a more accurate value of L. (z) can be derived if these 
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zeros are taken into account. 


G° Circles 
From Equation (5.18d) it is obvious that 65, (2) = 0 whenever 
the cofactor of element (j,i) of F(z) (or equivalently of F(z) or Q(z)) 
contains the equivalent of a row or column of zeros. This condition is 
obvious to the designer from the structure of Q(z) and can be defined 


formally as: 


m 
radius G° circles = Lam (22) = L ee feet Za) (923) 
# 


-6- 
Oo 
il 


0 when the cofactor of 95; (2) is zero or 1q55(2)| = 0 


1 when 65; (2) # 0 

In the general case when all the 054 (2) are non-zero the G° and Gersh- 
gorin circles are equivalent. On the other hand the G° circles are more 
accurate when some of the coefficients are zero and can be used whenever © 


F(z) is diagonally dominant. 


D7 Ss OmeCONGhUS TONS 

The most significant result of the preceding development of the 
direct Nyquist array method is the insight it gives to the physical 
problem. It is also important, from the point of view of the examples 
discussed in the following sections, that Equation (5.18) could be used 
as an "exact" Nyquist plot of hn. (z) as opposed to using the Ostrowski 


circles. 
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9-4 DESIGN OF A MULTIVARIABLE REGULATOR FOR THE THIRD ORDER EVAPORATOR 
MODEL 

This section presents the design of a multivariable regulator 
for the double-effect evaporator pilot plant described in Appendix A 
using the direct Nyquist array method. The model of the system is the 
Samenas, thesone Useds im SectionsS=3:o-ande423)and'it»consists of’ a discrete, 
three state, linear, time invariant model with normalized perturbation 
variables. The model was obtained from the tenth order model, derived 
by Wilson [77], using the Marshall's reduction technique [78] and is 
presented in Table A.3. 


The z-transform of the model is represented by the following 


equation: 
Wy (Z) S(z) 
Wo (Z) = | G(z) B, (z) (5.24) 
C,(z) Bo (z) 


where G(z) has been obtained from the state space model shown in Table 
A.3 by using an algorithm of Souriau-Frame-Faddeev [66]. G(z) is 
presented in Table 3.1. 

The outputs of the system are the first effect holdup, Wy 
the second effect holdup. Wo » and the product concentration, C,. 
Similarly the inputs are the steam flow to the first effect, S, the 
bottoms flow from the first effect, B,, and the bottoms flow from the 
second effect, B,. 

For reasons explained before, the design has been carried on 
using w-transfer function matrices instead of z-transfer function 


matrices so the procedure will appear to the user as identical to that 
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used for continuous systems. (w is a bilinear transformation of z such 


that z = {*W) 


/(1-w).) The plant w-transfer function matrix is presented 
in Table 3.2. The poles of the third order evaporator model (i.e. the 
zeros of the open-loop characteristic polynomial) are shown in Table 3.3. 
The system has two poles at the origin in the w-plane. However, if the 
Nyquist contour used to determine the stability of the system is the one 
shown in Figure 3.2 then the number of right-half-plane zeros (p,) of the 
open-loop characteristic polynomial is zero. 

The direct Nyquist array is more convenient to the user than 
its counterpart the inverse Nyquist array because of the information it 


provides. For example, the Nyquist array of G(w), shown in Figure 5.4, 


for the range of frequencies, 0.15 < w < 10.0, indicates that: 


1) The first effect holdup, Wy, is affected mainly by the bottoms 
Flow Bi; followed by the steam flow to the first effect, S. 

2) The second effect holdup, Wo » is affected mainly by the bottoms 
Flow By. However, the steam flow S and the bottoms flow Bo also 
affect Wo Significantly. 

3) The product concentration, Cos is mainly by the steam flow, S, 


followed by the bottoms flow B,. 
4) The bottoms flow B, does not affect C, and W. 


This information is not easily available from the inverse 
Nyquist array shown in Figure 3.3. The system in its present form jis 
highly interacting. This can also be observed from Figure 5.5 which shows 
the diagonal elements of G(w) with their respective Gershgorin bands. 
In each case the Gershgorin bands include the origin. (The labels in 


each diagram indicate the scale at the end of each axis.) 
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Step one: Pairing the output-input variables of the evaporator 


Interaction in the evaporator can be reduced by proper pairing 
of the variables. From the four points noted above, it can be concluded 
that the least interacting scheme is obtained when Wy is controlled by 
Bis Wo is controlled by B., and Cy is controlled by S. This means that 
column 1 must be made column 3, column 2 must be made column 1 and 
column 3 must be made column 2. This operation can be accomplished 


mathematically by post-multiplying the plant transfer function matrix by 


0 0 ] 
Ky ae i 0 0 (a) 
0 ] 0 


It should also be noted in passing that this first step of the 
Nyquist array approach can be regarded as an alternative to the method 
recommended by Bristol [8,9] to select the different control Joops in a 
multiloop scheme. 

The Nyquist array of the new open loop transfer function matrix 
Q, (w) = G(w) Ky is presented in Figure 5.6. The diagonal elements of 
Q, (w) with their respective Gershgorin bands are shown in Figure 5.7. 
Only the second column of the system Q, (w) is diagonally dominant. 


Therefore, it is necessary to proceed with another step in the design. 


Step two: Design of a pre-compensator to make the return-difference 


matrix, F(w), diagonally dominant 
The objective in this step is to design the simplest controller 


which will make the return-difference matrix, F(z), diagonally dominant. The 
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use of the Nyquist array plots is very helpful in this step. The Nyquist 


array of Q, (w) in Figure 5.6 suggests that: 


-column 1 can be made diagonally dominant by adding to this column 
2.1 times column 2. 
-column 3 can be made diagonally dominant by subtracting from this 


column 0.93 times column 2. 


These two elementary column operations can be done by post- 


multiplying Q, tw) by the following controller: 


120 0.0 O20 
Ke See lec cl ied) De93 (5325) 
0.0 0.0 ive] 


The new open-loop transfer function matrix is then equal to 
Qo(w) = Q)(w) Ky = Q(w) KyK, (Scie) 


Every column in Qo (w) is now diagonally dominant as can be 
seen from the Nyquist array of Qo (w) in Figure 5.8 and from the diagonal 
elements of Qo (w) with their Gershgorin bands in Figure 5.9. However, 
the diagonal dominance of column 3 can be improved further by subtracting 


0.4 times column 1 from this column 3, i.e. by post-multiplying matrix 


Qo (w) by: 


[oom 90.4 | 
Selo. ye 0 (5.28) 
Oe ef i0 


Figure 5.10 shows the Nyquist array of the new open Joop 
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transfer function matrix Q, (w) = Q, (w) Ke and Figure 5.11 shows its 
diagonal elements, with their respective Gershgorin bands. It can be 
seen from these figures that columns two and three are almost perfectly 
diagonally dominant (i.e. imply "non interaction") and that diagonal 
dominance can be improved in column 1 by adding 0.63 times column 3 to 


column 1, i.e. by using the controller: 


] 0 0 
Ky = fal ] 0 (5229) 
Oe OS meu ] 


The Nyquist diagrams of the diagonal elements of Q, (w) = Q3(w) Ky 
in Figure 5.12 indicate that Q, (w) is strongly diagonally dominant, i.e. 
the Gershgorin bands can hardly be noticed in this figure. 

Having obtained a controller matrix which makes the OLTFM 
diagonally dominant the stability of the closed-loop system can now be 
determined. Using tne Nyquist contour shown in Figure 3.2, Die 0 and 
hence the closed-loop system will be stable if none of the loci of the 
diagonal elements encircles the critical point (-1,0). In order to 
satisfy the stability criterion the sign in columns 1] and 2 must be 


changed. Thus it will be necessary to post-multiply Qy (w) by: 


= 0 0 
K, = 0 =| 0) (5230) 
0 0 ] 


The diagonal elements of Q. (w) = Q, w) Ke are presented in 
Figure 5.13 and it is obyious that none of them encircle the critical 


point. Therefore the system is stable. Thus the controller that will 
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Lay, 


make the OLTFM diagonally dominant and simultaneously make the closed- 


loop system stable is given by: 


Kp = Ky Ko Kg kg Ke 
Cie 
=0-030 0.0 0 
Kp AD ead Mats) 0.0 -0.400 


Als Seis Sig) OPRAH 


Step taree: Design of individual controllers for each loop 


(5S 1) 


(E32) 


Each loop in the system (i.e. the diagonal elements of the OLTFM) 


can be adequately represented by a first order transfer function. Thus the 


use of simple proportional, or proportional plus integral control in each 


loop should be sufficient to provide very good control. 


The gain margins of each loop can be estimated graphically from 


Figure 5.13 using the same approach as in the design of single variable 


control systems. The system will be closed-loop stable if 
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oop | O< kK, < 0.0406 24.6 
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oop 2 (Of) & Ky < 970203 AQ? 
i Todos. 
oa 0 < k, < 770360 = 27-8 


The final controller matrix will be given by: 


IS Kako 


where for proportional control Ke will be equal to 


Ko = diag ( k, ok» sk3) 


(5333) 


(5.34) 


(525) 


(5.36) 
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and for a proportional-integral control K. will be given by Equation (3.35). 
It should be noticed that the controller Kp and the stability 

margins given by Equations (5.32),(5.33),(5.34) and (5.35) are almost 

identical to the controller obtained using the characteristic locus method 

in Section 4.3. Also for the third order model it is possible to calculate 

a controller which will give a "perfectly" non-interacting system. This 

was done in Appendix E and the controller was almost identical to Equation 

(5.32). Thus in this example the "graphical" frequency domain design 

procedures led to essentially the same result as the analytical calcu- 


lation. Another interesting result is obtained by letting k, = 1.337, 


if 
ky = 1.0 and kz = 0.751 in Equation (5.37); the new controller matrix is 


then equal to: 


-0.842 0.0 Chad a 
Ke Keke eee COU 0.0 =U5 300 (5733) 
eo 12 Sia) aah eres 


and the new stability region will be given by: 
Loop 1 0 < ky < 357 = 18.4 (5.39) 
Loop 2 | Ok ko Soap = 49.2 (5.40) 
Loop 3 Bs Ten 0.751 > Sid (5.41) 


This is the same basic controller matrix with the same stability 
region obtained for the third order evaporator model using the inverse 
Nyquist array and the characteristic locus method in Sections 3.6 and 


4.3 respectively. Thus the simulated runs shown in Figures 3.11, 3.12 
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and 3.13 in Section 3.6 also applies for the controller matrix given by 
Equation (5.32). The application of the Ziegler-Nichols settings rule 

to select the final gains produced excellent control in the simulated 
runs but they were found to be too high when these gains were used in the 
pilot plant. For this reason a decision was made to tune the loops 
experimental ly on the proportional and proportional plus integral case. 
The selection of the matrix Ko by experimental tuning is presented in 


Chapter Six. 


9.5 DESIGN OF A MULTIVARIABLE REGULATOR FOR THE FIFTH ORDER EVAPORATOR 
MODEL 
In this example a multivariable regulator is designed using the 
direct Nyquist array for the double-effect evaporator plant described in 
Appendix A. The model used in this example is the same as the one used 
iMesect)0n 3./5 i.e. a discrete, fiye-State, linear, time invariant model 


as presented in Table A.2. The evaporator is represented in the z-domain 


by: 
Wy (z) 32) 
Wo (Z) = | G(z) B, (z) (5.42) 
C,(z) B, (z) 


where the plant z-transfer function G(z) is shown in Table 3.5. The 
design of the multivariable regulator, as in the previous cases and for 
reasons explained before has been done using a w-bilinear transformation 
of the plant z-transfer function matrix. Thus G(w) has been used instead 


of G(z), where z = (ltw)/(1-w). G(w) is shown in Table 3.6. 
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The poles of the plant transfer function matrix or the zeros 
of the open-loop characteristic polynomial are shown in Table 3.7. The 
plant has two poles at the origin in the w-plane. The Nyquist contour 
used to determine the stability of the system is shown in Figure 3.2. 
Using this contour the number of right-half-plane zeros (p,) of the 
open-loop characteristic polynomial is equal to zero. The Nyquist array 
of G(w) for a range of frequencies 0.1 < w < 5.0 is shown in Figure 5.14. 
(Different ranges of frequencies have been used in the different 
examples in this thesis in order to have an appropriate scale in each 


diagram.) It can be observed from Figure 5.14 that: 


1) the first effect holdup, Wy» is mainly affected by the bottoms 
flow B,. followed by the steam flow, S. 

2) the second effect holdup, Wo. is mainly influenced by the bottoms 
Flow B, and followed, at high frequencies, by the bottoms flow 
Bo and finally by the steam flow S. Note that at low 
frequencies the steam flow S has a greater effect than B. in 
controlling Wo. 

3) the product concentration Cy is affected almost equally by the 
steam flow S and the bottoms flow B,. 

4) the bottoms flow B, does not affect C, or W. 

This is essentially the same information obtained in the 
previous example which gives an indication of the accuracy of the third 
order model. This information also indicates that the present configura- 
tion of the double-effect evaporator is highly interacting. This can 
also be observed in Figure 5.15 which shows the Nyquist diagrams of the 


diagonal elements of G(w) with their corresponding Gershgorin bands. 
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Interaction can be reduced in this case, without the use of a compensator 


by properly pairing the input-output yariables. 


Step one: Pairing the evaporator input-output variables 


From the information obtained from the Nyquist array of G(w) 
in Figure 5.14 it can be concluded that the least interacting configura- 
tion is obtained when B, controls Wis Bo controls Wy and S controls C,. 

This first step in the design involves a renumbering in the 


columns of G(w). This is obtained by post-multiplying the OLTFM by: 


0 0 ] 
Ky =F ih 0 0 (5.43) 
0 1 0 


In Figures 5.16 and 5.17 the Nyquist array of Q, (w) = G(w) K 
and its diagonal elements with their Gershgorin bands are presented, 
respectively. Only column 2 is diagonally dominant. Consequently it is 
necessary to design a pre-compensator to make the system diagonally 


dominant. 


matrix F (w) diagonally dominant 
Step 2a) Design of constant pre-compensator 
This step involves the selection of the simplest controller 
that will make the return-difference matrix diagonally dominant, i.e. 
the selected controller must produce an open-loop system in which the 
Gershgorin circles of its diagonal elements do not enclose the critical 
point (-1,0). The use of the Nyquist array is again very helpful in this 


selection. 
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From Figure 5.16 it can be seen that column 1 can be made 
diagonally dominant if 2.0 times column 2 is added to column 1, i.e. if 


the following control matrix is used: 


le) 0.0 0.0 


K, = | 2.0 1.0 0.0 (5.44) 
0.0 0.0 1.0 


Gershgorin bands are shown in Figures 5.18 and 5.19 respectively. Columns 
1 and 2 of Q, (w) are now diagonally dominant. However the diagonal 
elements are unstable. These loops can be made stable by changing the 


signs in these two columns, i.e. by using the following controller matrix: 


= 0 0 
k=] 0 = 0 (5.45) 
0 0 1 


The new OLTFM is equal to 
Q3 (w) e Q,(w) K3 (5.46) 


and its Nyquist array and diagonal elements are shown in Figures 5.20 
and 5.21. Column 3 is still not diagonally dominant. However, it is 
obvious from Figure 5.20 that interaction on be reduced in this column 
by adding 0.4 times column 1 to column 3, i.e. by postmultiplying Q,(w) 


by the controller matrix: 


K, = | 0 1 0 (5.47) 
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Figures 5.22 and 5.23 present the Nyquist array of Qy(w) = 
Q3 (w) Ka and its diagonal elements with their respective Gershgorin 
bands. Diagonal dominance in column 1 can Slightly be improved by 
subtracting 0.6 times column 3 from column 1. This operation can be 


accomplished with the controller matrix: 


leo 0.0 0.0 
Ke = | 0.0 eo 0.0 (5.48) 
=O 36 0.0 ee0 


The Nyquist diagrams of the elements of Q. (w) = Q, (w) Ke are 
shown in Figure 5.24 and the Nyquist diagrams of the diagonal elements 
with their Gershgorin bands are shown in Figure 5.25. The open-loop 


transfer function matrix Q, (w) js not diagonally dominant in column 3. 
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A careful examination of Figure 5.24 suggests that the diagonal dominance can 


be improved in columns 1 and 3 but it will not be possible to make 
column 3 of the OLTFM, Q. (w) diagonally dominant using a constant 
pre-compensator. However it can be noticed in Figure 5.25 that the 


return-difference matrix, i.e. 
F(w) = I+ Q.(w) (5.49) 


is diagonally donenene (The Gershgorin circles do not enclose the 
Critical point (1,0); Thus it is possible, using a constant compensator, 
to continue in this case with the design. hts was not the case when the 
inverse Nyquist array method was used to destgn a static regulator for 
thts system. 

The pre-compensator used so far to make the return-difference 


matrix diagonally dominant is giyen by; 
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or 


-0.600 0.0 1.0 
Kn = 1704/60 0.0 -0.400 (51) 
1 eal) =i hot -0.800 


From the Nyquist array of Q., (w) in Figure 5.24 it can be 
noticed that diagonal dominance of F(w) could be improved a little in 
column 1 by subtracting 0.5 times column 2 from column 1. Similarly 
diagonal dominance in column 3 could also be improved by adding column 
2 to column 3. This operation can be accomplished by post-multiplying 


Q. (w) by matrix Kes where 


1 0 0 
Ke = [-0.5 1 1 (5.52) 
0 0 1 


The Nyquist array of Q6 (w) = Q (Ww) Ke» and the Nyquist 
diagrams of its diagonal element witn their respective Gershgorin bands 
are shown in Figures 5.26 and 5.27, respectively. The new pre-compensator 
matrix, Kp» used to make the return-difference matrix diagonally dominant 


is then given by: 


or 
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-0.600 0.0 Vay 
Kp I AO TAS 0.0 -0.400 (Sasa) 
2020 =| .0 -1.800 


The compensator Kp given by Equation (5.54) is approximately the 
Same as the one obtained in Section 5.3 using the direct Nyquist array 
method and the third order model. (This is another indication of the 
accuracy of the third order model.) 

The degree of diagonal dominance obtained in the return- 
difference matrix F(w) using the pre-compensator Ke is about the best that 
can be obtained when a constant pre-compensator is used. If better 
diagonal dominance is desired in the return-difference and in’ the open- 
loop transfer function matrix a dynamic compensator has to be designed. 

In order to investigate the effect of diagonal dominance in the perform- 
ance of the control system a dynamic pre-compensator was also designed 
as part of this example to improve the degree of diagonal dominance in 


the open-loop transfer function matrix. 


Step 2b: Design of a dynamic pre-compensator 
It was obseryed from the Nyquist array of the OLTFM Qo (w) and 


from the Nyquist diagrams of its diagonal elements in Figures 5.18 and 
5.19 respectively that columns |] and 2 of Q, (w) were diagonally dominant. 
A careful examination of the poles and zeros of the plant transfer 
function matrix and the Nyquist array of Qo (w) in Figure 5.18 suggested 
that column 3 could be made diagonally dominant if a phase-lead 


compensator is used and the following elementary operation is performed. 


column 3 = column 3 - 0.0559 (ee ) x column 2 (5.54) 
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This operation can be realized if the matrix Q (w) is post-multiplied by 


the following controller matrix: 


1 0 0 
] a ” w+ 7.400 
K3 (w) 0 ] 0.0559 G eraRaGa ) (5.55) 
0 0 ] 


The new open-loop transfer function matrix, Q3 (w) will be given 


by: 


Q3 (w) = Q,(w) Ky (w) (5.56) 


The Nyquist array of Q5 (w) and the Nyquist diagrams of its 
diagonal elements with their Gershgorin bands are presented in Figures 
5.28 and 5.29 respectively. It can be observed from Figure 5.29 that the 
three columns are diagonally dominant. However columns 1 and 3 are 
weakly diagonally dominant. The Nyquist array of Q5(w) in Figure 5.28 
and an analysis of the zeros and poles of this matrix suggests that column 
3 can be made almost diagonal if a phase-lead Rernensater is used and the 


following operation is performed: 


column 3 = column 3 - 0.024 (2 x column J (5.57) 


This operation can be accomplished if the matrix Q5 (w) is post- 


multiplied by 
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w + 7.400 
K'(w) = | 0 0 
mi (5.58) 
0 0 ] 


The Nyquist array of the new OLTFM 04 (W) = 05 (w) Ka (Ww) is 
presented in Figure 5.30. The Nyquist diagrams of the diagonal elements 
of a; (w) with their respective Gershgorin bands are shown in Figure 5.31. 
From this figure it can be observed that only column 1 of Q; (wr) does not 
have very good diagonal dominance. However an analysis of the Nyquist 
array of Qj (w) in Figure 5.30 and of the poles and zeros of this matrix 
Suggests that diagonal dominance in column 1 can be drastically improved 
if a phase-lag compensator is used and the following operation jis 


performed on aq (wv): 


column 1 = column 1 + 13.2 (ho x column 3 (5.59) 


This operation can be realized if the matrix Q, (w) is post- 


multiplied by: 


1 0 0 
Kr (w) =| ie 1 0 (5.60) 
w + 0.390 
“13.20, 4-7-4997) «= | 


The new open-loop transfer function matrix, a5 (w) is given by 


Qn(w) = Qg(w) Ke (w) (5.61) 


The Nyquist array of Q5 (Ww) and Nyquist diagrams of their 
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diagonal elements with their Gershgorin bands are presented in Figures 
5.02 and 5.33) respectively. From Figure 5.33 it can be observed) that 
Qs (w) has a high degree of diagonal dominance (it is.almost diagonal). 
However two of its diagonal elements are unstable. To correct this 


Situation it is necessary to post-multiply 05 (w) by 


ey Canben 0 
Keacalig0for Sle e000 (5.62) 
1 ae 


The Nyquist diagrams of the diagonal elements of the new open- 
loop transfer function matrix 5 (w) = Q5 (w) ke are presented in Figure 
DEo4. 

The dynamic compensator designed in this section to make the 


open-loop transfer function matrix diagonally dominant will be given by 


Ka(w) = KjKoKg(w) Ky(w) Ke(w) Kp (5.63) 
or 
“13.2 (WES) 0 1 
Kn (w) = -0.683 0 © -0.024 (AEF AO) | (5.64) 
-0.629 1 -0.1039 (MEF 


The z-transform of this controller is equal to 
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Zz - 0.4383 
“ele Weererye rage a ei 
res ; ; z + 0.7619 
Kh (z) 0.683 0 0.145 (2-7-5) | (5.65) 
A é s Zee LO Lg 
0.629 1.0 -0.628 (2—~7) 


Step three: Design of individual controllers for each loop 
Step 3a. Stability analysis for the constant compensator 

The controller matrix K, given by Equation (5.53) was designed 
to produce diagonal dominance in the return-difference matrix and is 
typical of the best results that can be obtained using a compensator of 
constants. The design can be continued by selecting the individual 
controllers in each loop such that the final controller matrix is given 


by; 
K = KRKe (5.66) 


where Ke is diagonal. It has the following form when only proportional 


control is desired: 


Ky 0 0 
Ke = 0 ko 0 (5:67) 
0 ky 


The main objective in this step is the selection of the constant 
ky ; kK, 3 k3 . In order to select these constants it is yery useful to 
first determine the stability region for which the closed-loop system will 
remain stable, i.e. the gain margin of each loop. A conservative 


stability region can be obtained from the Nyquist diagrams of the diagonal 
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elements of the OLTFM Q,(w) containing the respective Gershgorin bands in 
Figure 5.27. The closed-loop system in this case will remain stable for 


the following combination of gains. 


e - 1 = 

Loop le 0 < ky Ss 0.0693 aa Tana (5.68) 
: | i 

Loop 2: O0< kK, <TR 49.2 (5.69) 
e 1 = 

Loop 3: 0 < ky < 0.116 ae 8.62 (5570) 


The stability region given by Equations (5.68), (5.69) and (5.70) 
can be graphically represented by Figure 5.35. The closed-loop system 
will remain stable for any combination of gains inside the parallelopiped. 
Each side of this stability region can be expanded by making use of the 
Nyquist diagram of the exact transfer function of each loop when the 
rest of the loops are closed as given by Equation (5.18). However 
the determination of the expanded stability region involves a very 
complicated procedure and its complexity increases with the order of 
the system. The new region which is irregular has to be determined 
pointwise and each point has to be found by a trial and error procedure. 
In addition the use of Equation (5.18) will produce a discontinuous 
Stability region since it is only possible to find the expanded stable 
Space for the region which opposes each side of the paral lelopiped,: as 
shown in broken lines in Figure 5.35. This is because the use of the 
transfer function {h;(z),i=1,..m} given by Equation (5.18) only 
guarantees stability when the return-difference matrix is diagonally 
dominant,i.e. only one value of k. can be varied in each trial and error 


attempt and the values of the other gains must be within the region 
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defined by Equation (5.68)-(5.70). Because the benefits of obtaining 
an expanded stability region are marginal the calculation of the 
expanded gain space is not recommended and it was not done in this 
example. 

Once the stability region has been obtained it is possible 
to select the final gains in each loop. In this example a direct 
application of the Ziegler-Nichols controller setting rules is not 
feasible since the gain margins obtained do not represent the ultimate 
gain of each loop. Therefore the selection of the final gains has to be 
done in a somewhat arbitrary way. One logical way to make this selection 
is by choosing the gains based on the conservative gain margins given 
by Equations (5.68), (5.69) and (5.70) and by using higher percentages 
of these values when the uncertainty is greater i.e. when the Gershgorin 
circles are larger. Additional tuning of the loops may be required for 
these cases. 

An arbitrary selection in this example of a gain which is 34% 
of the gain margin for the first loop, 20% for the second loop and 65% 
for the third loop will produce the controller FDO320 shown in Table 
B.1 in the Appendix B. Controllers FD0310 and FD0330 can be obtained 
by using different percentages of the gain margins. Figures 3.36 and 
3.38 present in broken lines the simulated response of the evaporator 
using controller FD0320 for a 20% step change in the feed flow and 
10% step change in the setpoint of the product concentration, respectively. 

The experimental tuning of the compensator given by Equation 


(5.54) is presented in Chapter Six. 
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Step 3b. Stability analysis for the dynamic compensator 
Note that the dynamic compensator given by Equation (5.65) is 


very similar to the compensator for "perfect" non-interacting control 

(cf. Appendix E). The stability region can be determined using the Nyquist 
diagrams of the diagonal elements of 5 (w) in Figure 5.34. From this 
figure it can be established that the closed-loop system will be stable 


for the following combinations of gains: 


Loop 1: 0 < kK < D030 = 24.1 (5.74) 
: | = 

Loop 2: 0<k, < 0.0217 = 20-6 (575) 
- ] » 

Loop 3: 0<k3< 00475 «= 31-8 (Say 6) 


The final control matrix in this case will be given by: 
K'(z) = K(z) kK (5.78) 
B i ; 


where Ke is-also given by Equation (5.67) when only proportional control 
is desired. 

It should be noted that if Ko in Equation (5.67) is selected, 
such that k, = 1.464, ky = 1.0 and k, = 0.668, the controller k'(z) 
given by Equation (5.78) and its stability region will be approximately 
equal to the dynamic compensator and the stability region obtained in 
Section 3.7 using the inverse Nyquist array and the fifth order model. 

The selection of the final gains in this case was also done by 
experimental tuning. However the simulated responses of the evaporator 
for a 20% step change in the feed flow, when the gains used in each loop 
are 10% (controller NADY0510), a 20% (controller NADY0520) and 30% (controller 


NADY0530) of its corresponding ultimate gain, are also represented by 
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the continuous line in Figures 3.35, 3.36 and 3.37 respectively. In 
this case there is also an improvement in the control system when the 
gains in each loop are increased. This improvement is mainly reflected 
in the first effect level. 

Figure 3.38 presents the simulated response (continuous lines) 
of the evaporator when controller NADY0520 is used for a 10% step change in 
the setpoint of the product concentration, C,. Since the system is 
non-interacting when the dynamic compensator is used only the product 


concentration is affected in this case. 


Effect of Improving the Diagonal Dominance in the Evaporator 

In the previous example it is possible to invesigate the effect 
of improying the diagonal dominance in the system. Three compensators 
were obtained that make the open-loop and return-difference matrices 
diagonally dominant to different degrees. These compensators were Kn 
given by Equation (5.51), Kp given by Equation (5.53) and Ka given by 
Equation (5.64). The degree of improvement in the diagonal dominance of ° 
the system can be observed in Figures 5.25, 5.27 and 5.34 by comparing 
the radii of the Gersngorin circles. The compensator Kp makes the 
system almost diagonal while Kn gives the poorest diagonal dominance. 

Figure 5.36 presents the response of the open-loop compensated 
plant for a step change in the input corresponding to C,. As expected, 
improving the diagonal dominance of a system makes it less interacting. 

Figure 5.37 and 5.38 present the simulated response of tne 


evaporator for a step change in feed flow and in the setpoint of the 
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product concentration using controllers DNA0520, FD0320 and NADY0520. These 


controllers are based on the compensators Kn Kp and K respectively 
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and have gains which are of the same magnitude. The controller NADY0520 
was used as a base case and its gains are 20 % of the corresponding 
ultimate gains. It can be observed from these figures that the performance 
of the evaporator improves slightly by improving the diagonal dominance 
in the OLTFM not only for a setpoint change but also for a step change 

in a load disturbance. It should also be noted that the improvement is 
smaller for a step change in the setpoint of C5 than for a step change 

in the feed flow. 

The effect of improving diagonal dominance in the system is 
also shown in Figures 3.35, 3.36, 3.37 and 3.38 which compare the 
performance of the evaporator using constant versus dynamic controllers. 
It can be noticed that when the gains are approximately the same for the 
dynamic and the constant compensators the response when using the dynamic 
compensator is slightly better than when the constant compensator is used. 
This is the case not only for a step change in a setpoint but also for a 


step change in a disturbance variable. 


9.6 CONCLUSIONS 

The direct Nyquist array method has been found to be a yery 
practical and powerful design technique which arises intuitively by 
extending the conventional single-input, single-output Nyquist design 
procedure. This method has the advantage of requiring no new theory and 
all the concepts and mathematical operations involved are familiar to most 
control engineers. The direct Nyquist array method provides a systematic 
means of implementing control system design objectives such as: stabi- 
lity, reduction of interactions, high system integrity against component 


failures, good steady-state accuracy and acceptable transient responses. 
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The direct Nyquist array method is more practical than the 
characteristic locus method [4,42], because the Nyquist array plots 
provide more insight into the performance of the system and an indication 
of how to design a controller to produce a particular closed-loop 
characteristic rather than simply providing a test of whether or not the 
the desired result has been obtained. It must be noted that a procedure 
for the characteristic locus method has been developed recently by 
Kouvaritakis [508] to design a control system in a more systematic way. 

The direct Nyquist array method is more practical than the 
inverse [65] method because the designer specifies the controller elements 
directly (rather than specifying elements of the inverse of the control 
matrix and having the computer calculate the actual values later). Thus 
the designer has a better "feel" of the complexity of the actual controller 
he is developing. Also the direct Nyquist array consists of "standard" 
Nyquist diagrams that represent components of the actual physical system 
(rather than their inverses) so most users would get a better physical 
insight than with the inverse method. With the direct method it is also 
possible to relax the requirement of diagonal dominance required by ite 
inverse method and hence is applicable to a broader class of applications. 
From the computational point of view the direct Nyquist array method is 
more efficient than the inverse method since one extra operation is 
required in the latter method, i.e. the inversion of the plant transfer 


function matrix. 
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CHAPTER SIX 
EXPERIMENTAL RESULTS AND COMPARISON 
WITH OTHER METHODS 


6.1 INTRODUCTION 

This chapter discusses the experimental evaluation of the 
multivariable regulators, designed in Chapters Three, Four and Five using 
the inverse Nyquist array [65], the characteristic locus [4,42] and the 
direct Nyquist array methods respectively. All of the controllers were 
designed for, and applied to the double effect evaporator described in 
Appendix A. These methods are also compared with an optimal-quadratic 
multivariable controller and a multiloop control scheme that haye been 
applied previously to the same evaporator [27,53,77]. 

Each control system was eyaluated by introducing a step change 
in one of the disturbance variables of the system. Three different 
disturbances can be used to upset the pilot plant evaporator: feed flow, 
feed concentration and feed temperature. Only the first two have been 
used in this study. Feed flow has been found to be the most seyere 
disturbance in the system and this was the main reason for using it. The 
feed concentration disturbance has been used to evaluate the degree of 
interaction in the control system. The degree of interaction of the 
closed-loop system was also tested by introducing a step change in the 
setpoint of the product concentration in the double effect evaporator. 

The experimental runs described in this chapter were done in 
two sets separated by several montns during which the evaporator was 
operated almost continuously and underwent a major overhaul. The second 


set of runs was done primarily to evaluate a dynamic compensator that was 
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designed in the last stages of this investigation but also included some 
other runs for comparison. It should be noted that in the second set of 
runs a slightly lower (i.e. 80%) effective gain was used in the C. 
measurement transducer and also that there is slightly less noise evident 
in the system. For these reasons comparisons are made only between exper- 
imental runs from the same set. However both sets of runs lead to the 
same general conclusions. . 

The controllers were implemented on an IBM 1800 digital data 
acquisition and control computer, which is interfaced to the pilot plant 
evaporator, using a computer control package developed by Newell [53]. 

There were a number of factors that had to be considered when 
evaluating the different design methods and the performance of the 


controllers produced by each method. These included: 


1) Design metinod 
- inverse Nyquist array 
- direct Nyquist array 
- characteristic locus metnod 
- optimal quadratic regulator 
- conventional single yariable, multiloop control 
2) Model used as the basis for the controller design 
- third order state space evaporator model 
- fifth order state space evaporator model] 
3) Control modes 
- proportional 
- proportional plus integral 


- open-loop 
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4) Control performance 
- under load changes 
- under setpoint changes 
5) Augmented plant 
- compensator of constant elements 
- dynamic compensator 
6) Method of evaluation 
- simulated runs 


- experimental runs 


It was not feasible to investigate all possible combinations of these 
different factors but tne controllers evaluated as part of this study and 
the figure numbers that present the corresponding dynamic response of the 
evaporator variables, are summarized in Table 6.1. 

After the experimental results have been presented and 
discussed, the conclusions regarding the most important of the above 
factors are summarized. Unfortunately there is no single, quantitative 
performance criterion that adequately applies under ali conditions. 
Therefore in most cases the conclusions are derived by careful, qualita- 
tive comparison of the runs as grouped in Table 6.1. 

The first sections of this chapter discuss the experimenta} 
tuning that was required to determine the final gains in the controllers 
designed using the inverse Nyquist array, the direct Nyquist array and 
the characteristic locus methods. This is followed by the experimental 
results using an optimal multivariable controller and a multiloop control 


scheme. 
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6.2 EXPERIMENTAL TUNING OF THE CONSTANT COMPENSATOR OBTAINED FOR THE 
THIRD AND FIFTH ORDER MODEL USING THE INVERSE, THE DIRECT NYQUIST 
ARRAY, AND THE CHARACTERISTIC LOCUS METHODS 

For the third order evaporator model it was noticed that the 

Static compensator designed with the inverse Nyquist array method [65] in 

Chapter Three and with the characteristic locus method [4,42] in Chapter 

Four were practically the same as the one obtained using the direct 

Nyquist array metnod in Chapter Five. Furthermore it was pointed out also 

in Chapter Five that tne constant compensator obtained for the fifth 

order evaporator model using the direct Nyquist array method was almost 
equal to the one obtained for the third order model. Since the same 
basic compensator was obtained in all these cases this section covers 


the tuning of the constant compensator for all these cases. 


6.2.1 Selection of the Proportional Gains 

The use of the Ziegler-Nichols setting rules [12] to select 
the final proportional gains in each loop produce a set of gains which 
were found to be relatively high when they were applied in the field. 
The response of the evaporator was not as good as the simulated runs had 
indicated. Therefore the proportional gains in each loop were obtained 
by experimental tuning. This was done in a systematic way. Starting 
with the basic information obtained in Chapters Three, Four and Five. 
The tuning of the constant compensator for the third order model proceeded 
as follows ( For the fifth order model a more conservative stability 
region was obtained Sf different percentages than the ones indicated 
below have to be used to produce the results presented in this section.) 


Given the basic controller: 
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and the gain space for which the system remains stable: 


Loop 1: 0 < kK, <p 2/086 (6.2) 
Loop 2: OL < keaelo ae (6.3) 
Loop 3: 0 < kz < 27.8 (6.4) 


The objective was to find the matrix Ko. 


Ky 0 0 
Ke = | 0 ko 0 (6.5) 
0 0 Ky 
such that the control system 
NS Kako (6.6) 


gave the best response for a step change in the setpoints of the controll- 
ed variables or for a step change in the disturbance variables. 

The procedure followed to choose Ky > ko and Ky was very 
similar to the Ziegler-Nichols approach [12] commonly used to select the 
controller settings in single-loop systems. The chosen gain in each 
loop (ky Ko» kz) was always based on its respective ultimate gain. The 
gains used in each loop to control the evaporator were 10% (controller 
FD0310), 20% (controller FD0320) and 30% (controller FD0330) percent of 
the corresponding ultimate gain of each loop. The numerical values of 


these controllers are presented in Table B.1. 
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Figures 6.1, 6.2 and 6.3 present the experimental response of 
the closed loop system using the controllers FD0310, FD0320, and FDO330 
respectively for a step change of 20% in the feed flow. As expected by 
increasing the gains in each loop the offsets are reduced but the system 
becomes more oscillatory. The effect of increasing the gains can be 
clearly noticed in the three output variables which are being controlled: 
first effect level, Ws second effect level, Wo» and the product 
concentration, C.. When using controller FD0330 in this set of runs 
there was a sustained cycling in the system (Figure 6.3), indicating that 
the evaporator is on the verge of instability. 

From these runs it was concluded that a reasonable set of 
gains are the ones used in controller FD0320, which produced the results 
plotted in Figure 6.2. (A later run with a gain equal to 16% of the 
ultimate value for the C,-S loop and done when there was less noise evident 
in the evaporator system, is plotted in Figure 6.15 and has an even 
better response. ) 

The fact that the basic multivariable controller can be tuned 
experimentally using procedures similar to the ones used for single-loop 
controllers is probably one of the most important advantages of the 
frequency domain techniques studied here. With these methods if the 
process response is Rest ecco there is usually no need to completely 
redesign the controller. A satisfactory control system can usually be 
obtained by increasing or decreasing the gains in each feedback loop. 
This is not the case with many of the multivariable time domain methods, 
since if the design criterion is to be maintained an unsatisfactory controller 


has to be redesigned, re-evaluated and finally re-implemented experimentally. 
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Another important advantage of the frequency-domain methods 


discussed in this study is that with these methods it is always clear to 


the control engineer how the controller should be modified if it is not 


satisfactory. With other methods it is not always clear what must be 


modified, or how to obtain a satisfactory control system. 


The proportional controller FDO0320 was also tested for a 20% 


step change in feed concentration with the results shown in Figure 6.4. 


In this run there is no significant change in any of the output variables. 


Because the feed concentration affects only the product concentration 


this run is a measure of the interaction in tne controller. The system 


is almost non-interactive as expected. (An unexplained disturbance in the 
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steam flow, at T=104 minutes affected all the output and control variables.) 


The degree of interaction of the control system using the 


proportional controller FD0320 can also be observed in Figures 6.5 and 


6.6. These figures show the response of the eyaporator to a 10% step 


change in the setpoint of the product concentration, C,. 


the system is fast and the interaction between loops is low. For a 10% 


change in the product concentration, Cos there is a disturbance of only 


4.5% and 6.0% in the levels of the first and second effect respectively 


and they recover after a relatively short period of time. The system is 


a little more oscillatory for a step-up than for a step-down jn the 


setpoint. 

From the previous experimental runs 
ultimate gains obtained experimentally do not 
ones obtained from the Nyquist diagram of the 


models. The experimental ultimate gains (see 


are approximately 30% of the theoretical ones. 


it can be noticed that the 
agree with the theoretical 
third and fifth order 

run FD0330 in Figure 6.3) 


The difference is 


The response of 
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attributed to modelling errors. However the simulated results presented 
in Chapter Three (Figures 3.11, 3.13, 3.35-3.38) are qualitatively 
consistent with the experimental results. In both cases the main effect 
of increasing the gains in the controller is to reduce the offset in the 
first-effect level, Ws The increase of gains in each loop does not 
affect the second-effect level, Wo» and the product concentration, Cys 
very much. 


bec cape lecuioniotathesintegraly Constants 


In the previous section the most reasonable proportional gains 
were found to be the ones used in the controller FD0320. Consequently 
these gains were used as a Starting point in the selection of the 
proportional plus integral controller constants. In the third order 
model (cf. Tables 3.1 and 3.2) and to a lesser extent in the fifth order model, 
every open-loop transfer function in the system behaved approximately 
like a first-order system. Thus it was not possible, in this case, to 
use the crossover frequencies to estimate the integral constants. 
Therefore they were arbitrarily chosen to be 1/64 of the gains used in 
controller FD0320. This integral controller FD1320 shown in Table B.2 
was used to control the evaporator for a 20% step change in feed flow. 
Hes resDOnsewO het heesys CenpiseppesencedminnbrcurenOo/-erltecan be 
observed from this figure by comparing with the response of controller 
FD0320 in Figure 6.2 that the main effect of the integral action has been 
to eliminate the offsets. However in this case the maximum deviations 
are slightly higher than when using proportional control only. This 
indicated that the level of integral action was satisfactory and that 


the oscillatory behaviour of the evaporator was probably due to too high 
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a proportional gain. 

Since the system was on the verge of instability when the prop- 
ortional controller FD0330 was implemented on the evaporator, a new set of 
proportional gains was selected for the proportional-intearal controller 
FD1315, given in Table B.2. They were chosen to be 50% of the gains used 
in controller FD0330 since it is common practice to use 50% of the gain 
that produces sustained oscillations. The response of the pilot plant using 
controller FD1315 is shown in Figure 6.8 in which it can be observed that the 
oscillatory behaviour of controller FD1320 (Figure 6.7) has been almost 
eliminated. In controller FD1315.1 the integral constants were increased 
by a factor of 1.5 while keeping the same proportional constants as in 
controller FD1315. The response of the system with this controller is 
shown in Figure 6.9. The effect of increasing the integral action 
deteriorates the control. The maximum deviations are higher and the 
system is more oscillatory than in Figure 6.8. 

From these results it can be concluded that a reasonable set 
of values for a proportional-integral cmtroller are given by proportional- 


integral controller FD1315 defined in Table B.2. 


6.3 EXPERIMENTAL TUNING OF THE DYNAMIC COMPENSATOR FOR THE FIFTH ORDER 
MODEL OBTAINED USING THE INVERSE AND DIRECT NYQUIST ARRAY METHODS 
The selection of the proportional gains for the dynamic 
compensator was done in the same way as for the static compensator in 
the previous section. However since this compensator was design 
significantly later in this study than the other controllers a new set 
of runs was done to evaluate and select the proportional gains for this 


dynamic compensator. All the experimental runs presented in this section 
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are from this*recentisetsofiruns. 

The basic dynamic compensators obtained for the fifth order 
evaporator model using the inverse and the direct Nyquist array methods 
were essentially the same as the one given by Equation (3.56). The 
Stability margins for which the closed-loop will remain stable when 
using this compensator are given by Equations (3.59), (3.60) and (3.61). 

Figures 6.10, 6.11, 6.12 and 6.13 present the responses of the 
evaporator for a 20% step change in the feed flow when the gains used in 
each loop are 10% (controller NADY0510). 20% (controller NADY0520), 30% 
(controller NADY0530) and 50% (controller NADY0550) of their correspond- 
ing ultimate gain, respectively. The gain used in the concen- 
tration loop is actually 80% of the values indicated above. The numerical 
values of these controllers are presented in Table B.3. There is a low 
level of noise in the response of the system in each case. The experi- 
mental results are in this case, in close agreement with the simulated 
results (see Figures 3.35, 3.36 and 3.37). The effect of increasing the 
gains in each loop is reflected mainly in the reduction of the offset in 
tne first effect level, Wl. However the response of the system also 
becomes more oscillatory by increasing the gains. The product concentra- 
cons Cys and the second effect level, Wo» are not affected significantly 
when the gains in each loop are increased. 

From these results it was concluded that a suitable set of 
gains will be the ones used in controllers NADY0520 or NADY0530. 

Figure 6.14 presents the response of the evaporator for a 10% 
step change in the setpoint of the product concentration using 
controller NADY0520. Although the dynamic compensator makes the system 


model almost non-interacting the experimental response shows that there 
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is still some interaction in the system as evidenced by the small 
perturbations in Wy and Wo. 

The use of the dynamic compensator with only proportional 
control action produced very good results. Therefore integral action 
was not applied in this case. 

During the second set of runs it was observed that there was 
less noise in the system and a slightly lower (i.e. 80%) effective gain 
had been used in the Cy measurement transducer. For these reasons and 
in order to make a more fair comparison among different controllers some 
additional runs were performed using the constant compensator discussed 
in the previous section and the optimal multivariable control discussed 
iieoecu10ns.0. 571; 

In the recent set of runs the pilot plant evaporator was also 
operated using controllers FD0520 and FD0530 which are based in the 
Static compensator given by Equation 6.1. These controllers, defined in 
Table B.1, are almost identical to controllers FD0320 and FDO330, 
respectively. The only difference is that the effective gain in the 
loop corresponding to the product concentration, Cos is 80% of the yalue 
of the gain used in the latter controllers. Figure 6.15 and 6.16 present 
the responses of the pilot plant for a 20% step change in the feed flow 
when the controllers FD0520 and FD0530 are used. The performance of the 
system was in this case much better than With controllers FD0320 and 
FD0330. This was mainly due to the lower level of noise in the system 
and lower gain in product concentration measurement system. 

As in the previous runs when controllers FD0320 (Figure 6.2) 
and FD0330 (Figure 6.3) were used the main effect of increasing the gains 


in the loops is to decrease the offset in We The response of the 
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system was more oscillatory when controller FD0530 was used which confirms 
that the gains used in controller FD0520 or in FD0320 are the most 
Suitable ones. 

In Figure 6.17 controller FD0520 has been used to control the 
evaporator during a 10% step change in the setpoint of C,. Since the 
controller is only diagonally dominant (as opposed to non-interacting) 
some interactions are present in the responses of W, and Wo as they were 
in Figure 6.5. 

An optimal multivariable controller, OP0501, given by Equation 
6.19 and discussed in Section 6.5.1 was also used in this set of runs to 
control the evaporator during a 20% step change in the feed flow 
disturbance. The experimental response of the evaporator is presented 


in Figure 6.18. 


6.4 EXPERIMENTAL RESULTS WITH THE OPTIMAL REGULATOR CONTROL METHOD AND 
THE THIRD CRDER EVAPORATOR MODEL 
The Optimal quadratic regulator control method is probably the 
most widely known of the so-called multivariable time domain methods. 
This technique has been studied extensively at the University of Alberta 
by Newell [53]. The method can be summarized in the following way: 


Given the discrete state space model 
x(n + 1) = Ax(n) + Bu(n) + Dd(n) (6.7) 
y(n) = Cx(n) | (6.8) 


and the discrete linear quadratic performance index: 
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N-1 
J= } (x(n) - ipl Q(x(n) - X 4) + u(n - 1)! Ru(n - 1) (6.9) 
n=] 


where: 
x4 is the desired state vector 
Q is a positive semi-definite symmetric matrix 
R is a symmetric positive definite matrix (but since this is a discrete 
formulation it can be made arbitrarily smal1) 
N is a large integer such that at n = N the process has reached 
steady state the recursive relations which define 
the control matrices have converged to constant 
values. 
The optimal control method involves minimization of the quadratic 


performance index and produces a discrete law of the form: 
u(n) = Kx(n) (6.10) 


The matrix K is defined by a recursive relationship [73] which converges 
to constant values after a few interactions. Note that under this scheme 
the control is exercised over the system states rather than the outputs 
as in the frequency domain methods. 

Tne design freedom in this method is in the selection of the 
weighting matrices Q and R. In this study an exhaustive search was not 
made to determine the best values for these matrices. The values used 
in these experiments are the best ones found in previous work by Newel] 
[53] and Wilson [77]. 

For the third order model the state yariables are the same as 


the output variables. Consequently there is no need to estimate any of 
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the state variables. 


6.4.1 Proportional Control 
The values of the Q and R matrices used to calculate the 
optimal proportional control for the third order evaporator model are 


given by: 
Q = diag[10,10,100] (6.11) 
R = diagl0.05,0.05,0.05] (6m 12) 


The proportional controller matrix for a sampling interyal of 


64 seconds is given by 


4.904 -0.401 noe 
KoP0300 =) |) 5ah84 -].600 4.425 (6.13) 
4.093 9.655 B35, 


Figure 6.19 snows the response of the evaporator, using 
controller OP0300, for a step change of 20% in the feed flow disturbance. 
The gains of this controller seem to be high. There are small oscilla- 
tions in all the controlled variables but the offsets are very small. The 


response of the same controller for a 20% step change in feed concentration is 


shown in Figure 6.20. The control is very good in all the output variables. 
This proportional controller generates a system which is 
diagonally dominant and it is probably because of this property that the 


system is non-interacting. 
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6.4.2 Proportional-Integral Control 


The state vector was augmented, as described by Newell [53], to 


generate integral action on each output variable. 


The R matrix used in 


this case is the same as Equation (6.12) while Q matrix is given by: 


Oe=ediagthos 10,100,151 51 | 


For a sampling interval of 64 seconds 


control; OP1300, is equal.to 


9.490 
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The response of the evaporator, using 


20% step change in feed flow is shown in Figure 


(6.14) 


the proportional-integral 


437 (6.16) 
.090 


controller OP1300, to a 
6221.0 LWhereftfectaor 


adding integral action can be appreciated by comparing it with the 


proportional control OP0300 in Figure 6.19. 


With integral action the 


small offsets are eliminated, the system is a little more oscillatory 


and has larger overshoots. 


Other work [27,53] has shown that optimal 


proportional and proportional-integral controllers can give much better 


control of the evaporator than shown in Figure 6.19 and 6.21. 
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6.5 EXPERIMENTAL RESULTS WITH THE OPTIMAL REGULATOR CONTROL METHOD AND 
THE FIFTH ORDER EVAPORATOR MODEL 

In the fifth order model the number of states is greater than 
the number of outputs. The implementation of the fifth order optimal 
control requires the estimation of the unmeasured first effect concentra- 
con's C,. In previous studies this estimation was done using a Kalman- 
Bucy filter estimator [2,25] or a Luenberger observer [34,74]. However 
in this study a procedure programmed by Newell [53] was used, because of 


its simplicity and convenience. 


6.5.1 Proportional Control 


The proportional controller OP0500, was calculated for a 
Sampling interval of 64 seconds and using the following matrices in the 


quadratic performance index: 


1) 
i 


dvag (sl Gadi 1105100) (6.17) 


ye) 
iT] 


diag{0:05,0.05,0.05 | (6.18) 


The calculated optimal control, OP0500, is equal to 


6.511 -1.226 6, Oey = Ga0C7 ome laedc 

P — -_ 

Sonu © || AB2 0.368 0.689 ice) 9.707 | (6.19) 
3.054 1.056 0.0178 9.75 11.570 


The response of the evaporator using proportional controller 
OP0500 for a 20% step change in feed flow, is shown in Figure 6.22. The 
control of the product concentration, Co and the second effect holdup, 


Wo » is very good. There is only a small offset in the first effect 
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holdup, W- (During the most recent set of runs this controller was also 
used to control the evaporator. The response of this controller (0P0501) 
was presented in Figure 6.18. The results in Figure 6.18 are very 
similar to the one presented in Figure 6.22 and they are typical of the 
better results obtained with optimal quadratic controllers by previous 
workers [53,77]. It can be noticed also that there is evidence of less 


noise in Figure 6.18.) 


6.5.2 Proportional-Integral Control 


For this case the state vector was also augmented as described 
by Newell [53] to generate integral action on each output variable. The 
R matrix used in this run was given by Equation 6.18 and the Q matrix was 


defined by: 
OeaeO tad Osco. 00s) 's ist (6.20) 


The calculated control matrices for a sampling interval of 64 


seconds are: 


8.209 ext: -3.640 0.139 -15.44 

p ms 

Pamir 4.542 | 0.371 0.553 1.284 9.07 | (6.21) 
4.238 1.166 -0.059 12.26 14.22 
e277 0.027 -1.432 

I . zt e 

ee orcad | Wolk 0.298 0.895 (6.22) 

0.648 1.936 1.306 


The response of the evaporator, using the proportional-integral 


controller OP1500, for a 20% step change in feed flow is shown in Figure 
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6.23. This is more oscillatory than the response of the system with only 
proportional control. However the offsets are eliminated by the integral 


action. 


6.6 EXPERIMENTAL RESULTS USING A MULTILOOP CONTROL SCHEME 

The multiloop case was implemented by using diagonal matrices 
in the multivariable proportional plus integral control algorithms. This 
makes it equivalent to single variable control and also provides a "base 
case" for comparison with multivariable techniques. When the control 
parameters are the same as the diagonal elements in the multivariable 
controllers then a comparison of the two cases gives an indication of 
the effect of the otf-diagonal elements in the controller matrices. 

The pairing of the input and output variables for the multiloop 
evaporator control scheme was obtained by analyzing the direct Nyquist 
array (see Chapter Five). It was found that the configuration which is 
the least interacting is given by the following loops: 

- product concentration, Co controlled by steam flow, S. 
- first effect holdup, Ws controlled by bottoms flow B,- 
- second effect holdup, Wo» controlled by bottoms flow B.. 


Similar results were obtained by Newell [53] using Bristol's 


approach [8,9]. 


An exhaustive search for the best proportional and proportional- 


integral constants was not made. Instead the selection of the 
controller constants was made mainly to illustrate the effect on the 
evaporator response of the off-diagonal elements in the controller 
matrices discussed in the preceeding sections (i.e. to illustrate the 


advantage of multivariable versus multiloop control). 
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6.6.1 Multiloop Proportional Control 


The proportional gains used in the multiloop proportional 
control are shown in Table 6.2 for the run MLO200. These constants are 
equal to the corresponding diagonal entries of the best control] matrix 
(FD0320) found using the third order evaporator model and the inverse or 
direct Nyquist array method. In this way a direct comparison can be made 
between multivariable proportional control FD0320 and multiloop propor- 
tional control MLO200. The response of the latter controller for a 20% 
step change in feed flow is shown in Figure 6.24. Comparing this 
response witn the response of proportional control FD0320 in Figure 6.2 
it can be observed that the multiloop system produces higher offsets and 
higher maximum deviations in the three output variables. Figure 6.25 
Shows the response of control MLO200 toa 20% step change in feed 
concentration. Since this disturbance affects only the product concen- 
tration, Cos its effect in the control system is very small. These responses 


under multiloop control are typical of the best results obtained in other 
studies [27]. 


6.6.2 Multiloop Proportional-Integral Control 


The proportional and integral constants used in the multiloop 
proportional-integral controller were equal to the corresponding entries 
of the best proportional-integral matrix (FD1315) found using the third 
order evaporator model and the inverse or direct Nyquist array method. 
The controller constants used in this controller are the ones shown in 
Table 6.1 for run ML1200. The response to a 20% step change in feed 
flow is shown in Figure 6.26. The addition of integral action to the 
proportional controller eliminates the offsets. The negative aspects 


of introducing integral action are that some of the maximum deviations 
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are slightly higher especially for the product concentration, C The 


9° 
comparison of the response of the multiloop proportional-integral control 
ML1210 of Figure 6.26 with the response of the multivariable proportional- 
integral control FD1315 of Figure 6.8 shows that the multivariable 


control reduces considerably the maximum deviations in the output 


variables and makes the system respond faster. 


TABEESG 32 
CONTROL SETTINGS FOR MULTILOOP RUNS 


Control Integral 
Interval Proportiona| Constant 
Run Loop (Sec. ) Constant (Min-!) 


6.7 COMPARISON AND CONCLUSIONS OF THE EXPERIMENTAL RESULTS 


6.7.1 Proportional Versus Proportional-Integral Control 


The addition of integral action in the multiloop control scheme 
improved the response of the control system dramatically. The severe 
offsets associated with the proportional controller were eliminated with 
the integral action. An analogous improvement was not observed with 
the multivariable control systems. The addition of integral action to 


the multivariable controllers produced no significant improvement in 
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the response of the evaporator. The response of the system using only 
proportional control was very good and the offsets were very small due to 
the relatively high proportional gains. Thus there was little room for 
improvement. For the controller designed using the inverse and direct 
Nyquist array the only noticeable effect of adding integral action was 
to reduce the offset in the first effect level. When the optimal control 
method was used the performance of the system with integral action was 
more oscillatory than when only proportional action was used. However 
this could probably be corrected by proper choice of the Q and R matrices. 
Thus, it was concluded that a good multivariable design method 
would produce a proportional controller that would give at least as good 
control of the evaporator as a multiloop proportional plus integral 


controller. 


6.7.2 About the Control System Design Method 


When applied to the third order evaporator model tne inverse 
Nyquist array [65], the characteristic locus [4,42] and the direct 
Nyquist array method produce essentially identical proportional or 
proportional plus integral controllers. These methods lead naturally to 
similar controllers whenever a constant controller matrix exists which 
makes the system almost diagonal. The response of the double-effect 
evaporator with these controllers is very good. The main advantage of 
these methods is that they produce a controller which is very easy to 
tune in the field. Once the basic controller and gain stability region 
-has been obtained, the final controller can be chosen by tuning of the basic 
controller in the field in a manner similar to the tuning of single- 


loop systems. 


senogegi at woterogave sad) 40 senoqzes-an 
ey . fc 


. F 
to oft? bee Doon vtoy caw Tonto; anoregegeags 


« 


- 


ten leno! soqors Cort vievtgslay! dg 
Ne 


yolioriocs oft vt) caoOmevoTgE 


a 
— 
4 
> 
ma 
— 
é 


ing ang ars sl 


; l +9a336 anh) ofc ni Psetto aa subi ® 


) a 
ant he a RATIO T7104 +t hseu 26W hone 


a 
»e 
te 


y ‘re iy fotiw WAI eos Li ffj20 S70 
| besne rs ad Videoo TY bh ug> ota 
gs 2209 Seba! sno" cmd 20H) -) OE 
) Tagotrwodenm & e2utord), bf 
hem »5 sodesedava’ sit 0 fora 
ottoata 
— 
F itz 2 (oxen say guodA a) 


¢ ent od [ebo s oueya yalro bride oft ib Seq ges age” 
! . . , 4 i ai 
sont ae. tng [22,9] pool aiJerwroo meta ning su Kory chad 


ie 
46 [pnotivoqgetg [saidnebr \! ohinviens peed bottom Vex%6 yak rt : 


du 


159 


For the fifth order model the inverse Nyquist array method was 
unable to produce a constant compensator which would make the CLTFM 
diagonally dominant. Thus it was not possible to design a proportional 
controller for the fifth order model with this method. On the other hand 
the direct Nyquist array method did produce a proportional controller 
for the same model which was almost the same as the one obtained for the 
third order model. 

When dynamic elements were used in the compensator the inverse 
and the direct Nyquist array methods both produced essentially the same 
dynamic controller for the fifth order model. The response of the 
double effect evaporator with this dynamic controller was excellent. 

Another important feature of the inverse Nyquist array [65], 
the characteristic locus [4,42] and the direct Nyquist array method is 
that the control system design produced by these methods is capable of 
handling set point changes as well as rejecting any type of unwanted 
disturbances. This is not the case with the optimal control method 
where different control matrices are required for "“optimal' setpoint 
control and to regulate the system against undesirable disturbances. 

The optimal control method also produced very good control. 
The use of the fifth order model gave much better results than using the 
third order model. The critical step in this method is the selection 
of the matrices Q and R of the quadratic performance index and unfor- 
tunately the method includes no direct means of estimating what Q and 
R should be to produce a specified type of response. Another disadvan- 
tage of the method is that the control is based on the state yariables 
instead of the outputs. In many cases this requires the estimation of 


the states which increases the complexity of the control system. 
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Of all the methods studied here the multiloop design approach 
produced the worst control system and one would expect the performance to 
be even worse for processes with larger interactions, However this scheme 
is the simplest and most familiar approach for designing control systems. 

The controller which gave the best experimental performance 
was the tuned dynamic compensator obtained with the direct and the 
inverse Nyquist array methods using the fifth order model. The second 
best performance was obtained using the optimal control method and the 
fifth order model. (This was a 3 x 5 controller matrix.) The third best 
performance was obtained using the tuned constant compensator designed 
with direct Nyquist array method and the fifth order model or the third 
order model or using the inverse Nyquist array and the characteristic 


locus method with the third order model. 


Oey oO eeboutMthesOndervof, the: Model 

It is a significant fact that the same constant compensator 
was obtained for the third and fifth order models when using the direct 
Nyquist array method. In both cases the compensator represented the 
best that could be obtained when a constant compensator was uSed. 
However it was not possible to obtain the same high degree of diagonal 
dominance (i.e. non-interaction) with the fifth order model as had been 
obtained when using the third order model. This is obviously because 
the third order model is less complicated (and hence less representative 
of the real evaporator) than the fifth order model. With the inverse 
Nyquist array method [65] it was not possible to design a constant 
compensator that would produce a diagonally dominant system when using 


the fifth order model. This may indicate that the direct Nyquist array 
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method is less sensitive to reductions in the model order or accuracy 
and also less restrictive. 

In general, witn the inverse Nyquist array [65] the character- 
istic locus [4,42] and the direct Nyquist array methods the difficulty 
of designing a controller increases with increasing order of the model. 
This is because it is more difficult to visualize the type of operation 
or the type of controller which is required to obtain the specified 
design objectives. It also appears that the complexity of the controller 
matrix is increased (i.e. more dynamic elements are required) when the 
order of the model is increased even though the dimensions of the 
controller matrix remain the same. 

The experimental work has shown that there is not a direct 
relationship between the complexity of the controller and the performance 
of the system. There is not mucn difference between the performance of 
the pilot plant evaporator using a constant controller and using a 
dynamic controller which requires two phase lead and one phase lag 
compensator. (The performance of the system, however, is slightly better 


when the dynamic controller is used.) 


Ov As4 Whe ETfect of Diagonal Dominance 


The design of a control system for the fifth order evaporator 
model using the direct Nyquist array method, as described in Chapter 
Five, provided an opportunity to investigate the effect of diagonal 
dominance (of the open-loop and return-difference transfer function 
matrices) on the performance of the control system. In this example 
three compensators, Kn (Equation (5.52)), Kp (Equation (5.54)) and Kp 


(Equation (5.66)), were obtained that made the open-loop and return- 
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difference matrices diagonally dominant in different degrees. The vary- 
ing degree of diagonal dominance can be observed in the Nyquist diagrams 
of the diagonal elements of the respective open-loop transfer function 
matrix with their Gershgorin bands, i.e. Figures 5.25, 5.27 and 5.34 
respectively. 

Experimental runs using controllers based on compensators Kp 
and Kp have been discussed previously (i.e. controllers FD0320, FD0520, 
NADY0520). The controller based on compensator Ka given by Equation (5.52) 
was also implemented experimentally. The gains used with this compensator 
were chosen such that the elements of the resulting controller, DNAO520, 
were approximately of the same magnitude as the elements of controllers 
FD0320 and NADY0520. In this way a fair comparison could be made among 
controllers DNA0520, FD0320, FDO0520 and NADY0520. Controller DNA0520 
is presented in Table B.4. The responses of the double-effect 
evaporator, when using controller DNA0520, for a 20% step change in 
feed flow, feed concentration and for a 10% increase in the setpoint of 
the product concentration are presented on Figures 6.27, 6.28 and 6.29, 
respectively. The effect of improving diagonal dominance can be analyzed 
by comparing these figures with Figures 6.2, 6.4 and 6.5, respectively 
which present the response of controller FD0320 based on compensator Kp. 
All these ene are from the same set of runs and hence can be 
compared directly. 

For a 20% change in feed flow the improvement in the response 
of the system by improving the diagonal dominance of the controller is 
very small. The use of controller FD0320 gives smaller offsets and 


lower overshoots than when controller DNAO520 is used. For a 20% 
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decrease in the feed concentration there is no difference in the behaviour 
of these two controllers. The same applies for a step changed in the 
setpoint of the product concentration. In all three cases the controller 
FD0320 is a little more oscillatory wnich indicates that slightly higher 
effective gains were used in this controller. 

The effect of improving diagonal dominance can also be analyzed 
by comparing the responses of controller FD0520, based on compensator 
Kes in Figures 6.7 and 6.9 against the response of the system in Figures 
6.14 and 6.17 when controller NADY0520, based on compensator Kp» is used. 
As in the previous case there is no significant difference in the 
behaviour of these controllers for a 10% step change in the setpoint of 
the product concentration. For a step change in the feed flow there 
was a marginal improvement in the response of the product concentration 
when controller NADYO520 is used. It was surprising to find out in these 
experimental runs that a marginal improvement was noticed for load dist- 
urbances and not when a setpoint change was introduced into the evapora- 
tor because one would expect diagonal dominance to be more important in 
setpoint changes. 

From tnese observations it can be concluded that small 
improvements in the degree of diagonal. dominance in the open-loop and 
return-difference matrices produced only a marginal improvement in the 
closed-loop performance of the system. This means that a non-interacting 
system cannot always be justified and that Nyquist array design techniques 
are valuable because they allow the designer to trade off controller 
complexity versus the degree of non-interaction produced. 

The same general conclusions were obtained in Chapter Five 


based on simulated results. 
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CHAPTER SEVEN 
CONTROL SYSTEM DESIGN FOR MULTIVARIABLE PROCESSES 
CONTAINING TIME DELAYS 


7.1 INTRODUCTION 

Pure time delays are very common in process systems but can 
only be handled by a few design techniques. The presence of'a time delay 
in a transfer function makes it an irrational function and in most cases, 
unless they are approximated by rational functions, the basic Nyquist 


stability theorem [54] has to be modified to deal with irrational 


transfer functions [13,71]. Desoer and Wu [504] give sufficient 
conditions for the stability of multi-input, multi-output, linear, time- 
invariant, feedback systems which are a generalization of previous works 
on system stability that originally began with Nyquist. However, the 
design of multivariable systems with non-rational G(s) is still "not 
fully clarified [509] ". However transfer functions with time delays have 
the peculiarity that the basic Nyquist stability theorem can be applied 
to them directly without requiring any modification. One of the main 
reasons for the popularity of the single-input, single-output frequency 
domain techniques was due to this fact. In multivariable systems with 
pure time delays the direct Nyquist array design method is more convenient 
than the inverse Nyquist array method because, in the direct method, 
time delays cause the Nyquist locus to "Spiral into the origin" whereas 
in the inverse method the corresponding loci spiral out to infinity. 
The direct Nyquist design approach was implemented using 
numerical methods to perform operations such as matrix multiplication, 


inversion, etc. For example Q(s) = G(s) K(S) is calculated numerically 


164 


ava BIT ANO 
: r | ae 

93223099 SIGALRAVI TAM RO Bare as May ay? ORTOP 
eyAj3? OMIT OAIMIATS 


int mom 


a eve zeaogte, a} NORD Yay. Oar eyefad agtt sw a 
aq oT _esuptanon’ apizob wo 6 a patbaed od “S . 
: 1b 1 ptt (snntaavet po Ff tatsa M0 et SO Tetenew Ra 
bond oft . enol) am td BY 4d be! Smt MOAGNE , S16 vad: 229 
nakbewe ftw [ped of SertPhom ed o7 ten ise} moroodd yO TSE 
snabatytue svig (PO jan ” sanaeu ams? “stn oes rater 
ant eek isfum , Suge refun to yt titdedlé a "0? il . 
3) | pottes tis conmep, meas tote area sasdbes? S060% 
| saqyaxoH .s2tueyA nd lw negee vyiteatany li vtTl f0pie ae 
(i Sano tbe f0G HE Ie Sane: a bc ie id Pain FOR 7 
eyed ceptok caty néiw eaottonut IP ehay vovah,? 100 wus ol 


bs tTaae ad me marromrt? ee Hi tesa $ptooyh otaad eft Tank a trek Uae 


otsm wd. 6 ‘ant onata sot iibom, Aen ara yupan, fugdt IM, Yitosnth 


r soloution ota to ag Fane, ais ne ms 
; ” - li : 


| - —— er: Pate i " iv mre co a ni a ae Me ne 


a -” ‘ ‘ : : - py! i 7 a = } S: | 
oly "OC sais 6 ists | 
, ‘4 : bes Le ‘ P he >) mi: a 

2B 4 - ied 


Sh 
‘ eat a) as 


eal 


ig 


4- 
a 


_ 
7 


yak 


ue 


165 


at each frequency rather than trying to multiply G(s) and K(s) analyti- 
cally to produce a single matrix with each element being a complex 
function. It is relatively easy to incorporate time delays into the 
elements of G(s) and/or K(s) and since the calculations are done 
numerically the subsequent numerical multiplication step, plotting, 
etc remains unchanged. Further work would be required to determine 
how pure time delays could be incorporated into the inverse Nyquist 
array and the characteristic locus methods. 

The following section presents a design example using the 
direct Nyquist array method ¥ design a system, for a binary distilla- 
tion column which contains time delays in each element of its transfer 


function matrix. The results agree with previous work by Wood and Berry [79]. 


7.2 DESIGN OF A CONTROL SYSTEM FOR A BINARY DISTILLATION COLUMN MODEL 
In this example the direct Nyquist array method is applied to 
the design of a control system for a binary distillation column model. 
The model used was obtained by Berry [6] from an eight tray, nine inch 
diameter, pilot-scale distillation column by pulse testing. The column 
has a total condenser and a basket type of reboiler and is used to 
separate a binary mixture of methanol and water. The distillation 
column is located in the Department of Chemical Engineering and has been 
described in detail by Berry [6] and Wood and Berry [79]. The distilla- 
tion column used in this example can be represented by the following 


dynamic model: 
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= G(s) (deeb) 


where the output variables are the overhead product composition, Xn» and 
the bottom product composition, Xp The input variables are the reflux 
flow rate, R, and the steam flow rate, S. The plant transfer function 
matrix, G(s) is presented in Table 7.1. The design is carried out in 

the s domain, rather than the w domain used in previous examples to 
illustrate that both cases appear essentially the same to the user of 
the interactive design program. 

As can be observed from Table 7.1 the distillation column 
contains a time delay term in each element of its transfer function 
matrix. The Nyquist array of the distillation column is presented in 
Figure 7.1, for a frequency range of 0.1 to 3.0 radians/min. From this 
figure it can be observed that: 

- the reflux flow rate, R, and the steam flow rate, S, have almost 


the same effect on the overhead product composition, Xn 


the bottom product composition, Kes is affected more strongly 
by the steam flow rate, S, than by the reflux flow rate, R. 

This indicates that the present configuration is the best one 
to control the distillation column under a multiloop scheme. This fact 
is emphasized from observing Figure 7.2 which presents the Nyquist 
diagrams of the diagonal element of the distillation column with their 
Gershgorin bands. This figure shows that the plant transfer function 


matrix, G(s), is already diagonally dominant. Thus there is no 


UABicEm/ aa) 


DISTILLATION COLUMN TRANSFER FUNCTION MATRIX 


12.8 e 5 18.9 eic° 

Wace te ie ey 
G(s) = 

ie Sd ao 

noe ey i) detouael) 


Typical steady state operating conditions of the distillation 


column [79]. 
Stream 
Overhead* 
Reflux** 
Bottoms* 
Feed 


Steam** 


* Qutput variables: 


** Tnput variables: 


Flow Composition 
lb/min wt. % methanol 

bails 96.0 

Heo 5 96.0 

127 8 

2545 46.5 

lei 


overhead composition, bottoms composition 
reflux flow, steam flow 
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need to reorder the yariables. 

From Figures 7.1] and 7.2 it is also apparent that the diagonal 
elements in column 2 will be unstable under negative feedback. In order 
to correct this situation it is necessary to change the sign in this 


column by post-multiplying the plant transfer function matrix, G(s), by 


K = a2) 


The plant matrix G(s) is replaced in the subsequent operations 
by the OLTFM Q, (s) = G(s) K- The Nyquist array of the OLTFM Q, (s) is 
presented in Figure 7.3 and the Nyquist diagrams of the diagonal elements 
of Q, (s) are shown in Figure 7.4. The latter figure is plotted with an 
assumed loop gain of unity and it is clear from the plot of the second 
diagonal element that F(s) is NOT diagonally dominant. Howeyer F(s) can 
be made diagonally dominant by proper choice of the loop gains. Thus it 
is possible to design, in this case, a multiloop control system using the 
direct Nyquist array method. At this point there are two alternatives in 
the design. The first one is to proceed with the design of the multiloop 
control system and the second is to design a pre-compensator to improve 
the diagonal dominance of the system. Both alternatives have been 


evaluated in the following sections. 


Design of a Multiloop Control System 


Since the control of the present distillation column using a 


multiloop control scheme has been investigated by Wood and Berry [79] and 


168 


“ ° ‘ yf - 5 
7 . é TT a + 
a Se “Fa y : i 7) hey 
- > 
A 
{ 
> ai?) 
‘ - 
J 4 Z| q A ; 5 Ly as wu - ‘ 
- 

‘ : 
= ‘ ~ ‘o> 5 . +4 
Lares 6 A ‘ « ' ad 

’ ii 

e i, 
- — —— 

y 

~ i? a , & f ‘ if 7 iin F Sry i 
j Tah 7 ' - PMP t tf Ne er ¥s 
1 
= ss 
= 
~~ ‘ ’ q ° 
: 


een i . * . 7a 
' Mi ; wary 
; | 
be 3 
j 4 j ‘ <a 
} = gw Viet 
J . ¢ , 
_ ' 
> * a 
* 1? |. t . if fa 
@ 
“ ' » RS 6 a A 
’ ; ) ’ Pai a 3 i { 
‘ a a 
‘ nr Sse f OO} 
; : bi wa td ern f 
= 7 ‘ ' 8 a ‘ \- a . wit heal 
’ 4 
- al oh } i AA 
6G ; t26u wy is ropes ‘ 


r ¥ ve. ’ = 7 * r . Tae 4 5 rand 
t 


baad onten wiasey2 fovinos yop hilum B.o8so, 4h Or press 07.9 
7 | T leds 


* 


7 


4 
7 . : : : 7 errr i ae 
t govisenvedls e156. vient intog aids NM PFI YET 


ale 
© ~ 
sh 9 


= 


went 
it Fy ‘beso 1 
aw? 


. 
i ve ia 
4 


169 


the best controller constants have been obtained by experimental tuning 
by Berry [6] no attempt was made to select the final controller constants 
in this example. This section is presented to illustrate the information 
which the direct Nyquist array method provides. 

- From Figure 7.4 which contains the Nyquist plots of the 
Gershgorin bands for the diagonal elements it can be obseryed that the 
distillation column will remain stable with its loops closed for the 


following combination of gains: 


Loop 1: ae fe og = 0-926 (7.3) 
Loop 2: CRs aE 0.213 (7.4) 


This “region can be expanded by using the exact Nyquist locus 
of each loop when it is open and the rest are closed and a trial and 
error procedure. However no major benefits are obtained from the 
expanded stablity region. 

With the information provided by Equations (7.3) and (7.4) 
it is possible to select safely the proportional constants for each 


loop. 


Design of a Pre-compensator to Improve the Diagonal Dominance in the 
System 


The transfer function matrix of the distillation column shown 
in Table 7.1 can be considered to have the preferred structure ‘for a 
system which contains time delays. This is because the greatest time 
delays are in the off-diagonal positions and hence the augmented system 


can easily be made diagonally dominant by using a dynamic compensator. 
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This type of transfer function matrix is not uncommon. Hydraulic models 
have been found [3] to have the same type of transfer function matrix. 
The use of a static compensator does not significantly improve 

the dominance in this system. However using a dynamic compensator the 
following alternatives are possible: 

a) make the OLTFM an upper or lower triangular matrix 

b) improve the diagonal dominance of the OLTFM 

c) make the OLTFM a diagonal matrix (i.e., non-interacting ) 


Only the last two alternatives have been investigated in tnis study. 


1) Improving Diagonal Dominance in the OLTFM 


From the transfer function matrix, G(s) in Table 7.1 it can be 
observed that the elements in the same row have time constants of the 
Same magnitude. Thus it is possible to improve diagonal dominance by 
using only time delays in the compensator. From the matrix Q, (s) and 
its Nyquist array in Figure 7.3 it was found that diagonal dominance 
could be greatly improved in Q, (s) if the following elementary operations 
are performed: (With the present implementation, the use of pure time 
delays in the controller, K(s), requires that all the dynamic components 
be introduced at one time rather than in a series of steps. This 
limitation could be overcome by reprogramming but would reduce computing 


efficiency. ) 


2S 


column 2=.column 2 - 1.3.x e — xX colum 1 (725) 


and 


As 


column 1 = column 1 - 0.45 xe‘ x column 2. (7.6) 
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These operations are performed by using the following 


compensator: 


1 on eens 


Ko (s) = (Zen) 
{ipl 


The Nyquist array of the new OLTFM, Qo(s) = Q, (s)K,(s) and the 
Nyquist diagrams of its diagonal elements with their Gershgorin bands 
are presented in Figures 7.5 and 7.6 respectively. Figure 7.6 shows 
that the system is now strongly diagonally dominant and that a controller 
gain can be chosen that will make the system closed loop stable. 


The final compensator used in this case is given by: 


Ky = Ky Ky (7.8) 
or 
etre teas 
Ke = (7.9) 
Ors en] 


Padé Approximations of the Time Delays 


It should be pointed out here that in this example there have 
not been any approximations of the time delays present in the plant or 
the controller. However in order to investigate the use of approxima- 
tions for pure time delays the controller matrix K. was also generated 
using a second order Pade approximation [10]. The Nyquist array and 
Nyquist diagrams of the diagonal elements with Gersngorin circles for 


the OLTFM Qo(s)s in which a second order Pade approximation was used, 
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are presented in Figure 7.7 and 7.8. It can be noticed by comparison of 
these figures with Figures 7.5 and 7.6 that the second order Pade 
approximations compensate for the process time delays very well at all 
frequencies but the difference becomes apparent at high frequencies. 
Hence the use of the Pade approximation will lead to a smaller stability 


region. 


2) Designing a Non-interacting System 
The OLTFM Q, (s) (which is equal to G(s) but with the signs 


changed in the second column) can easily be made diagonal (non-interact- 
ing) by performing simultaneously the following column operations on 
Q, (s): 


1 Or! Sata | 


column 2 = column 2 - 1.478 x (57, 4-4) x ae 


xecolomniian (7.10) 


and 


Apsara 


column |] = column 1 - 0.34 x GG gesane Xx amas 


x column 2 Fee lie) 


The factors used in these operations were obtained by an analysis of the 
elements of matrix G(s). These operations are accomplished by 


post-multiplying Q, (s) by the following matrix: 


Garis flew aa 
2a 0Se 1) . 
ve (7.12) 


[Az4set | -4As 
“0.3479 -954-7) 


] -1.478( 


The Nyquist diagrams of the diagonal elements of the OLTFM, 
a5 (s) = Q, (s)K; are shown in Figure 7.9. The Gershgorin circles in this 


case have a zero radius. The final controller used to make the open- 
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loop transfer function matrix diagonal is then given by: 


jh eee 
Ka = Ky (7.13) 
or 
eds @ 4 oe 
ee (7.14) 
Bis oe 
Soils rrtorsars pole au 


Stability Analysis 


Controller matrix Kp given by Equation (7.9) and controller 
es given by Equation (7.14) produce a system with similar properties. 
From the Nyquist diagrams of the diagonal elements of Q,(s) and a5 (s) in 
Figures 7.6 and 7.8 the stability region for which the closed-loop 
system will remain stable can be obtained for each controller. In both 
cases the distillation column will] remain stable for the following 


combination of gains: 


: ere 
Loop 1: 0<k, < 0.73 7 1-28 (7.15) 


Loop 2: 0< ky < gag = 0.27 (7.16) 


From Figures 7.6 and 7.8 the crossover frequency of each loop can also be 
obtained. Again, in both cases the crossover frequency obtained for 


each loop is the same. 


LOOpel: crossover frequency : 1.5 radians/min 


Loopy 23 crossover frequency : 0.55 radians/min 


Non-interacting control of the present distillation column has 
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been studied previously by Wood and Berry [79]. The non-interacting 

control system was in this case derived directly from the transfer 

function representation of the distillation column and is the same as 

the one obtained by Wood and Berry [79]. (The same result could have 

been obtained by progressively increasing the diagonal dominance of the 
Nyquist array of Q(s).) Since they obtained the final controller 

constants by experimental tuning no attempt was made here to specify 
different constants. The following controller constants reported by Wood and 
Berry [79] for non-interacting control (which are the “best" values found in 
their experimental evaluation of conventional single variable control) were 


used in this study: 


Loop Xy - R: 
Proportional constant K = 0.35 
Integral constant Ky = 0.040 min’ 
Loop Kp = ge 
Proportional constant Ky an 2 
Integral constant K, = 0.014 min”! 


The complete distillation column including the-effect of the 
feed flow and feed concentration disturbances as modelled by Berry [6] 
was simulated using IBM's Continuous System Modelling Program (CSMP). 
Figures 7.10, 7.11, and 7.12 present the simulated response of the 
distillation column for a step change of approximately 14% (+0.34 1b/min) 
in the feed flow, when a multiloop, a diagonally dominant and 3 non- 
interacting control system is used, respectively. As expected from the 


previous Nyquist analysis the control system that makes the system 


. fi aan 
a a 


Pn 
oft . [052 soved Otte Duck ve Cowal am betbuda a 


~é , r 5 Sry iv 4 
iis to ' i a ae 


nwt ot mov qlee bav iit 06> otis at 26W. 4 


4) 
62 on’ ot bre, gatos * wiaekh heels gag To oot stieesget no ie 
bruco 3 mz aft) COST gered bee boo Ye pantsrd6 ane’ 
¢ val opath odd oahesayaat yfavizes1porg, vd pentatde 
for ot aft bentegdd Yadd erate | (2)) To ers te L 


saw tomes Of patows Teseembroqxe’ ye 2a n6d ene 


+ (At r - = 74 - | Eg ile ; =. 
4 9 iO 


taf fovanes pritwotlol SAF etanienod 302 


FTITIGI SND 


asta faetnor ont soe asits NOR 10? reve 
1. 


: - 
atone Tsnobeagdes to Horisurnrs rssndmisgxe ny 


7 Ag 
¢ ‘Svbuste zing nts 


Sa, “ | 


} -§ ~ . good me - 


sqypdanoo Tana} s2eqoxt 


: ito 
a on * GRBF2NOS 4 [sngedi 
Ny 
2 + ph tena 
sf .0 mae snatenco [anotiogeyl - ys oe 


Eyhey P1000. * 44 


175 


diagonally dominant (controller with only time delays given by Equation 
(7.9) behaves almost the same as the non-interacting control (given by 
Equation (7.9)). The only difference is that the response of the non- 
interacting control system leads the response of the system with the 
"approximate" controller slightly. The response of these two controllers 
is much more damped then the response of the multiloop control which is 
oscillatory. The overshoots in the case of the multiloop control are 
also higher. On an overall basis the response of the multiloop control 
system can be considered to be worse than the non-interacting and the 


diagonally dominant systems. 


fom CONCLUSIONS 

The direct Nyquist array metnod has been found to be a very 
powerful technique that can handle systems which contain time delays in 
the process system, G(s) and/or the controller, K(s), without requiring 
any approximations. The distillation column example presented in this 
Chapter illustrates this case. 

From the examples presented in this chapter and Chapter Five 
it can also be concluded that whenever a controller matrix exists which 
is physically realizable and makes the system non-interacting, it can 
be designed using the direct Nyquist array method. A very important 
advantage of the direct Nyquist array method is that with this method 
several alternatives can be investigated by compromising between the 
complexity of the controller and non-interaction in the system. In 
many cases, like in the aeesaniow ian oui example, it is possible to 
design a less complex system with almost the same characteristics as a 


non-interacting system. 
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CHAPTER EIGHT 
CONCLUSIONS AND FURTHER RESEARCH 


8.1 CONTRIBUTIONS 
The main contributions of this work are: 

1) Development and extension of the direct Nyquist array design 
method. 

2) A comparison of the inverse Nyquist array, characteristic 
locus and the direct Nyquist array design methods. 

3) Experimental implementation and evaluation of controllers 
designed using the inverse Nyquist array, characteristic locus, 
direct Nyquist array, optimal quadratic and conventional design 
techniques . 

4) Development, documentation and demostration of a set of inter- 
active computer programs to implement the multivariable 


frequency domain design techniques. 


8.2 COMPARISON OF FREQUENCY DOMAIN AND STATE SPACE METHODS 
The main advantage of the frequency domain techniques over 
the methods based on a state variable formulation are: 
1) The use of frequency domain methods obviates the need of 
having to determine a parametric model since frequency response 
data, which adequately defines the system, can be used directly. 
2) They take into account the sometimes conflicting design 
requirements for: stability, non-interaction, integrity and 
overall performance of the system. 


3) They produce a controller matrix which is easier to tune "on line” 
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because it includes a diagonal matrix of loop gains. 

The control system involves only the output variables, i.e. 

it does not require measurement of the state variables. 

Many control applications run into difficulty due to system ~ 
parameters which change as a function of time and/or due to 
errors in estimating parameters which mean that the model on 
which the control system design is based is different than the 
actual process on which the control system is finally implemented. 
Multivariable frequency domain techniques can deal with this 
problem in some cases. For example, if the change in system 
parameters (or measurement errors or noise) is equivalent to a 
change in one or more of the feedback gains {k;,i=1,..m} then 
the range of actual k values can be checked against the 
Stability limits on kK. generated during the frequency domain 
design. For each loop there is a range of values of kK. TOG 
which loop 7 is stable, and since the stability analysis is 
based on the Gershgorin or Ostrowski bands it is known that all 
the other loops will also be stable at least for the range of 
tk. ,i=1,..m} for which F(z) remains diagonally dominant. It is 
obvious that changes on any column of G(z),which when combined 
do. not produce any change in the diagonal dominance (i.e. radius 
of Gershgorin circles),would have no effect on the stability 
analysis and hence would not result in unstable system 
performance (The use of a "band" instead of a line locus during 
the stability analysis would normally result in conservative 


gain estimates). 
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The above result is a quantitative test. More subjectively 
it is known that the off-diagonal elements of G(z) have only a 
small effect on the off-diagonal elements of a diagonally 
dominant Q(z) = G(z)K(z). Thus one would expect that, at least 
for small changes in the off-diagonal elements, small parameter 
changes in G(z) would have a smaller effect on the closed-loop 
performance of diagonally dominant systems than if the loops 
were closed around G(z) directly. 

Thus the multivariable frequency domain design techniques 
are expected to give conservative controllers that are less 
sensitive to parameter changes than controllers designed by 


other methods. 


8.3 COMPARISON OF FREQUENCY DOMAIN DESIGN METHODS 

The characteristic locus methed [4,42] is. slightly more: 
general than the inverse [65] and the direct Nyquist array methods. 
However the use of characteristic values to determine the stability of 
the system makes the characteristic locus method a less practical 
procedure than the latter techniques, The characteristic loci diagrams 
are a very uselful tool to analyze control systems but they do not give 
adequate guidance about the type of controller which must be used to 
achieve the desired control specifications. Thus, with the characteristic 
locus method the selection of the controller cannot be done in a systematic 
way aS with the inverse and the direct Nyquist array methods when designing 
a constant compensator. However this problem is at least partially overcome 


in the characteristic locus method, when designing a proportional-integral 
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control matrix, if the procedure recently proposed by Kouvaritakis [508] 
is used. 

The inverse [65] and the direct Nyquist array methods are very 
similar. Their relationship is similar to the one that exists between the 
Single-variable frequency design procedures when using the inverse versus 
the direct polar plots. These methods have been found easier to use and 
more convenient than the characteristic locus method [4,42]. However the 
information provided by the Nyquist array diagrams is less exact than 
the information supplied by the set of the characteristic loci diagrams. 

The inverse and the direct Nyquist array methods produce 
basically the same controller when a high degree of diagonal dominance 
in the open-loop transfer function matrix is required. 

In spite of the similarities between the direct and inverse 
Nyquist array methods this study has found the direct Nyquist array 
method to be more convenient than the inverse Nyquist array technique 


for the following reasons: 


1) It involves the use of the familiar "direct" Nyquist diagrams 
that represent the mathematical relationship between physical 
input/output variables and hence are more familiar and mean- 
ingful to most users than the Nyquist plots of the transfer 
functions of the inverse of the open-loop transfer function 
matrix, a7! fz), used in the inverse method. This means that the 
information provided by the direct Nyquist array is more 
physically meaningful to the user than the information supplied 
by the inverse Nyquist array. Furthermore the DNA method can 


be interpreted as a direct and intuitive extension of conventional 
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sinte-input single-output (SISO) design techniques. Thus it 
relates more directly to the previous experience of most users. 

2) The direct Nyquist array method requires only that the return- 
difference matrix, F(z), be diagonally dominant (which is a 
physically meaningful characteristic) as compared to the inverse 
Nyquist array method which requires that both the inverse of 
the closed-loop matrix, ro), and the inverse of the open-loop 
matrix, 7! (z), be diagonally dominant (a condition which is more 
difficult to interpret physically). This makes the DNA a more 
convenient technique because it requires only that the (-1,0) 
point not be enclosed by the Gershgorin circles as opposed to 
the INA method which requires the same condition for the (-1,0) 
point and the origin. However if the plant transfer function 
matrix is known algebraically the requirement that the origin 
should not be enclosed by the Gershgorin circles in the inverse 
Nyquist array method can be relaxed by performing some additional 
calculations as explained in Section 3.3 

3) The DNA method can handle systems with non-square transfer 
function matrices directly (However, note that Q(z) is always 
square. ). 

4) The DNA method works with the control matrix, K(z), at all 
times (as opposed to K7!(2)) so the user is always aware of 
the form, complexity and magnitude of the elements of K(z) 
and hence it is relatively easier to be sure that the controller 
is practical and physically realizable. 


5) The intuitive derivation of the DNA method based on the 
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extension of conventional SISO techniques shows that when 

the return-difference matrix is diagonally dominant the 
Gershgorin circles and the Ostrowski circles are simply a means 
of approximating, or defining the limits of interaction 

terms in the analytical expression for h. (Zz) that relates the 
input/output pair Us cys. This analytical expression for h.(z) 
can be used to derive other limiting or approximate solutions 
that in many problems will permit a more accurate and convenient 
design than is possible with the Gershgorin or Ostrowski circles. 
For example the Gershgorin circles for a system with a "triangular 
OLTFM" will, in general, have non-zero radii whereas it is 
obvious from physical considerations that there are no inter- 
actions which influence system stability. The use of the ug? 
circles", which take into account any zero elements in the 
OLTFM, Q(z), show that h. (2) is equal to q, (Zz) when Q(z) 

is triangular. 

A final check on the stability and performance of the closed- 
loop multiple-input multiple-output (MIMO) system can be made 
using the Nyquist plots of {h,(z), Veta oHini, he PGP 1S 
diagonally dominant then for open-loop stable systems it is 
necessary and sufficient that each of the Nyquist plots of 

h, (z) satisfy the conventional Nyquist stability criterion. 

For the general case, the sum of the clockwise encirclements 

of the critical point (-1,0) by the Nyquist locus of 

fiz B= Tee must be equal to aise where Po is the number 


F 
of unstable poles of Q(z). (See Appendix D). 
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The DNA method provides a general procedure for pairing 
input/output variables of MIMO systems since the direct 
Nyquist array defines the relationship between each pos- 
Sible pair of input and output variables. In most 
practical applications the columns of Q(z) stiould be 
interchanged so that the "linkage" or "transmission" 
between VF and y, (represented by q.;(Z)) is greater than 
between y; and the other input variables. This procedure 
takes into consideration not only the steady state be- 
haviour of the system but the dynamic behaviour as well. 
The DNA is computationally more efficient than the inverse 
Nyquist array method since the latter method always requires 
an extra operation: the inversion of the plant transfer 
function matrix. 

It is easier to accomodate pure time delays in the DNA 
method since the Nyquist plot spirals in to the origin 
rather than out to infinity as it does in the INA method. 


The inverse Nyquist array method has also been found to 


have some advantages over the direct Nyquist array method, the most 


important of which are: 


1) 


The Ostrowski circles in the inverse Nyquist array method 
become small when the gains in the other loops are increased 
[69] . This means that the region of uncertainty associated 
with the inverse of the transfer function hy (2) is reduced 
when high gains are used. Note that in the limit as the 
radius of the Ostrowski circles approaches zero the Q54(Z_ 


approaches h.(z) (69] . In the DNA method, h, (z) and nyo (z) 
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must be calculated from Equation 5.18 and plotted (where 

he 7a((Z) is the value of h. (z) when all the gains, Kk, assume 
their maximun permissible values). The difference between the 
two loci shows the maximun effect on h, (z) that can be produced 
by increasing the other gains {k,, ioUeeen aii eemaimore 
exact indication of the effect of increasing the loop gains 
than implied by the radius of the Ostrowski circles in the INA 
method. Thus in the DNA method more exact information can be 
obtained by at the price of extra computer calculations. 

2) According to Rosenbrock [69] it seems that there is a tendency, 
in practice, for the off-diagonal elements of the inverse of 
the plant transfer function matrix, Gua to be relatively 
smaller than those of G. In this study this tendency has not 
been noticed. Also it is not clear whether this difference in 
"sensitivity" is an advantage or disadvantage. 

3) As in the single-variable case the inverse Nyquist array is. 
more convenient to use when the transducer feedback matrix is 
a nonconstant matrix. However it has been pointed out earlier 
that in these cases, it is always possible (and desirable) to 
make the transducer feedback matrix part of the plant transfer 


function matrix. 


8.4 FURTHER RESEARCH 
Some natural extensions and recommendations for future 
investigations in this area are: 


1) Development of a systematic procedure to design a control system 
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for non-square plants using the direct Nyquist array method (cf. 
work by Kontakos [507]). The research should be oriented toward 
the design of a non-square pre-compensator which in combination 
with the non-square plant will produce a square open-loop transfer 
function matrix that retains the most important characteristics 

of the plant. 

Further research is required into the analytical expression for 
h.(z) in order to derive other limiting or approximate solutions 
which will permit a more accurate and convenient design than is 


possible with the Ostrowski circles. | 


Serious consideration must be given to the combination of the 
direct Nyquist array method and the sequential return-difference 
method. The use of the direct Nyquist array procedure in the first 
stage of the design will minimize the basic problem of the sequential 
return-difference method which is the order of loop closing. 
Similarly the use of the sequential return-difference in the final 
stage of the design provides the direct Nyquist array method with 
a systematic and exact single-loop design procedure. 
Further work would be required to determine how pure time delays 
could be incorporated into the general direct Nyquist array, 
inverse Nyquist array, and the characteristic locus methods. 
In order to increase the flexibility of the computer-aided design 
package and make the design of a control system an easier task the 
following programs should be incorporated into GEMSCOPE: 

a) A minimal realization algorithm to obtain the state- 


space model from a transfer function matrix. 
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b) 


A program to obtain the transfer function matrix from 

the experimental frequency response of a system. Although 
the availability of a transfer function matrix is not 
required with the frequency domain techniques, the 


knowledge of this TFM could be helpful in the design. 
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Nyquist Contour in the z-Plane 
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te 1. US. 789 


Bicume gs Inverse Nyquist Array.of G(w) 
iG=3| 0.1 to 10.0 |G~!(w)| K = 1| 


Note: the code used to summarize the 
conditions used to generate this 
figure are defined in appendix B. 


Figure 3.4 


Diagonal Elements of G7! (w) with their Gershgorin 
Bands 
|G=3| 0.1 to 10.0] G-'(w)| K= 1| 
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Figure o%0 Inverse Nyquist Array % he 
Ky! 


|G=3| 0.1 to 10.0] Q, = (3.25)| 


Figure 3.6 


Diagonal Elements of a3! (w) with their Gershgorin 
Bands 
|G=3| 0.1 to 10.0] Q7'(w)| Kj! = (3.25)| 
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Figure 3.8 


Diagonal Elements of a3! (w) with their Gershgorin 
Bands 
|G=3| 0.1 to 10.0] Q5!(w)| k5'kj! = (3.27)(3.28)| 
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Figure 3.9 Inverse Nyquist nna if 93 (w wy 
|G=3| 0.1 to 10.0] Q5'(w wile Ke a (3552) 
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Figure 3.10 Diagonal Elements of a4! (w) with their Gershgorin 
Bands 


[G=3| 0.1 to 10.0] Q5'(w)| Kg'K3!KS'KD? = (3.32)| 
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Figure 3.14 Inverse Nyquist An of G(w) 
|G=5| 0.1 to 1.0] G- (R= 1 
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Figure 3.17 


Diagonal Elements of a7! (w) with their Gershgorin 
Bands , AS. 
|G=5| 0.1 to 1.0] Qo t(w)| Ky! = (3.39) 
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Figure 3.18 Inverse Nyquist acl of Qo(w) | 
[G=5| 0.1 to 1.0] Q5!(w)| Rlk,71 = (3.40) (3.39) | 
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Figure 3.19 Diagonal Elements of Canc with their Gershgorin 
Bands 
|G=5| 0.1 to 1.0] Q5'(w)| K5!Ky! = (3.40)(3.39)| 
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Figure 3.20 Inverse Nyquist oct) of Q3(w) 
}G=5| 0.1 to 1.0] Q3'(w)| KBhKs SUEY a erie, 42)(3-40)(3.39)| 
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Diagonal Elements of Q3! (w) with their Gershgorin 
Bands 
|G=5| 0.1 to 1.0] Q3!(w)| KgIKSIKGT = (3.42) (3.40) (3.39) | 
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Bigure 324 Diagonal Elements of aC with their Gershgorin Bands 


|G=5| 0.-01-to 1.0] Qe'w)| Kel... Ky! = (3.44) ... (3.39)| 
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Figure 3.25 Diagonal Elements of a5 w) W4 jth their Gershgorin Bands 
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Figure 5.20 Nyquist Array of Q3(w) 
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Figure 5.23 Diagonal Elements of Q,(w) with their Gershgorin Bands 
G=5| 0.1 to 5.0] Qg(wi| KyKokgkq = (5.43)(5.44) (5.45) 
BAT): 


24 


Nyquist Array of 
|G=5| 0.1 to 5.0| 


( 
“Bet 


) 


26/7 


ae. 


3 : T - - 
Pon. i —— = : pen dissen a —- 
; ; 
j : ~~? -: 
—) ' * i 
| ‘ . ay io 
oa 5 mm) 
' 
7 
4 ~ a 
J j i / 
a / 
vi ' ; Mi 
é P 4 
- ; a 
" ‘ 
u - 
e ' 
. , 
* ; al \ 
/ 
/ 
/ / 


268 
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NOMENCLATURE 


Alphabetic 


A nxn state plant matrix 

B nxg forcing function matrix 

B, Double-effect evaporator first-effect bottoms flow 

B. Double-effect evaporator second-effect bottoms flow 

C mxn output matrix 

Cy Double-effect evaporator product concentration 

d Disturbance vector 

d. Radii of the Gershgorin circles 

D Disturbance coefficient matrix 

e. mx] vector whose yen element is unity and the remaining 


elements are zero 


e(z) mxl error vector 

if. Scalar return-difference quantity 

f.(z) jth eigenvalue of matrix F(z) 

f, (2) (i,j) element of matrix F(z) 

F(z) mxm return-difference matrix 

g.(z) jth eigenvalue or diagonal element of matrix G(z) 

g.,(Z) jth column of matrix G(z) 

9; ;(2) (i,j) element of matrix G(z) 

G(z) mx2 or mxm plant transfer function matrix 

h.(z) Transfer function between UstY. pair when the other loops 
are closed 

H(z) mxm transducer feedback matrix 

i mxm identity matrix 

J Linear quadratic performance index 

k.(z) th diagonal element or jin eigenvalue of matrix K(z) 
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jth row of matrix K(z) 


2xm pre-compensator matrix 

Number of input variables 

mxm post-compensator matrix 

Interaction term for loop i 

Number of the plant outputs variables 

Number of the plant state variables 

Encirclements of the origin by the Nyquist locus of ir! (z)| 


Encirclements of the origin by the Nyquist locus of 
the jth eigenvalue of ro! (2) 


Encirclements of the origin by the Nyquist locus of 
the jth diagonal element of R7! (2) 

Number of the right-half-plane zeros of the open-loop 
characteristic polynomial 

(1,3) element of matrix Q(z) 


jth eigenvalue of matrix a7! (2) 


(i,j) element of matrix a7! (2) 
Quadratic performance index constant coefficient matrix 
mxm open-loop transfer function matrix 


mx] vector of setpoints or inputs 


jth eigenvalue of matrix R(z) or yin element of vector r(z) 


Input vector at the jth 


ith eigenvalue of matrix ae 


step 


(i,j) element of matrix Rl (2) 


Distillation column reflux flow rate, or quadratic performance 
index constant coefficient matrix 


mxm closed-loop transfer function matrix 


Double-effect evaporator steam flow to the first effect or 
distillation column flow rate 
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jth eigenvalue of matrix T(z) or Q(z) 


mxm return-ratio matrix 

2x] plant input vector 

Row of matrix V(z) 

Inverse of matrix W(z) 

w-bilinear transform variable 

mx eh eigenvector of matrix Q(z) or G(z) 

Matrix containing the eigenvectors, w.(Z), of Q(z) 
Double-effect evaporator first-effect liquid level 
Double-effect evaporator second-effect liquid level 
nxl state vector 

Distillation column bottom product composition 
Distillation column overhead product composition 
mx] plant output vector 

jth element of vector y(z) 


z-transform operator 


Greek Letters 


Subscripts 


Cc 


e 


Ostrowski coefficient factor 

Resulting function from mapping contour D (Nyquist locus) 
Radii of the Gershgorin circles 

Summation 

Interaction coefficient 


Frequency 


Closed-loop matrix or compensator matrix 


Error vector 
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fi Return-difference matrix 

i Running index, position or number of an element 
13 Diagonal element 

ery) Position of an element in a matrix 

I Integral control matrix 

j Running index, position or number of an element 
fe) Open-loop matrix 

P Proportional control matrix 

t Return-ratio matrix 

u Input vector 


Superscripts 


“a 


Element of the inverse of a matrix 


j Step number in the design procedure 

I Integral control matrix 

P Proportional control matrix 

i Matrix transpose 

] Second controller design for the same system and using 


the same design method 


Abbreviations 

CLCP Closed-loop characteristic polynomial 
CLTFM Closed-loop transfer function matrix 
DNA Direct Nyquist array method 

INA Inverse Nyquist array method 

MIMO Multiple-input multiple-output system 
OLCP Open-1loop characteristic polynomial 


OLTFM Open-loop transfer function matrix 


Sed 


SISO Single-input single-output system 


TFM Transfer function matrix 
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APPENDIX A 
PILOT PLANT EVAPORATOR AND ITS MODELS 


A.1 THE EVAPORATOR 

The pilot plant evaporator used in this work is a double 
effect unit with the two effects operating in series. This pilot 
plant has been described in detail by previous workers [27,53,77]. 

The major pieces of process equipment are shown in the schematic 
diagram in Figure A.1. The control loops shown in Figure A.1 
represent the multiloop control scheme applied to the evaporator in 
previous studies [27,53,77]. 

The first effect has natural circulation through its 18 inch 
long, 3/4 inch 0.D. tubes. It is heated by process steam. The second 
effect is a long tube vertical unit which is run in its forced cir- 
culation mode. It has three, six foot long, one inch 0.D. tubes. 

It is operated at a lower pressure than the first effect and is 
heated by the vapour produced in the first effect. 

The evaporator is fully instrumented and can be controlled 
by either Foxboro electronic controllers or under Direct Digital 
Control (DDC) from an IBM 1800 Data Acquisition and Control Computer - 
operating under MPX. Multiloop DDC can be applied directly using the 
computer control package; and advanced control schemes by user written 
programs which utilize a set of system programs to interface between 


the user and system control programs. 


A.2 THE EVAPORATOR MODEL 


The complete development of the double-effect evaporator 
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model has been presented by Newell [53] who derived fifth and tenth 
order non-linear models. Based on these models Wilson [77] calculated 
discrete five-state and three state, linear, time-invariant models 
using a linearization procedure and Marshall's aatall reduction method 
oie 
The linearized models in the discrete form can be represented 

by: 

X(n+1) = ox(n) + Au(n) + ed(n) (A-1) 


and 


y(n) = € x(n) (A-2) 


The elements of the vectors x, u, d, y are defined as norm- 


alized perturbation variables: 


where Wass is the normal steady state value of Ws 
nes VECLOVS =X, Und andsVeroy Chesil tctheorcerediScCretesmodel 
are defined in Table A-1. The coefficient matrices of this discrete 
time model, with a 64 second sampling interval, are shown in Table A-2. 
For the third order discrete model the state vector x is 
given by: 


C (A-4) 


» It 
X = (Ws W 9) 


2? 
The elements of the state vector are defined in Table A-1. The vectors 
u, d and y are equal to the ones defined in Table A-1. The coefficient 


matrices for this model, with a 64 second sampling interval, are shown 


in Table A-3. 
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DESCRIPTION OF THE EVAPORATOR VARIABLES 


SUCHE NEWS 8 


x! = EWI, Cl, HI, W2, C2] 

Wl First effect holdup 

Cl First effect concentration 
Hl First effect enthalpy 

W2 Second effect holdup 


C2 Second effect concentration 


Control Vector, u 


ul = [S, Bl, B2] 
S Steam flow 
BI] First effect bottoms flow 
B2 Second effect bottoms flow 


Disturbance Vector, d 


qd! = [F, CF, HF] 
F Feed flow 
CF Feed concentration 


HF Feed enthalpy 


Qutput Vector, y 


y' = [W1, W2, C2] 


Normal Steady State Value 


sis) Nae 
4.59% glycol 
T3952) BtU/ tbs 
Ale oe De 
10.11% glycol 
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156.9 BIU/1b. 
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FIFTH ORDER 


-0.0119 
0.0116 
0.0116 

-0.0138 
0.0137 


TABLE A-2 


DISCRETE EVAPORATOR MODEL 


0 -0.0008 -0.0912 
0.9223 0.0871 
-0.0042 0.4376 
-0.0009  -0.1052 
0.0391 0.1048 
=01031 Aen 
0 0 
0 0 
0.0848  -0.0406 
=0.0432 55.0 
1E.0 0 0 
0 0 0 
0 0 0 


ite) 


oo — oO) 2 © 


(Te-eo45sec.) 

0 

0 

0 

0.0001 

0.9603 
Onl eZ 0 -0.0050 
-0.035] 0.0785 0.0049 
-0.0136 -0.0002 0.0662 
0.0012 0 -0.0058 
0.0019 0.0016 0.0058 
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0.0033 
-0.0219 


TABLE A-3 


THIRD ORDER DISCRETE EVAPORATOR MODEL 


(Te=" 64sec) 
0 -0.0326 
Oi 0 (bys 020375 
0.9602 0.0529 
0 -0.0135 120 
0 -0.0157 C= 10 


-0.0400 0.0219 0 


-0.0811 
0.0854 
-0.0442 


0 
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-0.0406 
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APPENDIX B 
CONTROLLER MATRICES 


This appendix contains an explanation of the code used for 
the figures presented in this thesis. It also includes the numerical 
values of the control matrices that were implemented on the double- 
effect evaporator pilot plant. These matrices were designed based 
on the third or fifth order evaporator model and using the following 
frequency domain techniques: inverse Nyquist array, characteristic 


locus and direct Nyquist array method. 
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TABLE B-1 


PROPORTIONAL CONTROL MATRICES DESIGNED FOR THE 
THIRD AND FIFTH ORDER EVAPORATOR MODEL USING THE 
INVERSE NYQUIST ARRAY, THE CHARACTERISTIC LOCUS AND 
THE DIRECT NYQUIST ARRAY METHODS 


Relationship 
to the UIti- 
mate Gain, a 


Controller 
Designation 


FDO310 
FDO0320 Oa KY 4.92 9.84 isle -3.68 0.0 -2.22 
-4.85 -9.84 -9.84 
-4,65 0.0 8.34 
FD0330 O23 ky Tests VERTIS 8.34 sO.de OnU -3.33 
-7.27 -14.76 -14.76 
EDOD20% Ore KY 4.92 9.84 
EDUS30z% 03 KY Jetscme | 14516 


* The effective gain used in the C,, measurement transducer was 0.16 kK 


*x The effective gain used in the C., measurement transducer was 0.24 Kae 


....continued 


TABLE B-1 (Continued) 


Relationship 
Controller|to the k 
DesignatiomUltimate Gain 


a : 


Proportional Controller 


-7.75 


0 lag 


-9.20 


0 -5 .56 


-12.12 -24.6 -24 .6 
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TABLE B-4 


PROPORTIONAL CONTROL MATRICES DESIGNED FOR THE FIFTH 
ORDER EVAPORATOR MODEL USING THE DIRECT NYQUIST 
ARRAY METHOD 


Controller Ratt k k k Proportional 
Designation C lee 1 2 3 Controller 
ain 
-0.600 0 LEO 
DNA0500 oe 0 10 -0.760 0 -0.4 
-1.520 -1.0 -0.8 
-2.16 0 530 
DNA0520 OeZek B20 ee Sell -2.73 0 -2.0 


-5.47 -9.52 -4.0 
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TABLE B-5 


EXPLANATION OF THE CODE USED FOR FIGURES? 


Code 


G=5 
G=DISC 


Omitoe.0.0 


K, Ky= (5.43) (5.44) 


FDO0320 


FD 


NADY 


DNA 
INA + CL + DNA 


Explanation 


Third order evaporator model has been used to 
produce the Nyquist array, characteristic locus 
or the control matrix. 


Fifth order evaporator model has been used. 


Distillation column model has been used. 


Frequency range used to produce the Nyquist 
diagrams (rad/sec. or rad/min. ) 


The Nyquist array, or diagonal elements, of 
the open-loop matrix Q, (w) is displayed. 


The Nyquist array, or diagonal elements, of 
matrix Q5! (w) is displayed. 


Control matrices Ky and Ky has been used. They 
are given by Equations (5.43) and (5.44) respec- 
tively. 


Control matrix designation (see Tables in this 
Appendix). 


Frequency Domain Techniques (Inverse Nyquist 
Array, Characteristic Locus, and Direct Nyquist 
Array Method) were used to design controller. 


Dynamic compensator obtained using the Nyquist 
array methods. 


Direct Nyquist array method. 


Inverse Nyquist array or characteristic locus or 
direct Nyquist array method has been used to produce 
the control matrix under consideration, i.e. all 
produce the same controller. 


(1) Note: vertical bars are used to delimit each element of the code. 
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Table B-5 (continued) 


Code 
SIM 5 


Exp 


K=Table B-1 


+20%F 


+10% C., 


+20% FC 
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Explanation 


Fifth order evaporator model has been used to 
Simulate the results. 


Experimental results. 


The numerical value of the control matrix is 
contained in Table B-1. 


Proportional control only. 
Proportional plus integral control. 
Multi-loop control. 


Optimal control. 


A 20% step change in the feed flow was used in 
the run under consideration. 


A 10% step change (up and down) in the set point 
of the product concentration was introduced in 
the run under consideration. 

A 20% step change in feed concentration was used 


as the disturbance variable in the run under 
consideration. 


Part of first set of experimental runs. 


Part of second set of experimental runs. 
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APPENDIX C 


LIMITING VALUES FOR THE ¢..(z) FACTORS 


al 
This appendix derives the basic relationships that lead to 
the Gershgorin and Ostrowski bands and shows that they can be regarded 
as limiting values, or bounds, on the location of the Nyquist locus 
of the "exact" transfer function, h(z). 
When the return-difference matrix, F(z), is diagonally 
dominant, the limiting value of the 65, (2) factors in Equation (5.18) 
is given by 


Eyagyal ss (COR) 


a 
A specific case of a third order system has been chosen to show the 
validity of Equation (C.1). The second order case is trivial and the 
same approach can be used to prove Equation (C.1) for higher order 
systems. 

Equation (5.17) which gives the open-loop relationship 
between the input of the augmented plant Us and output Yo for a third 


order system when the other loops are closed can be written in the 


following general way: 


ho(Z) = Go9(Z) + o49(Z) Gyo(Z) + 99(Z) 439 (z) (C.2) 


where: 


(uta Jets ; (ute 

T¥Ky 1947 (2)/\T¥k3 593 3(2) 1¥kq 7944 (2) 

by9(2) eee ager z) (mia ) (Gea) 
I (ea eae lakeaag ie) 
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Given these Equations the problem can be stated in the following way: 
Show that 
|oyo(z)| <1 iC. 
and 
Idgo(z)| < 1 (cr 


When the return-difference matrix is diagonally dominant, then 


[ky y 1 €1a9q(2)1 + lag, (2) [3 


a, = ] (G 
J Se us ithe es 

ere kal fq, (z)| ¥ 195(z)|} fas (c 
& be k3933(Z)| : 

Equation (C.3) can be written in the following ways 
B 

ee (Cy 

i [ ; (is (3313 )| 
US Te OMEN lee 


where the z's have been dropped for convenience and g is defined by 


a! (es ee )- (eu & 
TK 1444/7 A T*K 33933 asi | 


Taking absolute values Equation (C.9) gives: 
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JI] =15 O14 On ® (CG; 


where Uys do are given by Equations (C.7) and (C.8) respectively and 
© is a factor greater than one defined by this equation. It should 


also be noted that following inequalities are valid: 


Oa, > | Gee 


and 


Od, > ] (Ge 


Replacing Equation (C.13) into Equation (C.12) we obtain 
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Equation (C.17) can be simplified by multiplying a, by a, and eliminating 


the common factors, then 


iseaher B (C.18) 
i jel 42 | SeSiea  ea | 
+k 94] | 1#k33933 4k1197] 
or 
11931 3 “33823 | K 171904 | 
+k 11977 | (14k 3 3933 14k 4974 
Idol < 7 i ee (C.19) 
ols 931 “338 23 peau “ Ly 
#494 | 1#k33933 k1197] 
But Oa, < @, therefore 
| K11931 | kK 39993 eel Ky 1924 | 
If an oe aie 108 
ld40! < 5 G (GE20)) 
(oe 1193) baades, oo K1192) | 
TH yay] | | T3353 K11974 
OG 
i 
Idy0 < yy <a (Geis) 


which shows the validity of Equation (C.5). Following the same approach 


it can be proved that 


|o35 , < 1 (GR22)) 
It should be noticed that the Ostrowski factor to determine the 
radius of the Ostrowski circles is given by: 


Ostrowski factor = max (o> On ) (Ge23)} 
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Equations (C.21) and (C.22) indicate that the region of 
uncertainty about where h,(z) is located, can be narrowed slightly 
by using the appropriate factor Oy and Oe instead of the maximum of 
these values. 
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APPENDIX D 
STABILITY THEOREM 


D.1 BASIC EQUATION 

Let th. (z), i=1,..m} be the transfer function between input i 
and output 7 when loop 7 is open and all the others are closed as 
defined by Equation (5.18) and as illustrated in Figure 5.2. Also let 
FI4 (2) be the cofactor of element (j,i) of the return-difference matrix, 
F(z) 

Equation (5.18) was derived independently but follows directly 
from the results obtained by Rosenbrock [68]. Using the notation defined 


above h, (z) can be written: 


J 
hy(z) = q,,(2) + #4, ___ (0.1) 


If loop i in Figure 5.2 is closed using the controller kK. 


the closed-loop characteristic polynomial for loop i will be given by: 


E ji 
Petals als 
& j#i 
iene ksh; (z) Sele k.q..(z) ct = (Ds23) 
Eoaatz4) 
or, bringing to a common denominator 
m 
(Tt kiq.s(z))F (2) 2 k 5954 (Z)F (2) 
JF 
as kh. (z) = JF (D3) 
eZ) 
Equation (D.3) then will be equal to 
1 + k«h. (z) = eee) aie (D.4) 
Peas) 


i, a+ f 
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where |F(z)| is the determinant of matrix F(z). 


becomes 


When all the loops {i=1,...m} are considered Equation (D.4) 


: felt) 
(1 + kh, (z)) = ip Fc (D.5) 


ae 


Ze 
me Ie) 
i=1 


Belletrutti and MacFarlane [4] have shown that, for the 


systems of interest in this work, the stability of the closed-loop 


system can be determined by the encirclements of the origin by the 


mapping of the contour D by |F(z)|.(See Chapter Four). In the 


following discussion it will be shown that, when F(z) is diagonally 


dominant the net encirclements of the origin contributed by the term 


in square brackets is zero and hence that the closed-loop system 


stability can be determined from the Nyquist plots of fh.(z), 1= 1) eelithe 


Dees TAB LELTY STHEOREM 


1) 


Let the return-difference matrix, F(z), be diagonally dominant 
and assume that the plant G(z) has Po unstable poles. 

Let D be the contour shown in Figure 3.1 and described in Section 
Bile o 

Let |F(z)| map D into Te encircling the origin Ne times clockwise. 
Let the diagonal elements, {f;.(z), ole Or E¢Z) map Di into 


T.. encircling the origin n.. times clockwise. 


11 
Let kh, (z) map D into Thi encircling the critical point (-1,0) 
Nha times clockwise. 


Let FIS (7) map D into Wa encircling the origin Me times clockwise. 
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Then if F(z) is diagonally dominant it has been shown [68] 


that 


, 
ee tes) ae OHS: 


Consequently since the minors FIJ (2) are also diagonally dominant 


m 
Meg = 421 "ay (D.7) 
4i 


Thus the number of encirciements of the origin by the mapping of D 


by |F(z)|™2 


will be equal to Ny» where 


m 

es & ie bP ae 

i = (D.8) 

The number of encirclements of the origin by the mapping of 

eee 

D by TIF (z) will also be equal to n,- Consequently from the principle 
of the argument it follows that the number of encirclements of the 


jie 


origin by |F(z) will be cancelled by the number of encirclements 


of the origin by 5 FG) and hence the stability of the closed-loop 

will be pe ne by the number of encirclements of the critical point 
(-1,0) by the mapping of D by {k.h.(z), j=1,..m} which is equal to the 
number of encirclements of the origin by the mapping of D by |F(z)]|. 


Thus the closed-loop system will be stable if and only if 
ch aaa (D.9) 


where Po is the number of zeros of the OLCP outside the unit circle. 
If the plant is asymptotically stable Baae O and then the 


system will be stable if and only if none of the Nyquist diagrams of 
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ksh (z) encircles the critical point (-1,0). 


APPENDIX E 


DIAGONALIZATION OF THE DOUBLE-EFFECT PILOT PLANT EVAPORATOR 
( NON-INTERACTING CONTROL ) 


E.1 DIAGONALIZATION OF AN OPEN-LOOP TRANSFER FUNCTION MATRIX. 

Rosenbrock [65] has shown that when the plant transfer function 
matrix, G(z), is non-singular and the poles of G(z) and the zeros of 
det G(z) lie inside the unit circle a stable control matrix, K(z), 
can be found, by performing elementary column operations on G(z), that 
will make the OLTFM, Q(z) = G(z)K(z), diagonal. Furthermore the OLTFM 
will have all its poles and zeros inside the unit circle in the complex 
plane. 

Although the third and fifth order evaporator models each 
have two poles at the (-1,0) point these models can be considered to be 
Stable and by drawing a small semicircle at z = 1 the poles can be 
considered to be located inside the unit circle. The determinants of 
the transfer function matrices also have their zeros located inside the 
unit circle. Thus it is possible in both cases to find a stable 
controller matrix that will produce a stable open-loop transfer function 
matrix. 

For the third and fifth order evaporator model a non-interacting 
controller has been designed that will mathematically cancel the off- 
diagonal elements. The controller and the OLTFM obtained for both cases 


are presented in the next sections. 


E.2 DIAGONALIZATION OF THE THIRD ORDER MODEL IN TABLE 3.1 
The non-interacting controller, K(z), designed for the third 


order model is given by: 
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| -0.62658 0 1.0 
K(z) = | -0.74844 0 -0 .40148 Lea) 
-0 .99229 1.0 -1.77564 


and the resulting open-loop transfer function matrix, Q(z)=G(z)K(z) is 


equal to: 


0.0811 ; F 
" 0.04062 
q(z) = 0 eT EGuae 0 (E.2) 
; ‘ 0.07045 
Zee] 


The controller obtained using the direct Nyquist array method 


(Equation 5.32) was approximately the same as Equation (E.1). 


E.3 DIAGONALIZATION OF THE FIFTH ORDER EVAPORATOR MODEL IN TABLE 3.5 
For the fifth order evaporator model the designed controller, 


K(z), which makes the OLTFM diagonal is given by: 


b 


- 3.1357 2— 0 1.0 
K(z) = | - 1.0 0 - 0.14577 a (Ee3} 
ee 25 OMG ey - 0.64340 a 
where: 
ee t20.7020 hae 22 = 029210 ol 2 - 0.8561 
z - 0.4384 z - 0.8898 Sen 
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and the resulting open-loop transfer function matrix is equal to: 


0.1190 e 0 0 
Q(z) = 0 0.4058 f 0 (ERs 4)) 
0 0 0.02202 g 
where: 
2 (z - 0.899) 
(ae) Wiz 08898) 
nie 0.4058 
z-1 


= (z+0.7628) (z - 0.8998) 
opm (zera0e9602)a(ze=5089216)m( ze 0004384 
The non-interacting controller is slightly more complicated 
than the dynamic compensator obtained using the direct Nyquist array 


method (cf. Equation (5.65)). 
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